
- •Комп’ютерний практикум 2014
- •1Аналіз статики і динаміки адсорбції
- •1.1Загальні положення
- •1.2Математичні моделі адсорбційних процесів
- •1.3Розрахункова частина
- •1.3.1Обробка експериментальних даних з рівноваги адсорбції в середовищі Excel
- •1.3.2Розрахунок ізотерми адсорбції, хвильової швидкості та розв’язок хвильового рівняння в середовищі MathCad.
- •1.4Додатковий теоретичний матеріал з динаміки адсорбції [3]
- •1.5Контрольні питання
- •1.6 Література
- •Загальні положення кінетики Ленгмюра-Хіншельвуда
- •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.3Розрахункова частина
1.3.1Обробка експериментальних даних з рівноваги адсорбції в середовищі Excel
Вихідні дані для статичних умов
Очищення газу від сірководню відбувається на адсорбенті оксид цинку при тиску 3,5 МПа і температурі 380оС. Змодельовані в статичних умовах наступні залежності адсорбованого сірководню q (г/100 г адсорбенту) від його концентрації в газі с (мг/м3):
c, мг/м3 |
0 |
8 |
16 |
24 |
32 |
40 |
48 |
56 |
64 |
72 |
80 |
q, г/100г |
0 |
6,48 |
10,29 |
12,8 |
14,58 |
15,91 |
16,94 |
17,75 |
18,42 |
18,98 |
19,44 |
Максимальна концентрація сірководню у газі, мг/м3, смакс=80.
Адсорбційна ємність оксиду цинку, г/100 г адсорбенту, qмакс=19,44.
Вихідні дані для динамічних умов процесу
Лінійна швидкість газу при робочих умовах: w=360 м/год.
Висота шару адсорбенту L=4,8 м.
Час роботи адсорберу у міжрегенераційний період τр=7200 год.
Концентрація сірководню у вихідному газі, мг/м3, свих=50.
Порозність адсорбенту p=0,5.
Насипна щільність адсорбенту, т/м3 ρн=1,3.
Алгоритм розрахунку
Ввести вихідні дані з таблиці в 2 стовпчика.
Побудувати діаграму 1 (точкова тут і далі) за наведеними даними: q=f(c).
Перерахувати в наступні 2 стовпчики вихідні дані у безрозмірні відносні значення: С=с/с0=с/смакс; Q=q/q0=q/qмакс. Всі обчислення виконуються в стовпчиках з використанням авто заповнення («протягування»).
Побудувати діаграму 2 за розрахованими даними: Q=f(С) (ряд 1).
Рівняння Ленгмюра (13) може бути представлено у вигляді (довести це):
. Підготувати відповідну таблицю і побудувати діаграму 3:
-1=f(
). Додати лінію тренда, показати лінійне рівняння на діаграмі (y=ax+b).
Відмітимо, що коефіцієнт b не дорівнює 0. Це пояснюється математичною неможливістю використати точку таблиці з координатами 0,0 для визначення k. Отже k=1/а. Розрахувати за рівнянням Ленгмюра (13) Qрозр=f(С) і нанести криву на діаграму 2 (ряд 2).
Рівняння обмінної адсорбції (13а) може бути представлено у вигляді (довести це):
. Підготувати відповідну таблицю і побудувати діаграму 3:
=f(
). Додати лінію тренда, показати лінійне рівняння на діаграмі (y=ax+b). Звернути увагу на значення коефіцієнта b. Розрахувати k=1/а та Qрозр=f(С) за рівнянням (13а) і нанести криву на діаграму 2 (ряд 3). Оцінити якість апроксимації в порівнянні з іншими рядами.
Розрахувати розподільче співвідношення Rat.
Розподільче співвідношення: є безрозмірним, тому необхідно привести свих і qмакс до формально однакової розмірності кг/м3: со= свих·10-6; q0= qмакс·10∙ρн, після чого розрахувати Rat.
Розрахувати (автозаповнення) значення хвильової швидкості Vwave(C) за рівнянням (18а) (константа рівноваги дорівнює значенню k, розрахованому в пункті 7):
з використанням значень С за пунктом 3. Рекомендується розрахунок провести в кілька стовпчиків з розбивкою рівняння на кілька складових. Побудувати діаграму 3 залежності Vwave(С)=f(С).
Створити стовпчик відносної висоти шару адсорбенту Z від 0 до 1 з кроком 0,1.
Розрахувати С=f(Z) (авто заповнення) – початковий розподіл відносних концентрацій вздовж адсорбера за рівнянням (21):
де К задається від 5 до 10 – значення К погодити з викладачем.
Розрахувати (автозаповнення) значення хвильової швидкості Vwave(Z) за рівнянням (18а) (константа рівноваги дорівнює значенню k, розрахованому в пункті 7):
,
з використанням значень С за пунктом 11. Рекомендується розрахунок провести в кілька стовпчиків з розбивкою рівняння (18а) на кілька складових. Побудувати діаграму 4
Vwave(Z )=f(Z).
Розв’язати хвильове рівняння методом характеристик (створити опис просторово-часового розподілу концентрації в адсорбері).
Незалежна змінна – час роботи адсорберу в безрозмірній формі
t= τр w/L.
Згідно цієї формули і вихідних даних tmin:=0; tmax = 7200∙360/4,8 = 540000
Для моментів часу t0 = 0; t1 =270000; t2 =540000 у три стовпчика розрахувати (автозаповнення) рівняння характеристики (19):
Zt = Z + t∙Vwave(Z,k),
при цьому значення Z використати за пунктом 10.
Побудувати діаграму залежності концентрацій вздовж осі адсорберу С=f(Z) від Zt0, Zt1, Zt2 (три ряди). Цифри 0, 1, 2 відповідають моментам безрозмірного часу t0, t1 t2. Значення С для трьох рядів діаграми використати за пунктом 11.
Індивідуальна самостійна робота
Згідно вище наведеного алгоритму виконати наступне завдання.
Очищення газу від сірководню відбувається на адсорбенті оксид цинку при тиску 3 МПа і температурі 350оС. Змодельовані в статичних умовах наступні залежності адсорбованого сірководню q (г/100 г адсорбенту) від його концентрації в газі с (мг/м3):
c, мг/м3 |
0 |
8 |
16 |
24 |
32 |
40 |
48 |
56 |
64 |
72 |
80 |
q, г/100 г |
0 |
4,17 |
7,14 |
9,38 |
11,11 |
12,5 |
13,64 |
14,58 |
15,38 |
16,07 |
16,67 |
Максимальна концентрація сірководню у газі, мг/м3, смакс=80.
Адсорбційна ємність оксиду цинку, г/100 г адсорбенту, qмакс=16,67.
Концентрація сірководню у вихідному газі, мг/м3, свих=40.
Всі інші згідно основного завдання.
Порівняти отримані результати з основним розрахунком і зробити висновки.
2) Згідно вище наведеного алгоритму, починаючи з п.4, виконати наступне завдання.
Очищення газу від сірководню відбувається на адсорбенті оксид цинку при тиску 3 МПа і температурі 250оС. Змодельовані в статичних умовах наступні залежності адсорбованого сірководню Q (відносна концентрація, частки) від його відносної концентрації в газі C (частки):
-
C
0,2
0,3
0,4
0,5
0,6
0,7
0,8
0,9
1
Q
0,111
0,176
0,25
0,333
0,429
0,538
0,667
0,818
1
Концентрація сірководню у газі, мг/м3, свих=30.
Адсорбційна ємність оксиду цинку, г/100 г адсорбенту, qмакс=15.
Всі інші згідно основного завдання.
Порівняти отримані результати з основним розрахунком і зробити висновки.