- •Міністерство освіти і науки україни
- •Скопійовано з оригіналу-макета, наданого автором
- •1. Математичне моделювання хіміко-технологічних процесiв
- •1.1. Загальні поняття
- •1.2. Класифікація моделей.
- •1.3. Фізичне моделювання
- •1. 4. Математичне моделювання
- •Знак ( –) при коефіцієнтах порівнянь означає, що потік спрямований з
- •1. 5. Засоби складання математичних моделей.
- •1. 5. 1. Емпiричний засіб
- •Загальна оцінка експериментальних засобів.
- •Приклад
- •1.5.2. Експериментально - аналітичний засіб
- •1. 5. 3. Теоретичний засіб
- •1. 5. 4. Зіставлення засобів побудови математичних моделей
- •5. 5. Вірогідність та простота моделі
- •5. 6. Рішення порівнянь математичного опису
- •5. 7. Перевірка адекватності та iдентифікація моделі
- •5. 8. Вибір математичної моделі
- •2.Моделювання Хімічних Реакторів
- •2. 1. Модель реактора ідеального змішування
- •2. 1. 1 Модель різ для опису стаціонарного режиму
- •2. 1. 2. Модель різ при протечі деяких реакцій
- •2. 1. 3. Дослідження моделі різ
- •2. 1. 4. Побудова q - t -діаграми і дослідження стійкості стаціонарних режимів різ
- •2. 1. 5 Умова стійкостi
- •2. 1. 6. Вплив вхідних параметрів на стаціонарні режими. Побудова статичних характеристик різ
- •2. 2. Реактор ідеального витиснення (рів)
- •2. 2. 1. Математична модель рів
- •2. 2. 2. Дослідження рів.
- •1.Зміна ступені перетворення при iзотермічному режимі
- •2.Зміна ступеня перетворення при адiабатичному режимі
- •2. 2. 3. Зіставлення різ та рів
- •2. 3. 5. Ячеєчна модель
- •2. 4. Дифузійна модель зподовжнім переносом речовини та тепла
- •2.5. Двохпараметрична дифузійна модель
- •3. Побудова математичнОї моделі каталітичного реактора
- •3. 1. Етапи побудови математичної моделі
- •3. 2. Структурний аналіз
- •3. 3. Моделювання процесу на одному зерні каталiзатора
- •3.4 Теоретична оптимiзація.
- •3. 5. Попередній вибір типу реактора .
- •3. 6. Моделювання процесу в шару каталiзатора.
- •4. Усталеність реакторних схем
- •4.1 Методи дослідження усталеності
- •4.2 Усталеність простих схем
- •4.3 Усталеність промислових реакторів.
- •5. Методи оптимізації технологічних процесів
- •5.1. Постановказадачіоптимізації
- •5.2. Цільова функція і її властивості
- •5.2.1. Нормалізація незалежних перемінних
- •5.2.2. Геометрична інтерпретація цільової функції
- •5.2.3. Особливі крапки і лінії цільової функції
- •5.2.4. Глобальний і локальний оптимуми
- •5.3. Методи рішення задач оптимізації
- •5.4.Аналітичні засоби
- •5.5. Загальна характеристика засобів рішення задач нелiнійного програмування
- •5.6. Градiєнтні методи рішення задач оптимiзації
- •5.6.1. Градієнт цільової функції
- •5.6.2. Обчислення похідних цільової функції
- •5.6.3. Засіб релаксації.
- •5.6.4. Метод градієнту
- •5. 6. 5. Засіб найскорішого спуска
- •5.7. Безградiєнтні методи рішення задач оптимiзації
- •5. 7. 1. Метод сканiрованiя
- •5. 7. 2. Метод локалiзації екстремума
- •5. 7. 3. Метод "золотого перетину"
- •5. 7. 4. Метод покоординатного спуска Гаусса - Зейделя
- •5. 7. 5. Метод Хука - Джiвса
- •5. 7. 6. Метод сканiрованiя
- •5. 7. 7. Симплексний метод
- •5.7.8. Метод Нелдера-Мида
- •5.8. Методи випадкового пошуку
- •5.8.1. Метод сліпого пошуку
- •5.8.2. Метод випадкових напрямків
- •5.8.3. Метод випадкових напрямків зі зворотним кроком
- •5.8.4. Одержання випадкових чисел
- •5.8.4.1. Метод добутків
- •5.8.4.2. Метод відрахувань
- •5.8.4.3. Одержання псевдовипадкових послідовностей з ірраціональних чисел
- •5.9. Порівняння різних методів рішення задач оптимізації методами нелінійного програмування
- •Література
3. 5. Попередній вибір типу реактора .
У якості прикладу розглянемо вибір типу апарату для процесу одержання формальдегiда із метанола. В наданому процесі протікає дві послідовні необоротні екзотермічні реакції, причому енергія активації другої побочної` реакції менш енергії активації основної реакції. В цьому випадку теоретично оптимальним є ізотермичний режим при максимально можливій температурі, що визначається умовами термостійкістi каталiзатора.
Моделювання було проведено для трьох різноманітних апаратів : трубчатий реактор з безупинним відведенням тепла (а), апарат з трьома адiабатично працюючими верствами каталiзатора та відведенням тепла між ними (б), апарат , що комбінувався послідовно, що становить послідовно розміщені трубчату частину і одну адіабатичну працюючу верству каталiзатора (в)(Риc. 3. 11).
Профілі температур для кожного із разгляданих апаратів приведені на Риc. 3. 12.
Як видно з приведених графіків, найкращим образом наблизизитися до теоретично оптимального режиму можна в апараті, що комбінувався. Параметрами, що визначаються вибір реактора, є жарактер і число фаз, потужність реактора, активність каталiзатора, величина адiабатичного розігріва, температурний інтервал допущеної праці каталiзатора, наявність дифузійного гальмування , зміна активності каталiзатора в часу та ін. .
3. 6. Моделювання процесу в шару каталiзатора.
При моделюваннi процесів в шару необхідно враховувати ряд особливостей :
1. Наявність двох фаз - твердої - зерна каталiзатора та газоподібної (жiдкой)- струм реакційної суміші.
2. Температура і концентрації компонентів в газовій і твердій фазах можуть розрізнятися надто значно.
3. Перенос тепла по шару здійснюється як по твердій фазі, так і потоку реагентов, а перенос речовини - тільки по потоку реагентів.
Для більшості промислових процесів в залежності від співвідношення висоти шара каталiзатора(l)до діаметру(d)стає можливим не враховувати оті чи інші елементарні процеси(рис. 3.13). Так, якщо L>d, оте істотне значення мають процеси переносу в радiальному напрямку і передачі тепла крізь стінку апарату, а при L>>d, оте, як правило, можна зневажити процесами переносу в продольному напрямку за рахунок дифузії і теплопровідностi. Процес в апараті такого типу може бути описаний двухпараметричною діфузіоною моделлю за умови, що D=0 .
Для апаратів, коли L<d, картина міняється: процесами радіального переносу і теплоотвода крізь стінку часто можна заневажити, а процеси подовжнього переносу грають надто помітну роль. Розглянемо понад докладно останній випадок.Якщо процеси переносу між зовнішньою поверхню зерна каталiзатора і ядром струму не виявляють помітного впливу на показники процесу, відбувається інтенсивний тепло і масообмін поміж фазами , оте для опису процесу в цьому випадку можна використовувати дифузійну модель з одним параметром. В разі достатньо більших лінійних швидкостей потоку та високих шар каталізатора укладення дифузійної , що складає , надто незначтний і для опису процесу в шару каталiзатора можна використовувати модель РІВ . Коли вплив зовнішнього тепло- і масообміна істотно, оте в цьому випадку використовується двохфазна модель для опису процесів в шару каталізатора. Схематично надану модель можна представити слідуючим чином(рис.3.13).
Шар ділиться на дві частини - тверда фаза (зерна каталiзатора)і вільний обсяг , по якому проходить весь струм реагентів . Хімічне перетворення діється тільки в твердій фазі. При висновку рівнянь математичної моделі прийняті наступні позначення:
- частка вільного обсягу шару та зерна; долі;
dF- поверхня обміну між фазами у виділеному елементарному обсязі шаруючи висотою dl, м2;
K - ефективний коефіцієнт теплопровідності по твердій фазі, вт/(м.К);
m,T - коефіцієнти массо- і теплообміну, м/с і вт/(м2.К), відповідно;
Сз ,Тз - концентрація і температура на зерні каталізатора;
Сk - теплоємність каталізатора, Дж/(м3.К);
Wн, W - швидкість хімічного перетворення, що спостерігається, віднесена до одиниці об'єму зерна каталізатору чи шару, моль/(м3.с).
Інші позначенняаналогичні дифузійної моделі.
Математичний опис процесу має наступний вид :
1. Матеріальний баланс в газовій фазі враховує процеси переносу маси за рахунок дифузійного і конвектівного струмів і масообмін між фазами. Розглянемо докладно висновок цього рівняння :
- Прихід речовини з струмом
-Витрата речовини з струмом
- Прихід речовини за рахунок dC потоком продольної діфузії
- Видаток речовини за рахунок продольной дифузії
- Передається каталiзатору за рахунок масообміну
m*dF*C
- Приймається від каталiзатора за рахунок масообміну
m*dF*Cз
Зміна кількості речовини в вільному обсязі шара * S * dl
Поділивши всі доданки на S dl і привівши подібні,одержимо в остаточному
виді дорівнення матеріального балансу для газової фази :
(3.19)
2. Тепловий баланс в газовій фазі враховує процеси переносу тепла за рахунок продольної теплопроводностi, конвектiвного потоку і теплообміну між фазами. Висновок аналогичен попередньому.
(3.20)
3. Матеріальнийбаланс в твердій фазі враховує перенос речовини між фазами зарахунокмасообміну та хімічне перетворення :
Позначивши W= та = одержимо:
(3.21)
4. Тепловийбаланс в твердій фазі враховує процеси переносу тепла за рахунок теплопроводностi і теплообміна між фазами , виділення тепла за рахунок хімічної реакції
(3.22)
Система рівнянь становить собою двохфазну математичну модель процесу в непорушному шарі каталiзатора. Надану модель використовують при моделюваннi процесів , для яких істотне значення має зовнішній тепло- і масообмін, а також для для розрахунку перехідних режимів.