книги / Моделирование химико-технологических процессов
..pdfРешения модели:
Импульсное возмущение
Для граничных условий |
Свх=0 |
и СВых = Сн = G/V при t = О |
|
dC |
С( dC |
‘,dt |
, С |
— = (СВХ-С)/?-». |
J — = -h r - » ln — —т —^ С = Сн ехр| |
||
at |
г С |
A t |
Cu |
где Сни G - начальная концентрация в потоке и общее количество индика тора соответственно.
Ступенчатое возмущение
Для граничных условий СВЬ1Х = 0 |
и Свх= const при t = О |
|||
</С |
= j * _ > in c ^ z £ ! |
: - г - > |
||
^ = (Свх- С ) / г |
|
|||
О^вх |
|
J |
7 |
t |
|
^ |
|
||
|
|
|
|
(4.3) |
C = CR l - |
ехр^ -у |
|
Основным параметром модели ИП является среднее время пребыва
ния t
Отклики модели на импульсное и ступенчатое возмущения приведе ны на рис. 4.2.
Рис. 4.2. Отклики модели ИП на типовые возмущения
4.2. Модель идеального вытеснения
Модель ИВ представляет теоретическую модель с идеализированной структурой движущегося потока. В соответствии с ней принимается поршневое течение потока без продольного перемешивания при равно мерном распределении концентрации вещества в направлении, перпенди кулярном его движению. Время пребывания всех частиц в системе одина ково и определяется как t = LjW или ( = F/ F, где L - длина аппарата, W -
скорость потока в аппарате.
Модели ИВ наиболее всего соответствуют процессы, протекающие в аппаратах трубчатого типа: трубчатые реакторы, теплообменники труба в трубе и другие аналогичные аппараты с отношением длины к диаметру L/d> 20 при Re > 2320. Принципиальная схема модели представлена на рис. 4.3.
- ^ 1 1 -
CBXWF |
CWF |
( с + — dzWF |
CBhlxWF |
У & ) |
|||
--------- ► |
|
W |
|
-------------------------- L --------------------------
Рис. 4.3. Принципиальная схема модели ИВ
Из решения системы следует, что любое возмущение на входе по вторяется на выходе через время, равное среднему времени пребывания I = L/W . Отклики модели на типовые возмущения приведены на рис. 4.4.
4.3. Ячеечная модель
Ячеечная модель впервые была предложена для каскада реакторов с мешалками и является одной из самых простых среди реальных моделей. В этом случае аппарат представляют состоящим из ряда п последователь но соединенных ячеек, через которые проходит поток вещества. При этом принимается, что в каждой из ячеек поток идеально перемешан, а между ячейками перемешивание отсутствует. Параметром ячеечной модели, ко личественно характеризующим продольное перемешивание, служит число ячеек полного перемешивания. При п —> оо, ЯМ -> ИВ, а при п -> 1, ЯМ -> ИП.
Ячеечная модель достаточно точно воспроизводит структуру потока в последовательно соединенных аппаратах с мешалками (каскад реакто ров), в массообменных колоннах с безпровальными тарелками и, частич но, в кипящем слое. При внесении соответствующих изменений в ЯМ (ячеечная с рециркуляцией) она может использоваться и для аппаратов с обратным перемешиванием потока: массообменные колонны с проваль ными тарелками, барботажные колонны, аппараты с кипящим слоем и т.д.
Принципиальная схема модели представлена на рис. 4.5.
Рис. 4.5. Принципиальная схема ячеечной модели
Математическое описание модели для случая, когда объемы ячеек и среднее время пребывания потока в каждой из них равны, имеет вид
|
_ ^вх Q |
|
|
dt |
t |
|
|
dCi _ C j-i-C j |
(4.9) |
||
dt |
t |
||
|
|||
dCn _ Q?-i ~ Q? |
|
||
dt |
t |
|
где t =VjV\ Vf- объем /-й ячейки.
В условиях стационарного режима и отсутствия каких-либо превра щений в аппарате Свх = СВЬ1Х = Сп. Отклики модели на типовые возмуще ния приведены на рис. 4.6.
Рис. 4.6. Отклики модели на типовые возмущения
Решения модели:
Импульсное возмущение
Для 1-й ячейки в соответствии с граничными условиями: Свх = 0 при
t = 0 и С\ ~ Сн на основании решения модели ИП |
|
С ,= С нехр(-г/Г). |
(4.10) |
Тогда для 2-й ячейки в соответствии с граничными |
условиями: |
Свх = Ci и С2 = 0 при t = 0 |
|
C2 =(t/i)C Hexp(-t/t). |
(4.11) |
Произведя аналогичные вычисления для всех ячеек, для «-ной будем
иметь
/7-1
(4.12)
с " ' 7 И ехрГ 7
Введя безразмерную концентрацию С(0) = Сп/Си и время 0 = г/г, функцию отклика (4.12) можно представить в виде
С(0) = япе'1" 1ехр (-и 0). |
(4.13) |
(и -1 )! |
|
Ступенчатое возмущение
При ступенчатом возмущении для случая скачкообразного уменьше ния концентрации до нуля аналогично получаем
(4.14)
С„
При ступенчатом возмущении для случая скачкообразного увеличе ния концентрации.
|
1 |
ехр |
t |
|
+ ...+ |
(4.15) |
|
Си, |
( » - # |
|
i |
Оценка параметра п ячеечной модели
Параметр п можно определить через моменты функции отклика на
импульсное возмущение: |
|
|
|
|
Начальный момент 2-го порядка |
м ; = |
1 + - * |
|
|
Центральный момент 2-го порядка |
l4 = M 2‘ - i 2 = — . |
(4.16) |
||
Безразмерный центральный |
0 |
2 |
1 |
|
момент 2-го порядка |
ia2 =CTe = - |
|
||
|
|
|
|
|
Рассмотрим пример определения параметров t |
и п ячеечной модели. |
Для моделирования процесса в аппарате было решено использовать ячеечную модель. Структура потока в аппарате исследована импульсным методом, результаты исследования и расчетов приведены в табл. 4.1. Тре буется определить целесообразность использования ячеечной модели.
Таблица 4.1
Результаты исследования структуры потока и расчета параметров модели
/, мин |
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
C3(t), г/л |
0 |
0,25 |
0,7 |
1,05 |
1,05 |
0,55 |
0,25 |
0,1 |
0,05 |
0,01 |
0 |
минт -1 0
С(0)
Cp {Q)
при п = 6
Решение.
Среднее время пребывания индикатора в потоке, мин:
ю
10 |
= 3,6. |
1 С Э/ |
|
/=1 |
|
Значения нормированной и безразмерной кривых
с ( 0 = # ^ ; с(е)=гс(/)=з,б -с(/).
/=1
Безразмерная дисперсия и количество ячеек
о е ^ & в, - - ! ) 2< ^ 9 = 0,172 |
и = - ^ = 5,8 . |
/=1 |
CT0 |
Принимаем п = 6. Тогда расчетные значения безразмерной кривой
при шести ячейках
^Л пЛ -1 -п в |
^ 6 Q 5 - 6 0 |
с /’(е)=- 9 е |
6 0 6 |
(и -1 )! 5!
|
|
Z C ? |
|
|
|
|
Се = 1 |
г |
= 0’36- |
|
|
Дисперсия относительно среднего значения безразмерной кривой |
|||||
|
ю t |
|
|
|
|
|
I |
( c |
? - c j |
|
|
_2 _/=1________ |
= 0,137. |
|
|||
ср |
|
10-1 |
|
||
|
|
|
|
|
|
Дисперсия адекватности |
|
|
|
|
|
|
l(c,9- c ,p)2 |
|
|||
_2 |
_/=]_______ |
= 0,005. |
|
||
|
|
10-1 |
|
||
|
|
|
|
||
Расчетное значение критерия Фишера |
|
||||
F = |
= 27,4 > F K= 3,25, |
|
|||
|
°ад |
|
|
|
|
где критическое значение |
F K взято |
по таблице |
(приложение 5) для |
||
5 %-ного уровня значимости. |
Вывод: |
F K < F - |
модель использовать |
||
целесообразно. |
|
|
|
|
|
4.4.Ячеечная модель с рециркуляцией
Ячеечная модель с рециркуляцией (ЯМ Р) является модификацией ячеечной модели. Необходимость создания этой модели заключалась в том, что ЯМ не всегда адекватно воспроизводила структуру потока в ре альном аппарате (например, в аппарате кипящего слоя, колонных аппара тах с провальными тарелками, насадочных аппаратах и т.д.).
Согласно ЯМ Р аппарат рассматривается как последовательность зон с сосредоточенными параметрами, причем каждая из зон эквивалентна ячейке идеального перемешивания. Далее предполагается, что между ячейками существуют обратные токи. Параметрами ЯМ Р являются сред нее время пребывания, количество ячеек и доля обратного тока. Принци пиальная схема модели приведена на рис. 4.7.
Математическое описание модели для случая, когда объемы ячеек, среднее время пребывания потока в каждой из них и количество обратного тока v соответственно равны, имеет вид
rdQ
^ |
= < x + v C 2 - ( F + v ) c l |
|
|
dt |
|
K ^ |
= (F + V) C ,_, + VC,+1 - { v + 2 v )q |
(4.17) |
v ^ = {v+v) c n_} - ( v +v) c n,
где Vi - объем одной ячейки, V- прямой поток вещества в ячейке, v - об ратный поток вещества в ячейке, С/ - концентрация вещества в соответст вующей ячейке.
Если величину v / V, равную доле обратного тока, обозначить как f
то разделив обе части уравнений системы (4.17) на F с учетом t = V/V,
получим
= +yc2-0 + /)c l
(4.18)
i ? f =(1 + /)См+УС,'+1_(1 + 2 / ) С '
Для |
системы |
выполняются |
следующие начальные |
условия: |
С\ = CiH, •••>Cj —CjH, |
С„ — С„„ при |
f —0. |
|
|
При п |
- > оо и / = const ЯМР —> ИВ; при л = const и / —> оо |
ЯМР -> |
||
ИП; при п |
оо и / -» оо ЯМР —> ДМ. |
|
|
В условиях стационарного режима и отсутствия каких-либо превра щений в аппарате Свх = С, = Свых = Сп. Отклики модели на типовые воз мущения представлены на рис. 4.8.
Рис. 4.8. Отклики модели на типовые возмущения
Оценка параметров п и/ ЯМР.
Параметры « и / можно определить, решая систему уравнений:
|
м\ |
4 ~ x 2) - 2 x ( l - x " ) |
(4.19) |
||
|
2(1 - х )2 |
|
|||
|
|
|
|
||
д/Э _ 1 2 |
6x(l + Зхп)+ 3H(I - |
х2) |
12х(\ + x)(l - |
хп) |
|
|
|
п\\- x f |
|
\ - x f |
(4.20) |
3 |
|
' |
’ |
где х =./7(1+Д а М® и - второй и третий начальные моменты безраз мерной С-кривой.