Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
Предоплата всего
Подписываем
PAGE 2
МИНИСТЕРСТВО ОБРАЗОВНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ
ВОЛГОГРАДСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ
УНИВЕРСИТЕТ
Л А Б О Р А Т О Р Н А Я Р А Б О Т А №-3
по предмету: «Моделирование химических процессов»
на тему: «Расчет на ЭВМ оптимального флегмового числа ректификационной колонны.»
Выполнил:
Студент гр. АВБ – 587
Жуков Б. Г.
Проверил(а):
___________________
Дата: «__»__________
Подпись ___________
Волгоград 2005 год.
Расчет на ЭВМ оптимального флегмового числа
ректификационной колонны.
1 CLS
5 INPUT "gf,xf1,xk1,xd1,a1,b,ma,mb"; gf, xf1, xk1, xd1, a1, b, ma, mb
10 gd = gf * (xf1 - xk1) / (xd1 - xk1)
15 gk = gf - cd
20 xf = (xf1 / ma) / ((xf1 / ma) + (1 - xf1) / mb)
25 xk = (xk1 / ma) / ((xk1 / ma) + (1 - xk1 / mb))
30 xd = (xd1 / ma) / ((xd1 / ma) + (1 - xd1) / mb)
35 f = (xd - xk) / (xf - xk)
40 ypf = EXP(a1 + b * LOG(xf))
45 rm = (xd - ypf) / (ypf - xf)
50 u1 = 10000!
55 rm1 = rm
60 r = beta * rm1
65 ny = 0: nu = 0
70 yp = xd
75 x = EXP((LOG(yp) - a1) / b)
80 dx = x - xf
85 IF dx > 0 THEN GOTO 105
90 x1 = yp * (r + 1 / r - (xp / r))
95 dny = (x1 - xf) / (x1 - x)
100 yn = ny + dny: GOTO 120
105 ny = ny + 1
110 y = r / (r + 1) * x + (xd / (r + 1))
115 yp = y: GOTO 75
120 xu = xk
125 yp = EXP(a1 + b * LOG(xu))
130 xu = ((r + 1 / (r + f)) * yp + ((f - 1) / (r + f)) * xk)
135 dx = xu - xf
140 IF dx <= 0 THEN 160
145 x2 = EXP((LOG(yp) - a1) / b)
150 dnu = (xf - x2) / (xu - x2)
155 un = nu + dnu: GOTO 165
160 nu = nu + 1: GOTO 125
165 z = un + yn
170 u = (r + 1) * z
175 d = u - u1
180 IF d < 0 THEN 200
185 r0 = r / beta
190 PRINT "yn="; un; "un="; un; "r0="; r0; "rm="; rm; "f="; f; "ypf="; ypf
191 PRINT "gk="; gk; "gd="; gd; "z="; z; "u="; u
195 GOTO 210
200 u1 = u
205 rm1 = r: GOTO 60
210 END
Yn=0,79 Un=4,397 R0=2,24
Rm=2,03 F=1,379 YPF=0,61
Gk=11121,5 Gd=6878,505 Z=10,1
U=34,44
Таблица 1
|
Бензол-толуол |
Произв. GF, кг/ч |
Концентрация легколетучего компонента, % |
||
|
XF |
XD |
XW |
||
|
18000 |
35 |
98 |
1,7 |
Число питания
Равновесная концентрация паровой фазы, соответствующая исходной смеси XF=0, 388 YF=0,61
Минимальное флегмовое число
Из графика YF* = 0,61, следовательно, можно определить min флегмовое число:
Число целых теоретических тарелок в укрепляющей части NY=6
Число целых теоретических тарелок в исчерпывающей части Nu=4
Общее число теоретических тарелок в исчерпывающей части UN=4,397
Общее число теоретических тарелок в колонне Z=10,1
Расчетное флегмовое число
R=1,3*Rmin+0,3=1,3*1,68+0,3=2,484
Необходимо перевести массовые концентрации в мольные:
Исходные, расчетные и справочные данные.
|
Наименование параметра |
Размерность |
Обозначение |
Величина параметра |
|
|
Исходные данные |
||||
|
Производительность по исходной смеси |
кг/час |
GF |
GF |
18000 |
|
Абсолютные массовые концентрации легколетучего компонента : В исходной смеси В кубовой жидкости В дистилляте |
Мольные доли |
XF XK Xd |
XF1 XK1 XD1 |
35 1,7 98 |
|
Коэффициенты степенного уравнения, аппроксимирующего равновесие в бинарной смеси |
a |
A B |
0,128 0,73 |
|
|
Молекулярные веса компонентов; Легколетучего Труднолетучего |
кг/моль |
MA MB |
MA MB |
78 92 |
|
Расчетные параметры |
||||
|
Производительность по дистилляту |
кг/час |
Gd |
GD |
6878,505 |
|
Производительность по кубовой жидкости |
кг/час |
Gk |
GK |
11121,5 |
|
Число питания |
F |
F |
1,379 |
|
|
Равновесная концентрация паровой фазы, соответствующая исходной смеси XF |
YPF |
0,61 |
||
|
Минимальное флегмовое число |
RM |
RM |
2,03 |
|
|
Число целых теоретических тарелок в укрепляющей части |
Ny |
NY |
6 |
|
|
Число целых теоретических тарелок в исчерпывающей части |
Nu |
NU |
4 |
|
|
Доля теоретических тарелок в укрепляющей части |
DNY |
|||
|
Общее число теоретических тарелок в укрепляющей части |
yN |
YN |
||
|
Доля теоретических тарелок в исчерпывающей части |
DNU |
|||
|
Общее число теоретических тарелок в исчерпывающей части |
UN |
UN |
4,397 |
|
|
Общее число теоретических тарелок в колонне. |
Z |
Z |
10,1 |
|
|
Функция минимизации |
U |
U |
34,44 |
|
|
Расчетное флегмовое число |
R |
R |
2,484 |
|
|
Оптимальное флегмовое число |
R0 |
R0 |
2,24 |
|
|
Абсолютные мольные концентрации легколетучего компонента: В исходной смеси В кубовой жидкости В дистилляте |
XF XK XD |
XF XK XD |
0,388 0,02 0,983 |
GFXF
Ф
GWXW
GDXD