Оптимизация химико-технологических процессов органического синтеза. Бочкарев В.В - 29 стр.

UptoLike

ВУЗ: 

ТПУ | Томск

Составители: 

Рубрика: 

29
1. Примем значения подинтегральных функций
A
R
и 1/r
A
при
X
A
= 0:
A
R
(0)=1; 1/r
A
(0) = 1.
2. Делаем шаг по степени конверсии X
A
: X
A
(i) = X
A
(i-1)+
X
A
. Значение
X
A
принимаем в пределах 0,001
0,01.
3. Вычисляем последовательно значения C
A
, C
R
, C
S
, C
B
, r
1
, r
2
, r
A
,
R
, 1/r
A
по вышеприведенным формулам.
4. Вычисляем значения интегралов
R
A
A
X
dX
A
0
и
1
0
r
dX
A
A
X
A
, используя
метод трапеций.
5. Рассчитываем выход продукта R X
R
и объемную скорость подачи
исходного сырья W
0
.
6. Проверяем, достигнут ли максимум по X
R
. Если “да”, то запоминаем
оптимальные значения X
R
, X
A
, W
0
.
7. Проверяем условие завершения расчетов X
A
1. Если “нет”, то
возвращаемся к п. 2.
8. Выводим на печать результаты расчетов.
Ниже приведена программа на языке Бейсик, реализующая
приведенный алгоритм вычислений.
CLS
K1 = .0009: K2 = .0011: Ca0 = 1: Cb0 = 1: Vr = .9
dFR = 1: FR = 1: integralF = 0: integralR = 0
interval = .05: maxCR = 0
PRINT " XA CA CR CS "
FOR XA = .001 TO .999 STEP .001
CA = Ca0 * (1 - XA)
CR = Ca0 * XA * FR
CS = (Ca0 - CA - CR) / 2
CB = Cb0 - CR - CS
R1 = K1 * CA * CB: R2 = K2 * CA * CR: Ra = R1 + R2
dFR1 = (R1 - R2) / Ra
integralF = integralF + (dFR + dFR1) / 2 * .001
FR = integralF / XA
dFR = dFR1
CR = Ca0 * XA * FR
CS = (Ca0 - CA - CR) / 2
CB = Cb0 - CR - CS
R1 = K1 * CA * CB: R2 = K2 * CA * CR: Ra = R1 + R2
integralR = integralR + .001 / Ra
W0 = Vr / (Ca0 * integralR)
xR = XA * FR
IF XA < interval GOTO 1
Nitro PDF Trial
www.nitropdf.com

Страницы