- •Введение
- •Цели и задачи моделирования:
- •Требования, предъявляемые к модели
- •Классификация моделей.
- •Виды моделирования
- •Физическое моделирование
- •1. Геометрическое подобие.
- •2. Физическое подобие.
- •Математическое моделирование
- •1. Теоретический метод
- •2. Метод аналогии.
- •Эмпирический метод
- •При определении реакции объекта на стандартное возмущение на вход подается стандартный сигнал. Это может быть одиночный импульс, ступенчатое либо синусоидальное изменение входного параметра:
- •Экспериментально-аналитический метод
- •Математическое моделирование с использованием эвм
- •Рассмотрим содержание каждого из этапов.
- •Моделирование химических процессов и аппаратов
- •Типовые математические модели структуры потоков в аппаратах
- •М одель аппарата идеального смешения
- •Начальные условия (ну).
- •Ограничения.
- •М одель аппарата идеального вытеснения.
- •Положения и определения кинетики химических реакций
- •1. Для одностадийной необратимой реакции
- •Тогда для гомогенной простой необратимой реакции вида
- •Показатели эффективности химических реакций
- •Устойчивость химических процессов
- •Оптимизация химико-технологических процессов
- •1. Составить математическую модель объекта.
- •2. Выбрать критерий оптимальности.
- •3. Установить ограничения.
- •4. Выбрать оптимизирующие факторы.
- •5. Записать целевую функцию.
- •6. Выбрать метод оптимизации.
- •Аналитический метод оптимизации
- •Оптимизация простой реакции
- •Оптимизация параллельной реакции Рассмотрим реакцию
- •Поиск оптимума численными методами
- •Методы одномерного поиска
- •Методы многомерного поиска
- •Метод покоординатного спуска.
- •Метод градиента
- •Метод случайного поиска
- •Метод динамического программирования
- •Принцип оптимальности
- •Решение комбинаторной задачи
М одель аппарата идеального вытеснения.
Условия физической реализуемости этой модели выполняются в случае поршневого потока, когда предполагается, что в направлении его движения смешивание полностью отсутствует, а в направлении, перпендикулярном движению, происходит идеальное смешение. Время пребывания всех частиц в системе одинаково и равно отношению объема системы к объемному расходу жидкости, τ=V/Q.
Уравнение, описывающее изменение концентрации в зоне идеального вытеснения, имеет вид:
,
где S – площадь поперечного сечения, м2; Z – координата длины, м; U – скорость потока, м/с.
Здесь левая часть – изменение количества вещества в выделенном элементарном слое. Первое слагаемое правой части – количество вещества, вносимого с потоком в выделенный объем, второе слагаемое – количество вещества, уносимого с потоком.
Раскрыв скобки и сократив на dz, получим:
,
где Q=US – объемный расход, м3/с.
Уравнение теплового баланса:
или после упрощения:
Для АИВ, в котором протекает химическая реакция и имеется теплообменник, уравнения материального и теплового балансов имеют вид:
,
,
где FУД = F/V – удельная поверхность теплообмена, м2/м3.
Уравнение кинетики в общем виде записывается:
,
где - вектор концентраций; {С1, С2,…Сn}, k - константа скорости.
В математическую модель АИВ входят также граничные и начальные условия и ограничения.
Граничные условия для АИВ - условия протекания процесса на входе в слой:
z=0, C(0,τ)=CВХ, T(0,τ)=TВХ
Начальные условия устанавливают значения концентраций и температуры в реакторе в начальный момент времени:
τ = 0, C(z,0) = C(0), T(z,0) = T(0)
Ограничения:
,
В случае стационарного режима:
- уравнение материального баланса:
- уравнение теплового баланса:
Поскольку τ = V/Q = S·L/Q откуда L=τ·Q/S и dz = (Q/S)·dτ, делая подстановку, получаем:
(*)
(**)
Поскольку переменная τ имеет размерность времени и обозначает время пребывания элементарного объема реагирующей смеси в зоне реакции, то уравнения (*) и (**) можно рассматривать так же как математическое описание АИС периодического действия.
Положения и определения кинетики химических реакций
1. Для одностадийной необратимой реакции
Согласно закону действующих масс скорость реакции равна константе скорости умноженной на произведение концентраций исходных веществ в степенях, равных стехиометрическим коэффициентам:
,
где k – константа скорости реакции, Сi – концентрация i-ого исходного вещества в реакции, νi – стехиометрический коэффициент, стоящий перед i-ым исходным веществом в реакции, n - количество исходных веществ в реакции.
Тогда скорость реакции, показанной для примера, будет определена как
Зависимость константы скорости k от температуры определяется законом Аррениуса:
,
где k0 – предэкспоненциальный множитель;E () – энергия активации реакции, Дж/моль; R = 8.31451 Дж/(мольК)– универсальная газовая постоянная.
По определению скорость превращения j-ого компонента ωj выражается изменением количества этого компонента в единицу времени в единице реакционного пространства. Под реакционным пространством подразумевается зона протекания реакции – объем V или поверхность контакта F.
Тогда
скорость превращения j-ого
компонента в замкнутой системе при
протекании гомогенной реакции:
,
при
протекании гетерогенной реакции:
.
С учетом того, что Сj = nj/V получим:
.
В
случае если V=const
имеем:
.
С другой стороны скорость превращения j-ого компонента ωj выражается соотношением:
,
где
m
- число реакций, в которых j-ое
вещество принимает участие, νi
– стехиометрический коэффициент,
стоящий перед j-ым
веществом в i-ой
реакции, ri
– скорость i-ой
реакции, знак «+» ставится перед
произведением
,
если в i-ой
реакции j-ое
вещество получается и знак «-» ставится
перед произведением
,
если в i-ой
реакции j-ое
вещество расходуется.
