- •Міністерство освіти і науки україни
- •Скопійовано з оригіналу-макета, наданого автором
- •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. Порівняння різних методів рішення задач оптимізації методами нелінійного програмування
- •Література
2. 1. 1 Модель різ для опису стаціонарного режиму
Стаціонарний режим характеризується значеннями параметрів, які установились, тому для його опису можна користуватися статичною моделлю РІЗ, в якій доданки, що враховують зміну речовини і тепла в часу (праві частини рівнянь) дорівнюються нулю.
Тодірівняння балансів приймуть такий вигляд:
матеріальний баланс :
Q (CBX – C)– VW = 0 (2. 7)
енергетичний баланс:
QCP (TBX – T)+ VW(–H) –S(T – Tx) = 0 ( 2. 8 )
Інші рівняння моделі не змінюються.
Таким чином,статична модель РІЗ становить систему алгебраїчних рівнянь.
2. 1. 2. Модель різ при протечі деяких реакцій
В цьому випадку матеріальний баланс - це система рівнянь, кількість яких дорівнює кількості ключових компонентів.
Якщо ключоваречовина бере участь в деяких реакціях, то швидкість його зміни буде дорівнена сумі швидкостей реакцій, в яких бере участь дана ключова речовина. Рівняння енергетичного балансу тільки одне. Виділення тепла в результаті протечі хімічних реакцій розраховується як сума творів швидкостей реакцій на тепловий ефект відповідної реакції.
_ _
Рівняння кінетики: W1= k1f1(C1) ….. Wn= knfn(C)
2. 1. 3. Дослідження моделі різ
Поводження об'єкту визначаєтьсясукупністю значень параметрів, які входять у математичний опис.
Проведемо класифікацію параметрів моделі РІЗ.
Вхідні параметри :
– технологічні : Tвх, Cвх, Qвх, Tх*;
– конструктивні: V, S*;
– фізичні: a, CP, r, CA;
– параметри, які характеризують реакцію: (–DH), E, B.
Вихідні параметри :
– технологічні: Tвих, C*;
– час :t.
Для знаходження всього трьох вихідних параметрів необхідно задати 13 вхідних параметрів. Для спрощення дослідження моделі і скорочування кількості параметрів застосовують запис рівнянь моделі в безрозмірної формі. Вводиться ряд безрозмірних параметрів (таблиця 2.1).
Таблиця 2. 1
Безрозмірні
|
Розмірні
|
Збільшення
|
Ступінь перетворення Х = (С0– С) / С0 |
Концентрація C= (1–X) С0 |
DC = – C0 dx |
Час t`=t / |
Час t= t` |
dt = dt` |
Температура = (T – T0)/ (b T0) |
Температура T = T0 (1 + b) |
DT = bT0d |
Швидкість реакції |
Швидкість реакції |
|
Константа швидкісті реакції |
Константа швидкісті реакції |
|
Внаслідок одержуємо модель РІЗ в безрозмІрному виді :
- Матеріальний баланс :
(2.9)
- Єнергетичний баланс:
(2.10)
- Рівняння кiнетикi :
(2.11)
де
– відношення усередненої теплоємністi апарату до теплоємністі реакційноі суміші ("безрозмірна теплоємність") ;
- адіабатичний розігрів, тобто величина, що характеризує на скільки градусів розігріється реакційна суміш при повній мірі перетворення;
–адіабатичнийрозігрів в безрозмірной формі;
- безрозмірний параметр тепловідвода;
b= RT0/ E - критерій Аррейнiуса;
X - міра перетворення;
T0- опорна температура, К;
0- константа швидкісті реакціі при опорной температури;
R - газова постійна;
- час перебування (контакту);
С0 – концентрація на вході;
t - віношення поточного значення часу до часу контакту (" безрозмірний час ").
Підрахуємо кількість параметрів в безмірной моделі ;
а) вихідні параметри: - технологічні, X
- час t`.
б) вхідні параметри: - технологічні XBX , BX, X.
в) фізичні-,.
г) параметри,що характеризують реакцію-.
В безмірної моделі 8 параметрів. Таким чином число параметрів скоротилося за рахунок зникнення конструктивних параметрів. Це є важливим достоїнством безмірноі моделі РІЗ.