ФЭА / АИТ / Пособие "Математическое моделирование и оптимизация химико-технологических процессов с применением моделирующей программы CHEMCAD " 4
(автор - student, добавлено - 29-04-2014, 16:28)
Математическое моделирование аппаратов
Каждая математическая модель аппарата в СhemCade связана окном редактирования (одним или несколькими) параметров математической модели. Вызов окна ввода (редактирования) параметров выполняется в режиме моделирования следующими способами: – двойным щелчком на изображении аппарата; – командой контекстного меню Edit Unit Op Data (Редактировать параметры аппарата). – соответствующие команды меню Specificatons (Спецификации). Команда Specifications/Select UnitOps (Спецификации/Выбор оборудования) позволяет выбрать отдельные единицы оборудования. Их выбор выполняется аналогично выбору потоков. Команда Specifications/All UnitOps (Спецификации/Все оборудование) позволяет автоматически выбрать все оборудование технологической схемы. Окна для ввода параметров появляются на экране последовательно, в соответствии с ID номерами оборудования. Вид окна ввода параметров определяется типом оборудования и используемыми параметрами оборудования, которые заложены в его модулях расчета. Окно может содержать один и более разделов. Ниже представлены окна редактирования параметров некоторых, наиболее часто используемых математических моделей аппаратов в ChemCade.
В ChemCad'e представлены модули расчета теплообменников Heat exchanger (HTXR) с одним или двумя входными потоками. При одном входном потоке модуль служит как нагреватель или как охладитель. Если у теплообменника два входных потока, то доступны более сложные режимы. Рассмотрим ввод параметров для двустороннего теплообменника. Окно ввода параметров содержит три раздела (Рис. 23)
В разделе Specifications (Спецификации) представлены основные данные, такие, как перепады давления, температуры потоков и т.д. Значения Pressure Drops: (defaults=0) (Перепад давления (по умолчанию=0)) вводятся в поля Stream ID (ID Потока). Это положительная величина или 0. Calculation Modes: (Режимы расчета:) позволяют выбрать варианты расчета. В списке Backcalc Mode (for Autocalc): (Режим обр. счета (для автомат.расчета)) содержатся варианты для обратного счета. Он используется только для режима автоматического расчета. По умолчанию установлено 0 No back calculation (0 Не включен). Последующие четыре варианта позволяют пересчитать температуру и расход входных потоков. Список Utility Option: (Вспом. опция) используется только для двусторонних теплообменников и позволяет так определить параметры второго выходного потока, чтобы он соответствовал требуемой (рассчитанной) тепловой нагрузке аппарата. По умолчанию установлено 0 Utility Option Off (0 Вспом. опция выключ.) При выборе другого варианта необходимо задать либо выходную температуру, либо долю пара в каждом выходном потоке. В соответствующие поля Enter one specification only (Введите только один параметр) вводятся значения следующих параметров: Temperature Stream ID (Температура потока ID) – температура первого Vapor fraction stream ID (Доля пара потока ID) – доля пара первого (вто Subcooling stream ID (Недогрев потока ID) – температура недогрева первого (второго) выходного потока. Задается температура ниже точки кипения; Superheat stream ID (Перегрев потока ID) – температура перегрева первого (второго) выходного потока. Задается температура выше точки росы, Соответствующие поля Delta temperature specifications: (Значения допустимых температурных погрешностей:) используются для задания: Min delta temperature (Минимальный перепад температур) – Hot outlet-cold inlet и Hot inlet-cold outlet (Горячий выходной – холодный Stream ID – stream ID (Поток ID – поток ID) – разности температур между Stream ID – stream ID (Поток ID – поток ID) – разности температур между Heat Duty (Тепловая нагрузка) – значение тепловой нагрузки. Задается для теплообменника с одним или двумя входными потоками. Для теплообменника с Heat Transfer Coeff. and Area specification: (Значения коэффициента теплопередачи и поверхности теплообмена:) используются для задания значений: Heat Transfer Coeff. (U) (Коэффициент теплопередачи (U)) – коэффици Area/shell (Поверхн./кожух) – площади поверхности теплообмена. Раздел Misc. Settings (Необязательные параметры (настройки)) используется для задания информации по кожухам, трубам, трубному пространству и др. параметрам (Рис. 24).
Здесь же выводятся рассчитанные значения: тепловой нагрузки теплообменника, среднелогарифмической разности температур, поправочного коэффициента и др. В разделе Cost. Estimations (Стоимостные оценки) приводятся значения рассчитанных стоимостных оценок (Рис. 25).
Ректификация в ChemCad'e представлена набором модулей: приближенного расчета ректификации (Shortcut Column (SHOR)), строгого расчета ректификации и абсорбции (TOWR Distillation Column), строгого расчета ректификации и абсорбции нефтяных смесей (Tower plus), расчета ректификации с химической реакцией (SCDS Column) и др.
Модуль TOWR Distillation Column широко используется при моделировании технологических схем. С его помощью рассчитывается колонна любого типа –ректификационная, абсорбционная, абсорбционная с кубом и стриппер. Модуль имеет множество спецификаций, что делает его чрезвычайно гибким в использовании. Окно ввода параметров модуля TOWR имеет пять разделов. Раздел General (Общая конфигурация) содержит General Model Parameters (Общие параметры мат.модели) (Рис. 26Рис. 29).
В списке Condenser type (Тип конденсатора) выбирается тип конденсатора: – 0 Total or no condenser (0 Полный или нет конденсатора) (установлен по умолчанию); – 1 Partial condenser (Парциальный конденсатор); – 2 Total with water decant (Трехфазный полный конденсатор со сливом воды); – 3 Partial with water decant (Трехфазный парциальный конденсатор со сливом воды). Если конденсатор отсутствует, то указывается тип, равный 0. Для полного конденсатора дистиллят будет жидким. Для парциального конденсатора дистиллят будет паром, даже если есть жидкий продукт, выходящий из конденсатора, который будет определен как жидкий боковой продукт. Если выбран 2 или 3 тип конденсатора, то используется метод расчета коэффициентов фазового равновесия в системах с несмешивающимися с водой компонентами. Следующие поля используются для задания: Subcooled temp. (Температура недогрева) – значение температуры переохлаждения (задается только для полного конденсатора). Top pressure (Давление вверху) – давление в конденсаторе или вверху Cond pressure drop (Перепад давления в конденсаторе) – перепада давления в конденсаторе, задается как положительное число. Colm pressure drop (Перепад давл в колонне) – перепад давления в колонне, В поле No. of stages (Число тарелок) вводится значение числа тарелок, включая куб и конденсатор. Тарелки нумеруются сверху вниз. Если есть конденсатор, он рассматривается как 1-я тарелка. Если есть куб, он рассматривается как N-я тарелка. Максимальное число тарелок – 100, минимальное – 2. Feed stages: (Тарелки питания:) используется для задания номеров тарелок питания. Они вводятся сверху вниз, и их количество должно быть равно числу потоков питания, подводимых к колонне. Боковые нагреватели/охладители задаются как энтальпийные потоки (потоки с энергией, но без расхода компонентов); их положение задается так же, как и положение тарелок питания. Положительная энтальпия соответствует нагреванию, отрицательная – охлаждению. В разделе Specifications (Спецификации) задаются параметры, используемые при расчете Heat And Material Balance Specifications (Теплового и материального балансов) (Рис. 27).
В списке Condenser mode: (Выбор режима конденсатора:) выбирается режим работы конденсатора: – -1 No condenser spec (see tray spec) (Есть конденсатор, но нет задания ха – 0 No codenser (Нет конденсатора); – 1 Reflux ratio(R/D) (Флегмовое число); – 2 Condenser duty, negative (Тепловая нагрузка конденсатора, отрицательная); – 3 Distilate temperature (Температура конденсатора); – 4 Distilate mole flowrate (Общий мольный расход дистиллята); – 5 Distilate rate of one component (Мольный расход компонента в дистилляте); – 6 Distilate component mole fraction (Мольная доля компонента в дистилляте); – 7 Distilate comp fraction recovery (Отношение количества одного компонента в дистилляте к его количеству в потоках питания); – 8 Distilate total fraction recovery (Отношение количества потока питания – 9 Two components molar flowratio (Отношение мольных расходов 2-х – 10 Volume flow rate (crude flowrate) (Объемный расход дистиллята); – 11 Distilate mass flowrate (Массовый расход дистиллята); – 12 Weight fraction of one component (Массовая доля одного из компонентов в дистилляте); – 13 Reserved for TPLS module (Резервировать для модуля TPLS); – 14 Reserved for TPLS module (Резервировать для модуля TPLS); – 15 Molecular wеight (Молекулярная масса дистиллята); – 16 API gravity of the distilate (Удельная плотность дистиллята по стандартам Американского института нефти); – 17 Volume fraction of the component (Объемная доля одного компонента в – 18 Reflux mole flowrate (Мольная доля флегмового числа); – 19 Reflux mass flowrate (Массовая доля флегмового числа). По умолчанию установлен режим 0 No condenser, указывающий на отсутствие конденсатора (например, в абсорбере). Любой другой режим больше 0 показывает, что конденсатор есть, и его характеристики должны быть заданы в поле Specification (Значение). Величина тепловой нагрузки должна быть отрицательным действительным числом. Для режима 0 ничего указывать не нужно. Для режимов 5, 6, 7, 9,12, 17 нужно задать Component (Компонент), а для режима 9 задать еще Divided (Делитель). Для режима -1 нужно задать Stage (Номер тарелки), в списке Select optional tray mode (Выбор режима тарелки) выбрать нужный режим с 0 по 28, указать Specification (Спецификации) тарелки. Для режимов 3, 4, 6, 7,10 нужно задать Component. В списке Select reboiler mode: (Выбор режима куба:) выбирается режим работы кубового кипятильника: – -1 No reboiler spec (see tray spec) (Есть куб, но нет задания характеристик); – 0 No reboiler (Нет куба); – 1 Vapor/Bottoms mole flow ratio (V/B) (Кубовое число); – 2 Reboiler duty, positive (Тепловая нагрузка куба); – 3 Bottom product temperature (Температура куба); – 4 Bottom mole flowrate (Общий мольный расход кубового продукта); – 5 Bottom mole rate of one component (Мольный расход одного компонента кубового продукта); – 6 Bottom component mole fraction (Суммарная мольная доля компонен – 7 Bottom comp fraction recovery (Отношение количества одного компонента в кубе к его количеству в потоках питания); – 8 Bottom total fraction, recovery (Отношение расхода потоков питания к – 9 Two components molar flow ratio (Отношение мольных расходов двух – 10 Volume flowrate (crude flowrate) (Объемный расход кубового продукта); – 11 Botton mass flowrate (Массовый расход кубового продукта); – 12 Weight fraction of one component (Массовая доля одного компонента в – 15 Molecular weight of bottom product (Молекулярная масса кубового продукта); – 16 API gravity of bottoms (Удельная плотность кубового остатка по стандартам Американского института нефти); – 17 Volume fraction of one component (Объемная доля одного компонента Режим 0 показывает, что куба нет (например, в абсорбере). Любой другой режим больше 0 показывает, что куб есть, и его характеристики должны быть заданы в поле Specification (Значение). Для режима -1 задаются спецификации аналогично режиму -1 конденсатора. Если в колонне присутствует отбор бокового продукта, то в Side Product Specifications (Спецификации бокового продукта) для каждого отбираемого потока задаются Stage (Тарелка) – номер тарелки отбора, Side product mode (Режим бокового отбора) и Specification (Спецификация). Раздел Convergence (Сходимость) содержит параметры Estimates and Convergence Parameters (Параметры оценки и сходимости) (Рис. 28).
General estimates (Общие оценки) содержат следующие параметры: Dist rate (Расход дистиллята) – оценку расхода дистиллята. Ее рекомендуется задавать в мольных единицах. Reflux rate (Расход флегмы) – оценку расхода флегмы. Оценочное значение расхода флегмы также вводится в мольных единицах. Если расход дистиллята и флегмы не известен, то эти поля остаются пустыми, т.к. модуль TOWR задает собственные оценки. Temperature estimates (Оценки температуры) используются для задания оценок: – Т top (T наверху) – температуры в верхней части колонны. Задается оценочная температура конденсатора или верхней тарелки. Если значение не – Т bottom (T внизу) – температуры в нижней части колонны. Задается оценочная температура куба или последней тарелки. Если температура не известна, то поле остается пустым; – Т2 – температуры на тарелке 2. Задается, если существует большая разность температур между конденсатором и верхней тарелкой. В соответствующих полях Side product estimations (Оценки бокового продукта) задаются оценки мольных расходов продуктов отбора. Если они не известны, то поля остаются пустыми. При задании Convergence (Сходимости) указывается: Iterations (Итерации) – максимальное число итераций; Initial flag (Признак начального профиля) – начальные оценки для расчета ректификации: – 0 None – установлено по умолчанию. Для этого требуется минимум – 1 Reload column profile (Перезагрузка профиля колонны) – перезагружает профиль колонны, вычисленный в предыдущем расчете (базовом варианте). Число тарелок и компонентов должно оставаться тем же, что и в базовом варианте; – 2 Temp, profile only (Оценка температурного профиля); – 3 Temp.&vapor profile (Оценка температуры и расхода пара); – 4 Temp., vapor, liquid profile (Оценка температуры, пара, жидкости); – 5 Temp., vapor, liquid pres (Оценка температуры, пара, жидкости и профиля давления). Режимы 2, 3, 4, 5 используются, если TOWR не сходится из-за чрезвычайно нелинейного характера температурного профиля пара или жидкости в колонне. Следующие поля используются для задания: – Tolerance (Погрешность) – допустимой погрешности; – Damping factor (Демпфирующий фактор) – демпфирующего отношения – Water Quench (Подача воды) – признака подачи воды. Если сверху колонны подается большое количество воды, и она присутствует на каждой тарелке, то сходимость может быть достигнута, когда признак установлен. Calculated results (Вычисленные величины) содержат значения (Рис. 29).
Condenser duty (Нагрузка конденсатора) – тепловой нагрузки конденсатора (отрицательная величина); Rebolier duty (Нагрузка куба) – тепловой нагрузки куба (положительная Reflux mole (Флегмовое число в мольных единицах) – флегмового числа Reflux mass (Флегмовое число в массовых единицах) – флегмового числа Reflux ratio (Флегмовое число) – флегмового числа. В разделах Cost Estimation (Оценка стоимости) 1 и Cost Estimation 2 представлены опции для оценки стоимости ряда конструктивных элементов тарельчатых и насадочных колонн, конденсатора, кипятильника (Рис. 30).
ChemCad предоставляет средства для решения большого количества задач, связанных с реакторами, начиная с простых стехиометрических реакций и кончая множественными кинетическими реакциями.
Модуль Stoichiometric reactor (REAC) моделирует стехиометрический реактор при наличии набора стехиометрических коэффициентов, ключевых компонентов и степеней превращения. Реактор может быть адиабатическим, изотермическим или с подводом/отводом тепла. Окно ввода данных параметров модуля содержит два раздела (Рис. 31). В разделе General Specifications (Общие спецификации) представлены опции для задания общих параметров и условий.
В области Specify Thermal Mode: (Задание теплового режима:) выбирается один из тепловых режимов работы реактора: 1. Adiabatic (Адиабатический) – адиабатический; 2. Isothermal (Изотермический) – изотермический; 3. Heat Duty (Тепловая нагрузка) – с заданной тепловой нагрузкой. В списке Key component (Ключевой компонент) определяется позиция ключевого компонента. Предполагается, что ключевой компонент является реагентом. Это обязательный для ввода параметр. В следующих полях задаются значения: Frac. Conversion (Степень превращения) – степени превращения ключевого компонента (значение от 0 до 1). Это обязательный для ввода параметр. Heat of Reaction (Теплота реакции) – теплота реакции. Она задается при Reactor Pressure (Давление в реакторе) – давление в реакторе. Если задано значение 0, то давление в реакторе будет равно давлению входного потока. Calc H of Reac. (Расч. значение теплоты реакции). В полях Stochiometric Coefficients: (Стехиометрические коэффициенты:) задается последовательный набор стехиометрических коэффициентов: отрицательных – для реагентов, положительных – для продуктов и нулевых – для веществ, не участвующих в реакции. При необходимости задание стехиометрических коэффициентов можно продолжить, выбрав следующих раздел-вкладку More Components (Другие компоненты) (Рис.32).
У модуля REAC может быть один вход и до 3-х выходов. Если есть несколько выходов, то первый содержит пар (если таковой присутствует), а второй и третий содержат жидкости (если таковые присутствуют).
Модуль Equilibrium reactor (EREA) используется для расчета равновесных реакторов. Он может работать с двухфазными системами, но реакция будет происходить только в одной фазе, которая задается пользователем. Окно модуля имеет два раздела. Раздел General Specifications (Общие спецификации) содержит параметры общего назначения (Рис.33).
В поле Number of reactions (Число реакций) задается число реакций. Можно одновременно задать до 20 реакций. В Pressure drop (Перепад давления) задается перепад давления в реакторе. Область Reactor Model (Модель реактора) содержит опции для определения модели реактора и фазы реакции. Модель реактора выбирается в списке Specify reactor type: (Типы моделей реактора:): – General equilibrium reactor (Общий равновесный реактор); – Shift reactor (Конвертор СО); – Methanation reactor (Метанатор). EREA может использоваться для моделирования любого набора реакций. Для общего равновесного реактора данные по равновесию и стехиометрии задаются в меню данных реакций. Для реакции конверсии СО и для реакционной системы получения метана в ChemCad'e имеются требуемые равновесные данные. Предполагается, что реакция происходит в одной фазе, которая задается пользователем: – Liquid only (Только жидкая) – фаза жидкая (по умолчанию); – Vapor only (Только пар) – фаза паровая; – Liquid reaction, Mixed phase (Жидкая реакция, смешанная фаза) – реакция протекает в жидкой фазе, смешанная фаза; – Vapor reaction, Mixed phase (Паровая реакция, смешанная фаза) – реакция протекает в паровой фазе, смешанная фаза. Любой равновесный реактор может быть адиабатическим, адиабатическим с заданной тепловой нагрузкой или изотермическим. Выбор теплового режима реактора выполняется в области Thermal mode: (Тепловой режим:): – Adiabatic (no heat exchange) (Адиабатический (без изменения тепла)) – – Isothermal (specify temp) (Изотермический (с заданной темпер.)) – для – Specify heat duty (Заданная тепловая нагрузка) – это адиабатический Задание режимов расчета выполняется в области Specify calculation mode: (Задание режима расчета:): – Reaction conversion (Конверсия реакции) – при задании конверсии реакции расчет равновесия не обязателен, а тепловые и материальные балансы Approach delta T (Температурное приближение) – задается температурное приближение в поле Temperature delta в зависимости от выбранного Approach Fraction (Степень приближения) – степень приближения к равновесию используется вместе с равновесной конверсией. Степени приближения к равновесию задаются в экранах реакций, которые выводятся на экран после определения числа редактируемых реакций. Окна содержат опции Раздел More Specifications (Другие спецификации) содержит дополнительные данные для ввода (Рис.34):
Number of iterations (Число итераций) – задается допустимое число итераций для обеспечения сходимости уравнений. По умолчанию установлено Tolerance (Погрешность) – задается погрешность для расчета равновесия. Edit reaction number (Редактирование номера реакции) – указывается, 0 – все реакции; -1 – ни одной реакции; N – реакция с номером N. Для выбора единиц измерения используются соответствующие списки области Reaction Engeneering Units (Единицы измерения параметров реакции): Temperature Units: (Единицы температуры:), Pressure Units: (Единицы давления:), Heat of Reaction Units: (Единицы теплоты реакции:), Molar Flow Units: (Мольные единицы расхода:). В области Convergence Method (Метод сходимости) выбирается метод сходимости. По умолчанию установлен метод, рассматривающий потоки компонентов как независимые переменные и, таким образом, отвечающий за обратимость реакции при указании приближенной доли. Другая опция позволяет применить метод, использующий превращение компонентов потоков питания в качестве самостоятельных переменных и, таким образом, не отвечающий за обратимость при указании приближенной доли. В поле Temperature reference for heat of reaction: (Контрольная температура теплоты реакции:) задается значение контрольной температуры реакции. В поле Calculated overall heat of reaction Rxn (Рассчитанная общая теплота реакции) выводится рассчитанное значение общей теплоты реакции. У модуля EREA один вход и может быть до трех выходов. Если есть два или три выходных потока, то первый выходной поток – пар (если есть), а второй и третий выходной поток – жидкость (если есть).
Модуль Kinetic Reactor (KREA) служит для поверочного и проектного расчетов кинетических реакторов идеального вытеснения RFR (РИВ) и реакторов идеального смешения CSTR (РИС). Каждый из реакторов (РИВ или РИС) может быть жидкофазным или газофазным. Допускаются также двухфазные реакторы, но реакция может иметь место только в одной фазе. Окно модуля имеет два раздела. Раздел General Specifications (Общие спецификации) содержит параметры, которые используются как для реактора РИВ, так и для реактора РИС (Рис.35).
В поле Number of reactions (Количество реакций) вводится число реакций. Допускается до 20 одновременных реакций. В поле Reactor pressure (Давление в реакторе) задается давление на входе в реактор. Если оно не задано, то используется давление во входном потоке. В Pressure Drop (Перепад давления) определяется давление, при котором протекает химическая реакция в реакторе. Для РИС эта величина равна нулю. Для РИВ профиль давления равномерно распределяется между входом и выходом по числу шагов интегрирования. В области Reactor Model (Модель реактора) определяется тип реактора и фаза протекания реакции. Тип реактора выбирается в списке Specify reactor type: (Определение типа реактора:): Cont Stir Tank Reac (CSTR) – реактор идеального смешения (РИС). Plug Flow Reactor (PFR) – реактор идеального вытеснения (РИВ). Следующие опции позволяют задать фазу протекания реакции: Liquid (Только жидкость) – (по умолчанию); Vapor only (Только пар); Liquid reaction, Mixed phase – реакция протекает в жидкой фазе, смешанная фаза; Vapor reaction, Mixed phase – реакция протекает в паровой фазе, смешанная фаза. В области Thermal mode: (Тепловой режим:) выбирается тепловой режим (вариант расчета): – Isothermal (specify temp) (Изотермический (определить температуру)). Для него рассчитывается тепловая нагрузка. В этом случае в поле напротив необходимо ввести значение температуры. В противном случае используется температура входного потока. Adiabatic (no heat exchange) (Адиабатический (без изменения тепла)).– Specify heat duty (Заданная тепловая нагрузка). Для этого режима в случае РИС с заданной тепловой нагрузкой итерационно определяется температура. В поле напротив задается количество тепла, добавляемое или отнимаемое от реакции. Используются глобальные единицы измерения. В случае РИВ тепловая нагрузка равномерно распределяется на каждом шаге интегрирования и затем используется для расчета температурного профиля. Spec PFR temp, profile (later) (Заданный температурный профиль для Specify PFR utility U (Заданное условие вспом. потока для РИВ). В этом В области Specify calculation mode (Определение режима расчета) выбирается режим расчета: Specify Volume, Calculate conversion (Задан объем, рассчитать степень Specify conversion, Calculate volume (Задана степень превращения, рассчитать объем). Если выбран режим расчета Specify Volume, Calculate conversion, то в поле Reactor Volume (Объем реактора) для РИС задается объем. Для РИВ вводится либо объем, либо объем как произведение длины, диаметра труб и числа труб. В режиме Specify PFR utility U объем задается только в виде произведения. При режиме Specify Volume; Calculate conversion объем реактора используется для вычисления степени превращения реагентов, а также выходных составов и условий. Если выбран режим расчета Specify conversion, Calculate volume, то объем используется в РИВ для определения размера реактора. В списке Key Component (Ключевой компонент) выбирается номер ключевого компонента, относительно которого определяется объем реактора. Расчет связан со степенью превращения одного (ключевого) компонента независимо от числа реакций. Если выбран режим расчета Specify conversion, Calculate volume, то в поле Conversion (Степень превращения) задается степень превращения ключевого компонента. При задании этого параметра определяется объем реактора. Для РИС объем определяется временем пребывания, для РИВ пока не будет достигнута требуемая степень превращения ключевого компонента. Раздел More Specifications (Другие условия) содержит геометрические параметры и численные параметры для выполнения расчета (Рис.36).
Топология модуля KREA зависит от наличия вспомогательного потока. Если этот поток отсутствует (термические режимы 1-4), то у реактора может быть множество входов и три выхода, где 1=пар, 2=первая жидкость, 3= вторая жидкость. Если задан вспомогательный поток (термический режим=5), то у реактора может быть два входа и два выхода. Первые вход и выход содержат потоки процесса, а вторые вход и выход – вспомогательные потоки. После ввода данных разделах модуля KREA и нажатия на кнопку OK появляется следующее окно ввода параметров химических реакций (Рис.37):
Похожие статьи:
|
|