О САЙТЕ
Добро пожаловать!

Теперь вы можете поделиться своей работой!

Просто нажмите на значок
O2 Design Template

ФЭА / АИТ / Лабораторная работа №2 «Оптимизация технологического процесса в каскаде идеального смешения»

(автор - student, добавлено - 28-09-2017, 17:32)

Скачать: id-smesheniya.rar [213,33 Kb] (cкачиваний: 8)

 

Лабораторная работа №2

«Оптимизация технологического процесса

в каскаде идеального смешения»

 

 

1. Цель и содержание работы

 

Объектом изучения в данной лабораторной работе является установившийся процесс протекания реакции в группе аппаратов идеального перемешивания.

 

Цель работы заключается в следующем:

 

·закрепление полученных в теоретической части курса "Моделирование технологических процессов" знаний о постановке, формализации и методах решения оптимизационных задач для каскада реакторов идеального смешения;

·получение практических навыков решения реальных оптимизационных задач на ЭВМ;

·получение количественных представлений о влиянии важнейших режимных параметров процесса на совокупность оптимизирующих факторов.

 

2. Сведения из теории

 

Математические модели химических реакторов строятся на основе гидродинамических моделей, учитывающих характер распределения времени пребывания частиц потока реагирующей смеси в данном реакторе с добавлением уравнений химической кинетики.

 

2.1. Основные понятия химической кинетики

 

Скорость химической реакции - изменение числа молей реагентов в результате химического взаимодействия в единицу времени на единицу объема:

 

. (2.1.1)

 

Согласно уравнению (2.1.1), вводя концентрацию с, получим:

 

, (2.1.2)

или . (2.1.3)

 

Знак "+" указывает, что в реакции вещество накапливается, знак "-" - что концентрация вещества снижается.

Константа скорости характеризует процесс, протекающий на микроуровне взаимодействия молекул, и зависит от вида молекул, вступающих в реакцию веществ, и от температуры.

 

.

 

Порядок реакции представляет собой сумму показателей степеней при концентрациях реагирующих веществ. Скорость реакции многих компонентов аппроксимируется уравнением:

 

, (2.1.4)

 

в котором , где - порядок реакции по компоненту , - по компоненту , - общий порядок.

Степень превращения: xA=(cA0 - cA)/cA0, где cA0 и cA – начальная и конечная концентрации реагирующего вещества.

2.2. Типовые математические модели структуры потоков в аппаратах.

 

 

Модель идеального смешения.

 

Условия физической реализуемости этой модели выполняются, если во всем потоке происходит полное смешение частиц потока. Любое изменение концентрации вещества на входе потока в зону идеального смешения мгновенно распределяется по всему объему. Уравнение изменения концентрации имеет вид:

 

, (2.2.1)

где:

- скорость потока, поступающего в зону идеального смешения, м3/ч;

- объем зоны идеального смешения, м3;

Выходные кривые при ступенчатом и импульсном возмущении изображены в табл. 2.2.1. Модели идеального перемешивания соответствует апериодическое звено.

 

 

Модель идеального вытеснения.

 

В соответствии с этой моделью принимается поршневое течение без перемешивания вдоль потока при равномерном распределении вещества в направлении, перпендикулярном движению. Время пребывания в системе всех частиц одинаково.

Математическое описание модели имеет вид:

 

, (2.2.2)

где:

- концентрация вещества;

- время пребывания частиц в реакторе;

- линейная скорость потока;

- координата.

Модели идеального вытеснения соответствует звено чистого запаздывания.

 

 

Диффузионная модель.

 

Различают однопараметрическую и двухпараметрическую модели.

Однопараметрическая модель. Ее основой является модель вытеснения, осложненная обратным перемешиванием. Основным параметром является коэффициент продольного перемешивания .

 

 

(2.2.3)

 

Двухпараметрическая модель. В этой модели учитывается перемешивание потока модели в продольном и радиальном направлениях. Модель характеризуется коэффициентом продольного и радиального перемешивания.

 

 

(2.2.4)

 

 

Каскадная модель.

 

Она представляет собой последовательное соединение аппаратов идеального перемешивания.

Математическое описание включает линейных дифференциальных уравнений:

 

 

, . (2.2.5)

2.3. Закон Аррениуса.

Зависимость константы скорости реакции от температуры представляется в виде закона Аррениуса:

 

. (2.3.1)

 

Модель реактора идеального перемешивания состоит из двух уравнений: гидродинамического и теплового. Однако можно ограничиться одним гидродинамическим уравнением, которое в установившемся режиме имеет вид:

 

,

- время пребывания в реакторе.

 

Для простой химической реакции уравнение может быть записано в следующем виде:

 

, (2.3.2)

 

где:

- температурный коэффициент (константа скорости реакции);

 

 

- абсолютная температура реакции;

- расчетное значение абсолютной температуры;

- постоянный сомножитель температурного коэффициента;

- энергия активации;

- универсальная газовая постоянная.

 

 

Из рис. 2.3.1 видно, что в диапазоне допустимого изменения температуры кривая закона Аррениуса может быть аппроксимирована прямой линией вида:

 

(2.3.3)

 

 

С учетом сказанного выше, обозначив через общее фиксированное время пребывания в каскаде, можем описать процесс в каскаде следующей системой уравнений:

 

(2.3.4)

 

 

 

Рис. 2.3.1. Определение оптимальной температуры.


 

 

 

 

 

Таблица 2.2.1

Типовые модели гидродинамических процессов

 

Типовые модели

Схема потока

Математическое описание

Характер отклика

 

 

 

ступенчатое возмущение

импульсное возмущение

1. Модель идеального вытеснения

2. Модель идеального смешения

 


Продолжение таблицы 2.2.1

 

Типовые модели

Схема потока

Математическое описание

Характер отклика

 

 

 

ступенчатое возмущение

импульсное возмущение

3. Диффузионные модели вытеснения

а) однопараметрическая

б) двухпараметрическая

 

 

 

 

4. Ячеечная модель

 


Таблица 2.2.2

Передаточные функции идеальных типовых моделей

 

Типовые модели

Передаточная функция

Примечание

Модель идеального смешения

, где:

- объем аппарата

- объемная скорость жидкость

Модель идеального вытеснения

- транспортное запаздывание аппарата, где:

- длина аппарата

- линейная скорость жидкости

m аппаратов идеального смешения равного объема в ряд

___

m аппаратов идеального смешения неравного объема в ряд

Модель идеального вытеснения и идеального смешения в ряд


2.4. Источники массы и тепла в потоках.

 

Для учета изменения концентрации вещества в результате химических реакций и массообмена уравнения динамики должны учитывать величину, имеющую смысл интенсивности источников вещества .

Величина может быть положительной или отрицательной, в зависимости от того, образуется или расходуется -ое вещество в потоке.

Частные выражения для источников вещества и тепла:

а) Химическая реакция в потоке: интенсивность источника вещества в уравнениях, характеризующих изменение его концентрации, равна скорости образования этого вещества , т.е.

 

. (2.4.1)

 

б) Массообмен потока с окружающей средой:

 

,

где - скорость массопередачи.

 

, (2.4.2)

 

где:

- поверхность массообмена к единице объема, м23;

- коэффициент массопередачи, м/ч;

- равновесная концентрация вещества в среде, кмоль/м3;

 

в) Теплообмен потока с окружающей средой:

 

,

 

где - скорость теплопередачи.

 

,

 

где:

- поверхность теплообмена, отнесенная к единице объема, м23;

- коэффициент теплопередачи, Дж/(м2×град×ч);

- температура среды, град;

3. Порядок выполнения работы.

 

Решение задач, связанных с отысканием оптимальных условий проведения химических реакций, имеет важное значение при выборе эксплуатационных параметров химических реакторов.

Под оптимальными условиями понимаются оптимальное время пребывания реагентов в реакторе и оптимальная температура в реакторе , обеспечивающие максимальное или минимальное значение заданного критерия Rопт (концентрация продукта на выходе реактора), т.е. для оценки оптимума необходимо выбрать критерии оптимизации (выходной параметр). На основании выбранного критерия оптимизации составляется целевая функция (функция выгоды), представляющая собой зависимость критерия оптимизации от параметров, влияющих на его значение.

Задача оптимизации сводится к нахождению экстремума (в данном случае максимума) целевой функции (x2).

3.1. Постановка оптимальной задачи для каскада

идеального перемешивания

 

Необходимо найти значение и , обеспечивающие экстремум критерия оптимальности при наличии ограничений в виде гидродинамических уравнений при условии:

 

.

Рис. 3.1.1. Каскад идеального перемешивания

 

 

 

3.1.1. Аналитический метод нахождения

оптимального времени пребывания частиц

 

Для реакции типа найдем состав смеси на выходе реактора идеального смешения.

Аналитические выражения для скоростей образования реагентов и имеют вид:

, (3.1.1)

 

. (3.1.2)

 

Стационарный режим определяется системой гидродинамических уравнений (при t1=t2=t):

, (3.1.3)

 

. (3.1.4)

 

В результате решения, которой, находим:

 

, (3.1.5)

 

. (3.1.6)

 

В частном случае из уравнения (3.1.6) получим:

 

. (3.1.7)

 

Критерием оптимальности является максимальная концентрация (степень превращения) продукта Р на выходе реактора. Для критерия оптимальности:

 

, (3.1.8)

=

×

-

K

T

k

e

E

RT

(

)

1

01

/

1

1

где

 

 

=

 


- константы скорости реакции;

t - время пребывания частиц в реакторе;

R – универсальная газовая постоянная;

T1 – температура реагента А;

T2 – температура реагента P.

 

Дифференцируем выражение (3.1.8) по и , находим систему уравнений, определяющих оптимальные условия реакции:

(3.1.9)

 

. (3.1.10)

 

Эта система эквивалентна следующей:

 

(3.1.11)

 

. (3.1.12)

 

Из уравнения (3.1.11) выводим формулу для нахождения оптимального времени пребывания в реакторе при заданном значении температуры:

 

. (3.1.13)

 

Подставляя выражение (3.1.13) в уравнение (3.1.12), имеем:

 

. (3.1.14)

 

Решение уравнения (3.1.12) для заданного времени пребывания имеет вид:

 

. (3.1.15).

Отсутствие абсолютного экстремума у критерия оптимальности (3.1.8) означает, что его наибольшее значение следует искать на границе области допустимых значений переменных и .

3.1.2. Характеристика режима и построение графиков.

 

 

Рис. 3.1.2. Изменение концентрации в зависимости от времени

для реакции при различных значениях концентрации

продуктав исходном сырье.

 

Примечание:

При малом значении функция имеет максимум - рис. 3.1.2. а).

 

На рис. 3.1.2 показано изменение концентраций в зависимости от времени пребывания для различных значений . Величину , соответствующую максимальной концентрации , можно найти двумя методами:

1) аналитическим методом;

2) методом Лагранжа.

Из вида формулы (3.1.13) следует, что с возрастанием температуры оптимальное время уменьшается, и наоборот.

Необходимо выяснить, когда возрастанию выхода продукта способствует повышение, а когда - снижение .

Подставим выражение (3.1.13) в соотношение для критерия оптимальности (3.1.8). В результате получим:

 

. (3.1.16)

 

Продифференцируем (3.1.16) по и приравняем производную к нулю:

 

.

 

Знак производной совпадает со знаком разности .

Если , то , следовательно, выход продукта увеличивается при повышении температуры .

Если , то , следовательно, выход продукта увеличивается при повышении температуры .

Таким образом, если заданы ограничения на температуру , то оптимальная температура будет:

при ;

при .

рассчитывается по формуле (3.1.13).

 

Если заданы ограничения , то оптимальное время пребывания будет:

при ;

при .

Оптимальная температура определяется решением уравнения (3.1.12).

 

3.1.3. Метод Лагранжа.

Пусть требуется найти экстремум функции:

 

, (3.1.17)

 

которая зависит от - переменных , на которые накладываются ограничения:

; (3.1.18)

 

, m<n.

 

Для решения этой задачи вводится вспомогательная функция:

 

, (3.1.19)

где:

- неопределенные множители Лагранжа.

 

Для решения задачи отыскания экстремума функции при ограничениях на переменные (3.1.18) необходимо решить систему уравнений, где функция определяется выражением (3.1.19), совместно с ограничениями (3.1.18):

. (3.1.20)

 

3.2. Краткие сведения из теории.

Выражения (2.3.2) и (2.3.3) играют роль ограничений типа равенств, следовательно, можно вместо условного экстремума функции искать безусловный экстремум обобщенной целевой функции вида:

 

. (3.2.4)

 

В явном виде выражение для концентрации ключевого реагента на выходе из каскада реакторов (xp) получается путем преобразования (3.2.1) и имеет вид:

 

(при t1=t2=t). (3.2.5)

 

Прологарифмировав (3.2.5), получим:

 

. (3.2.6)

 

Введем обозначения:

(3.2.7)

 

С учетом (3.2.7) выражение (3.2.6) примет вид:

 

. (3.2.8)

 

Выражение для ограничения по времени пребывания реагента в каскаде следует из (3.2.7):

 

; (3.2.9)

 

или что, то же самое:

 

, (3.2.10)

 

где:

; ; ; . (3.2.11)

 

С учетом вышеизложенного задача отыскания оптимальных значений переменных в данном случае выглядит следующим образом:

 

.

 

При ограничениях:

Для дальнейших расчетов принять следующие числовые значения:

 

(3.2.12)

; ,

 

где:

; ; ; ; ; ;


3.3. Порядок выполнения работы

Записать обобщенное (через множители Лагранжа) выражение для целевой функции в виде концентрации ключевого реагента на выходе каскада с учетом ограничений типа равенств:

, (3.3.1)

где:

; .

 

· записать обобщенное выражение целевой функции (через множители Лагранжа);

· записать необходимые условия существования экстремума обобщенной целевой функции;

· с помощью ЭВМ решить полученную систему уравнений и найти оптимальные значения переменных, рассчитать оптимальные значения искомых переменных;

· изменить концентрацию выходного потока на каскад xA и решить задачу ранее приведенными способами;

· построить графики зависимости оптимальных значений tопт и Топтот концентрации выходного потока (степени превращения);

· написать отчет о проделанной работе.

 

4. Результаты работы.

Согласно пункту 3.1. по методу Лагранжа:

·решить систему уравнений на ЭВМ;

·определить оптимальные параметры (и ), которые обеспечивают экстремум критерия оптимальности и максимальную концентрацию вещества на выходе;

·построить графики зависимости концентрации от оптимальных значений и ;

·найти и аналитическим методом, сравнить результаты с расчетами, полученными по методу Лагранжа;

·привести сравнительные характеристики оптимальных режимов.

· написать отчет о проделанной работе.

Исходные данные:

; ; ; ;

; ; ; ;

Допустимый предел повышения температуры .

Program laba2;

uses crt;

const T1=312; T2=220; r=8.32; k01=69; k02=55.34; e1=5636*4.2 ;E2=1112*4.2; x0=0.8;

var y1,y2,y3,y4,y5,a0,a1,a2,a3,a4,K1,K2,max:real;

tayMin,tayMax,ctr:integer;

tay,tay1,tay2,ty1,ty2:integer;

begin

ctr:=1;

repeat

clrscr;

if ctr>1 then writeln('Neverniy diapazon!');

write('Vvedite nijniy diapazon vremeni prebivaniya chastici v reaktore: tayMin=');

readln(tayMin);

write('Vvedite verhniy diapazon vremeni prebivaniya chastici v reaktore: tayMax=');

readln(tayMax);

ctr:=ctr+1;

until (tayMin+tayMax=840) and (tayMin>0) and (tayMax>tayMin);

K1:=60*k01*exp(-E1/(R*T1));

K2:=60*k02*exp(-E2/(R*T2));

A0:=ln(x0);

A3:=1/K1;

A4:=1/K2;

A2:=ln(K1);

Max:=0;

for tay1:=tayMin to tayMax do

for tay2:=tayMax downto tayMin do

begin

A1:=ln(tay1+tay2);

y2:=ln(1+tay1*K1);

y3:=ln(1+tay2*K2);

y4:=exp(y2);

y5:=exp(y3);

y1:=a0+a1+a2-y1-y2;

if y1>max then

begin max:=y1;

ty1:=tay1;

ty2:=tay2;

tay:=round(a3*y4+a4*y5-a3-a4);

end; end;

writeln('OptimaLnoe vremya prebivaniya chastic v 1 reaktore:',ty1:6,'sec');

writeln('OptimaLnoe vremya prebivaniya chastic vo 2 reaktore:',ty2:6,'sec');

writeln('OptimaLnoe vremya prebivaniya reagenta v kaskade:',tay:6,'sec');

writeln('Koncentraciya reagenta na vixode:',exp(max):6:3,'moli/l');

readkey;

end.

ОТВЕТ для диапазона 300-540

Оптимальное время пребывания частиц в 1 реакторе: 300 сек

Оптимальное время пребывания частиц во 2 реакторе: 540 сек

Оптимальное время пребывания реагента в каскаде: 840 сек

Концентрация реагента на выходе: 2,223 моль/л.

 

Ключевые слова -



ФНГ ФИМ ФЭА ФЭУ
Copyright 2018. Для правильного отображения сайта рекомендуем обновить Ваш браузер до последней версии!