- •Моделирование химико-технологических процессов
- •Оглавление
- •Введение, основные понятия и определения
- •1. Методы моделирования и область их применения
- •1.1. Физическое моделирование (фм)
- •1.2. Математическое моделирование (мм)
- •2. Основные принципы и направления при построении и решении математических моделей
- •2.1. Составление математической модели
- •2.2. Нахождение решения математической модели
- •2.3. Проверка моделей на адекватность
- •3. Математическое описание структуры потока в аппарате – основа построения моделей
- •3.1. Методы исследования структуры потоков
- •3.2. Основные характеристики функции распределения потока по времени пребывания в аппарате
- •4. Типовые модели структуры потока
- •4.1. Модель идеального перемешивания
- •4.2. Модель идеального вытеснения
- •4.3. Ячеечная модель
- •4.4. Ячеечная модель с рециркуляцией
- •4.5. Диффузионная модель
- •4.6. Сравнение аппаратов соответствующих моделям ип и ив
- •5. Методы статистического анализа эксперимента
- •5.1. Основные характеристики случайных величин
- •5.2. Равномерное распределение
- •5.3. Нормальное распределение
- •5.4. Доверительные интервалы и доверительная вероятность, распределение Стьюдента
- •5.5. Определение общей дисперсии для серии параллельных опытов
- •5.6. Оценка дисперсии нормально распределенной случайной величины
- •5.7. Проверка однородности результатов измерений
- •6. Планирование эксперимента
- •6.1. Основные понятия и определения
- •6.2. Выбор области проведения эксперимента
- •6.3. Полный факторный эксперимент (пфэ)
- •6.4. Дробный факторный эксперимент
- •7. Оптимизация эксперимента
- •8. Планы второго порядка
- •8.1. Центральное композиционное планирование
- •8.2. Ортогональный план второго порядка
- •8.3. Ротатабельный план второго порядка
- •Список литературы
- •Приложение 1. Квантили нормального распределения
- •Приложение 2. Квантили распределения Пирсона
- •Приложение 3. Значения параметра для различных уровней значимости и степеней свободы
- •Приложение 4. Квантили распределения Стьюдента
- •Приложение 6. Пример использования модели ип для описания процесса непрерывной массовой кристаллизации
- •Приложение 7. Примеры использования типовых моделей для описания процесса теплообмена
- •Лицензия лр № 020370
4.2. Модель идеального вытеснения
Модель ИВ представляет теоретическую модель с идеализированной структурой движущегося потока. В соответствии с ней принимается поршневое течение потока без продольного перемешивания при равномерном распределении концентрации вещества в направлении, перпендикулярном его движению. Время пребывания всех частиц в системе одинаково и определяется как или , где L – длина аппарата, W – скорость потока в аппарате.
Модели ИВ наиболее всего соответствуют процессы, протекающие в аппаратах трубчатого типа: трубчатые реакторы, теплообменники труба в трубе и другие аналогичные аппараты с отношением длины к диаметру L/d > 20 при Re > 2320. Принципиальная схема модели представлена на рис. 4.3.
СвхWF
СвыхWF
Рис. 4.3. Принципиальная схема модели ИВ
Получим математическое описание модели аппарата ИВ из уравнения материального баланса относительно элементарного объема, равного произведению площади поперечного сечения на толщину элементарного слоя потока Fdz.
При наличии каких-либо превращений или в условиях нестационарного режима изменение количества вещества в элементарном объеме будет равно разности между количеством вещества, поступающего в данный объем и выходящего из него
, (4.4) |
где t – текущее время; z – координата, вдоль которой перемещается поток; F – площадь поперечного сечения аппарата.
Раскрыв скобки после некоторых преобразований, получим
(4.5) |
Для стационарных условий и при отсутствии каких-либо превращений в аппарате Свх = Свых.
Решение модели, удовлетворяющее начальному условию,
C(0, z) = Cн(z) при t = 0 , 0 < z < L; (4.6) |
и граничному условию
C(t, 0) = Cвх(t) при z = 0, t > 0, (4.7) |
есть
. (4.8) |
Рис. 4.4. Отклики модели на типовые возмущения
Из решения системы следует, что любое возмущение на входе повторяется на выходе через время, равное среднему времени пребывания . Отклики модели на типовые возмущения приведены на рис. 4.4.
4.3. Ячеечная модель
Ячеечная модель впервые была предложена для каскада реакторов с мешалками и является одной из самых простых среди реальных моделей. В этом случае аппарат представляют состоящим из ряда n последовательно соединенных ячеек, через которые проходит поток вещества. При этом принимается, что в каждой из ячеек поток идеально перемешан, а между ячейками перемешивание отсутствует. Параметром ячеечной модели, количественно характеризующим продольное перемешивание, служит число ячеек полного перемешивания. При n , ЯМ ИВ, а при n 1, ЯМ ИП.
Ячеечная модель достаточно точно воспроизводит структуру потока в последовательно соединенных аппаратах с мешалками (каскад реакторов), в массообменных колоннах с безпровальными тарелками и, частично, в кипящем слое. При внесении соответствующих изменений в ЯМ (ячеечная с рециркуляцией) она может использоваться и для аппаратов с обратным перемешиванием потока: массообменные колонны с провальными тарелками, барботажные колонны, аппараты с кипящим слоем и т.д.
Принципиальная схема модели представлена на рис. 4.5.
Свх,
Рис. 4.5. Принципиальная схема ячеечной модели
Математическое описание модели для случая, когда объемы ячеек и среднее время пребывания потока в каждой из них равны, имеет вид
|
- - - - - - - - - - - - - - - - - |
(4.9) |
- - - - - - - - - - - - - - - - - |
|
где – объем i-й ячейки.
В условиях стационарного режима и отсутствия каких-либо превращений в аппарате Свх = Свых = Сn. Отклики модели на типовые возмущения приведены на рис. 4.6.
Рис. 4.6. Отклики модели на типовые возмущения
Решения модели:
Импульсное возмущение
Для 1-й ячейки в соответствии с граничными условиями: Свх = 0 при t = 0 и С1 = Сн на основании решения модели ИП
. (4.10) |
Тогда для 2-й ячейки в соответствии с граничными условиями: Свх = С1 и С2 = 0 при t = 0
. (4.11) |
Произведя аналогичные вычисления для всех ячеек, для n-ной будем иметь
. (4.12) |
Введя безразмерную концентрацию С() = Сn/Сн и время , функцию отклика (4.12) можно представить в виде
. (4.13) |
Ступенчатое возмущение
При ступенчатом возмущении для случая скачкообразного уменьшения концентрации до нуля аналогично получаем
. (4.14) |
При ступенчатом возмущении для случая скачкообразного увеличения концентрации
. (4.15) |
Оценка параметра n ячеечной модели
Параметр n можно определить через моменты функции отклика на импульсное возмущение:
Начальный момент 2-го порядка . |
|
Центральный момент 2-го порядка . (4.16) |
|
Безразмерный центральный момент 2-го порядка |
. |
Рассмотрим пример определения параметров и n ячеечной модели.
Для моделирования процесса в аппарате было решено использовать ячеечную модель. Структура потока в аппарате исследована импульсным методом, результаты исследования и расчетов приведены в табл. 4.1. Требуется определить целесообразность использования ячеечной модели.
Таблица 4.1
Результаты исследования структуры потока и расчета параметров модели
t, мин |
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
Cэ(t), г/л |
0 |
0,25 |
0,7 |
1,05 |
1,05 |
0,55 |
0,25 |
0,1 |
0,05 |
0,01 |
0 |
С(t), мин–1 |
0 |
0,062 |
0,175 |
0,262 |
0,262 |
0,137 |
0,062 |
0,025 |
0,012 |
0,0025 |
0 |
|
0 |
0,278 |
0,555 |
0,833 |
1,111 |
1,389 |
1,666 |
1,944 |
2,222 |
2,499 |
2,777 |
C() |
0 |
0,225 |
0,629 |
0,943 |
0,943 |
0,494 |
0,225 |
0,090 |
0,045 |
0,009 |
0 |
при n = 6 |
0 |
0,121 |
0,734 |
1,053 |
0,838 |
0,483 |
0,227 |
0,093 |
0,034 |
0,012 |
0,004 |
Решение.
Среднее время пребывания индикатора в потоке, мин:
.
Значения нормированной и безразмерной кривых
; .
Безразмерная дисперсия и количество ячеек
.
Принимаем n = 6. Тогда расчетные значения безразмерной кривой при шести ячейках
.
Среднее значение безразмерной кривой
.
Дисперсия относительно среднего значения безразмерной кривой
.
Дисперсия адекватности
.
Расчетное значение критерия Фишера
,
где критическое значение Fк взято по таблице (приложение 5) для 5 %-ного уровня значимости. Вывод: Fк < F – модель использовать целесообразно.