- •Комп’ютерний практикум 2014
- •1Аналіз статики і динаміки адсорбції
- •1.1Загальні положення
- •1.2Математичні моделі адсорбційних процесів
- •1.3Розрахункова частина
- •1.3.1Обробка експериментальних даних з рівноваги адсорбції в середовищі Excel
- •1.3.2Розрахунок ізотерми адсорбції, хвильової швидкості та розв’язок хвильового рівняння в середовищі MathCad.
- •1.4Додатковий теоретичний матеріал з динаміки адсорбції [3]
- •1.6 Література
- •Загальні положення кінетики Ленгмюра-Хіншельвуда
- •1.5Контрольні питання
- •2.1Математична обробка експериментальних даних
- •1.1 Механізм 1: взаємодія адсорбованих компонентів а та в з урахуванням адсорбції продукту реакції с
- •2.2Механізм 2: взаємодія адсорбованого компоненту в із компонентом а із газової фази (ударний механізм Ріділа)
- •2.3Механізм 3: взаємодія адсорбованого компоненту а та адсорбованого за дисоціативним механізмом компоненту в з урахуванням адсорбції продукту реакції с
- •2.4Механізм 4: взаємодія адсорбованого за дисоціативним механізмом компоненту в та компоненту а із газової фази
- •3Порядок виконання роботи
- •4Література
- •5Розрахунок фізичних властивостей суміші газів
- •5.2Розрахунки
- •1. Густина газової суміші
- •2. Теплоємність газової суміші
- •3. Динамічна в’язкість газової суміші
- •4. Кінематична в’язкість газової суміші
- •5. Теплопровідність газової суміші
- •6 Кінематичний коефіцієнт дифузії азоту у газовій суміші
- •5.3Індивідуальна самостійна робота.
- •5.4Література
- •6Макрокінетика гетерогенно-каталітичних процесів. Зовнішньо-дифузійна область: Розрахунок окиснення аміаку
- •6.1Загальні положення
- •6.2Алгоритм розрахунку
- •1. Попередні розрахунки
- •2. Розрахунок кількості каталізатора
- •3.4 Розрахунок числа сіток (використати функцію округлвверх):
- •6.3Індивідуальна самостійна робота
- •6.4Література
- •7Макрокінетика гетерогенно-каталітичних процесів. Внутрішньо-дифузійна область: розрахунок колони синтезу метанолу під середнім тиском
- •7.1Алгоритм розрахунку
- •7.2Індивідуальна самостійна робота
- •7.3Контрольні питання
- •7.4Література
- •8Розрахунок теплообміну в колоні синтезу аміаку
- •8.1Мета роботи
- •8.2Короткі теоретичні відомості
- •8.3Опис середовища розробки
- •8.4Розрахунок
- •Початкові температури газових потоків (за рисунком 1 та 2), °c:
- •8.5Алгоритм розрахунку з коментарями
- •8.6Література.
1.4Додатковий теоретичний матеріал з динаміки адсорбції [3]
Динаміка адсорбції може бути описана за допомогою моделі динаміки рівноважної ізотермічної адсорбції за умовами.
Приймемо, що між концентраціями адсорбтиву та адсорбенту в кожен момент часу та у кожній точці шару існує рівноважне співвідношення. Також виділення теплоти у шарі нескінченно малі, а температури потоку та адсорбенту однакові і процес відбувається за відсутності повздовжньої дифузії.
У відповідності до зроблених припущень процес динаміки сорбції може бути описаний за допомогою системи рівнянь.
|
(24) |
Система (24) має два розв’язки: один за наявності кінцевого градієнту концентрацій у колонці (випадок 1), інший – за, так званого, «обривного» фронту (випадок 2) – фронт такої гіпотетичної форми, за якої концентрація адсорбтиву стрибкоподібно змінюється з с = с0 до
с = 0. Ці розв’язки записуються наступним чином:
Випадок 1:
|
(25) |
Випадок 2:
|
(26) |
де Ui
– швидкість руху концентраційної точки
і шаром адсорбенту,
– похідна ізотерми адсорбції у точці
і,
U – швидкість руху усіх точок «обривного» фронту;
с0 – концентрація адсорбтиву у потоці;
q0 – рівноважна їй величина адсорбції.
Рівняння (25) носить назву закону Вікке, рівняння (26) – закону Вільсона. Дослідимо дію цих законів для різних ізотерм адсорбції.
Опукла ізотерма адсорбції
Ізотерма адсорбції q = f(c) та криві розподілу концентрацій q = q(z) для різних проміжків часу зображені на рисунку 5.
Нехай початковому розподілу відповідає лінія τ0. Оскільки фронт не обривний, то кожна точка цього початкового розподілу буде переміщуватись шаром зі своєю характерною швидкістю, що задається рівнянням (25).
а
– опукла ізотерма; б – криві розподілу
концентрацій у шарі адсорбенту.
Рисунок 5 – Опукла ізотерма адсорбції та розподілення концентрацій у шарі адсорбенту при адсорбції.
Виділимо на ізотермі дві точки (1 та 2), зміщення яких по шару від початкового положення τ0 буде предметом аналізу. Для наведеного на рисунку 5 виду ізотерми (dq/dc)1 > (dq/dc)2, тому, згідно з (25), U1 < U2. Це свідчить про те, що відстань між точками 1 та 2, що початково задана розподілом τ0, зменшується – точки зближуються. Таким чином розподіл τ0 трансформується у розподіл τ1, τ2 і т.п. з поступово зростаючим градієнтом концентрацій. За обривного фронту дія закону Вікке зупиняється та всі точки розподілу (τ5 та т.п.) поширюються з однаковою швидкістю U, що визначається законом Вільсона (26): U1 < U < U2.
Лінійна ізотерма адсорбції
Для лінійної ізотерми адсорбції рівняння (25) співпадає з рівнянням (26). Тому, незалежно від напрямку руху потоку та виду початкового розподілу точки переміщуються шаром з постійною швидкістю U1 = U = U2 та будь-який початковий розподіл зберігає свою конфігурацію.
Увігнута ізотерма адсорбції
Ізотерма та трансформації розподілів наведені на рисунку 6. Приймемо, що початковий розподіл такий самий, як і у попередньому випадку. Згідно з ізотермою, (dq/dc)1 < (dq/dc)2, що, у відповідності з рівнянням (25), дає U1 > U2. Таким чином, у умовах, що розглядаються розподіл концентрацій τ0 має тенденцію до нескінченного «розмивання».
а
– увігнута ізотерма; б – криві розподілу
концентрацій у шарі адсорбенту
Рисунок 6 – Увігнута ізотерма адсорбції та розподілення концентрацій у шарі адсорбенту при адсорбції
Наведений розгляд найпростішої динаміки сорбції дозволяє зробити ряд важливих висновків:
1. Форма ізотерми адсорбції, в залежності від початкових умов та напрямку руху потоку, є фактором, яких спричинює затиснення або розширення фронту сорбції, що переміщується зернистим шаром.
2. Переміщення фронту сорбції зернистим шаром може бути явищем стадійним (наприклад «розмитий» розподіл перетворюється в «обривний»), при чому на кожній стадії проявляється дія своїх характерних законів.

.