
- •Моделирование химико-технологических процессов
- •Оглавление
- •Введение, основные понятия и определения
- •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.4. Ячеечная модель с рециркуляцией
Ячеечная модель с рециркуляцией (ЯМР) является модификацией ячеечной модели. Необходимость создания этой модели заключалась в том, что ЯМ не всегда адекватно воспроизводила структуру потока в реальном аппарате (например, в аппарате кипящего слоя, колонных аппаратах с провальными тарелками, насадочных аппаратах и т.д.).
Согласно ЯМР аппарат рассматривается как последовательность зон с сосредоточенными параметрами, причем каждая из зон эквивалентна ячейке идеального перемешивания. Далее предполагается, что между ячейками существуют обратные токи. Параметрами ЯМР являются среднее время пребывания, количество ячеек и доля обратного тока. Принципиальная схема модели приведена на рис. 4.7.
Рис. 4.7. Принципиальная схема ЯМР
Математическое описание модели для случая, когда объемы ячеек, среднее время пребывания потока в каждой из них и количество обратного тока v соответственно равны, имеет вид
|
- - - - - - - - - - - - - - - - - - - - - - - - - - - - |
|
- - - - - - - - - - - - - - - - - - - - - - - - - - - - |
|
где Vi –
объем одной ячейки,
–
прямой поток вещества в ячейке, v –
обратный поток вещества в ячейке, Сi
– концентрация вещества в соответствующей
ячейке.
Если величину
,
равную доле обратного тока, обозначить
как f, то разделив обе части уравнений
системы (4.17) на
с учетом
,
получим
|
- - - - - - - - - - - - - - - - - - - - - - - - - - - |
|
- - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
Для системы выполняются следующие начальные условия: С1 = С1н, ..., Сj = Сjн , ..., Сn = Cnн при t = 0.
При n и f = const ЯМР ИВ; при n = const и f ЯМР ИП; при n и f ЯМР ДМ.
В условиях стационарного режима и отсутствия каких-либо превращений в аппарате Свх = Сi = Свых = Сn. Отклики модели на типовые возмущения представлены на рис. 4.8.
Рис. 4.8. Отклики модели на типовые возмущения
Оценка параметров n и f ЯМР.
Параметры n и f можно определить, решая систему уравнений:
|
|
|
где x = f/(1+f),
а
–
второй и третий начальные моменты
безразмерной С-кривой.
Связь между числом ячеек n, долей
обратного тока f и безразмерной
дисперсией
можно описать уравнением
|
ЯМР наиболее адекватно описывает структуру потока в насадочных и секционированных колонных аппаратах, где наблюдается заброс вещества в сторону, противоположную направлению основного потока.
4.5. Диффузионная модель
В основе диффузионной модели лежит допущение, что структура потока описывается уравнением, аналогичным уравнению молекулярной диффузии, но в отличие от коэффициента молекулярной диффузии здесь используется коэффициент продольного перемешивания или турбулентной диффузии Dl. Можно также сказать, что основой ДМ является модель ИВ, осложненная обратным перемешиванием. Параметром модели является коэффициент продольного перемешивания Dl. Принципиальная схема модели представлена на рис. 4.9.
Рис. 4.9. Принципиальная схема диффузионной модели
Для получения математического описания
модели составим уравнение материального
баланса для элемента аппарата длиной
dz. В рассматриваемый элемент поступают
конвективный поток FWC
и поток, вызываемый турбулентной
диффузией
,
а покидают его – конвективный поток
и поток, вызываемый турбулентной
диффузией
.
В соответствии с законом сохранения массы разность между входящим и выходящим потоками должна быть равна накоплению вещества в рассматриваемом элементе,
(4.22)
– |
Раскрывая скобки и сокращая подобные члены, получим
|
Данное уравнение является основным уравнением диффузии. Здесь F – площадь сечения аппарата, м2; w – скорость потока, м/с; t – время, с; C – концентрация вещества, кг/м3; Dl – коэффициент продольного перемешивания, м2/с.
При Dl 0, ДМ ИВ. При Dl , ДМ ИП.
На настоящий момент не существует точных аналитических методов определения коэффициента диффузии. С достаточной степенью точности его можно определить только экспериментально, например, через центральный момент второго порядка кривой отклика на импульсное возмущение и критерий Пекле,
|
где
– критерий Пекле.
При значениях
,
можно принять
|
Уравнение (4.24) является основным уравнением для определения критерия Пекле по экспериментальным данным. Отклики модели на типовые возмущения представлены на рис. 4.10.
Рис. 4.10. Отклики модели на типовые возмущения
ДМ наиболее адекватно описывает структуру потока в аппаратах трубчатого и колонного типа, где наблюдается заброс вещества в сторону, противоположную направлению основного потока.