Методы анализа нелинейных динамических моделей. - Холодниок М.
Скачать (прямая ссылка):
4.3.1.1. Задача 11. Система "реакция - диффузия" для кинетики типа
"брюсселятор"
Если использовать модельную кинетику типа "брюсселятор"* (см. задачу 7),
то функции / и g в уравнениях (4.3.7) принимают вид
f{x, у) = А - (В+1)х + х2у, (РП-1)
g (х, у) = Вх - х2у.
Тривиальное стационарное решение (4.3.11) в этом случае записывается как
х(г)*=А, y(z) = B/A. (Р11-2)
Роль параметров в этой задаче играют Dx, Dy, L, А, В.
4.3.1.2. Задача 12. Система "реакция - диффузия", случай SW-кинетики
Если рассматривать распределенную систему с кинетикой, описываемой SH-
моделью (см. задачу 4), то функции fug в уравнениях (4.3.7)
представляются в виде
f{x, y) = a(v0 + ху)/{\ +хЧ) - х(1 А-У), (Р12-1)
g(x, у) = х (Р -f у) - дх. (Р12-2)-
Решение системы уравнений f - g = 0, в отличие от модели типа
"брюсселятор", здесь может быть найдено только численно (при заданных
значениях параметров а, б, vo > 0, р, у >• > 1). При этом в определенном
диапазоне изменения параметров можно получить несколько решений. Как и
выше, задачу 12 можно рассматривать с граничными условиями всех трех
типов. Отметим, что данная задача имеет восемь параметров: Dx, Dy> L, a,
6, v0, P, у.
4.3. Задачи с распределенными параметрами
115
¦4.3.1.3. Задача 13. Система "реакция - диффузия" в случае модели
Майнхардта
Одной из наиболее известных моделей морфогенеза в настоящее время
является модель системы "реакция-диффузия" с кинетикой типа активатор-
ингибитор, предложенная Майн-хардтом [4.44-4.47]. Эта модель описывает
пространственную дифференциацию ткани в процессе эмбриогенеза при почти
симметричных начальных условиях.
Модель описывается системой двух уравнений "реакция- .диффузия" типа
(4.3.7). Функции f и g в данном случае имеют (r)ид
f (х, у) = РоР + срх2/у - рх, (Р13-1)
g (х, у) = с'р'х2 - vy. (Р13-2)
Величины р и р' характеризуют здесь плотность источников для веществ,
действующих как активатор и ингибитор, а ро, с, р,
¦с', v - положительные параметры. У этой системы существует
единственное тривиальное решение х, у:
х = РоР/М- + cpv/c'p'p, (Р13-3)
у = c'p'x2/v. (Р13-4)
Для задачи 13 можно использовать граничные условия всех трех типов - ГУ1,
ГУ2 или ГУЗ. В данной задаче имеется 10 параметров: Dx, Dy, L, р, р', р0,
с, р, с', v.
4.3.2. Задача 14. Трубчатый неизотермический реактор с аксиальным
перемешиванием
Модель трубчатого реактора с аксиальным перемешиванием используется как
для гомогенных, так и для гетерогенных (каталитических) реакторов
([4.48], [4.49]). В последнем случае полезна упрощенная, псевдогомогенная
модель, основанная на предположении, что гетерогенную систему
"катализатор-реакционная смесь" можно заменить гомогенной средой с некими
эффективными характеристиками. Такая псевдогомогенная модель и
формулируется в задаче 14.
Рассмотрим трубчатый реактор с теплопередачей через стенку (см. рис.
4.10). Будем считать, что реакционная смесь полностью перемешивается в
радиальном направлении; тем самым мы будем рассматривать только
продольные градиенты концентраций компонент и температуры. Далее,
предположим, что плотность потока компонент в продольном направлении
описывается соотношением вида ~-De{dc/dt), аналогичным
116
Глава 4
закону Фика, где De - эффективный коэффициент диффузии* с - концентрация
компоненты.
Аналогично будем предполагать, что плотность потока тепла в продольном
направлении задается соотношением вида -ke(dT/dl), где ke - эффективный
коэффициент теплопроводности.
Рис. 4.10. Трубчатый реактор с аксиальным переносом тепла и массы;
граничные условия типа Данквертса.
Предположим, далее, что в реакторе протекает реакция первого порядка,
описываемая выражением для скорости реакции вида kooCexpi-E/RT) (Е -
энергия активации, R - газовая постоянная) и с тепловым эффектом
(энтальпией реакции) -АЯг. Мы будем считать, что плотность теплового
потока через стенку трубки определяется выражением 4U(T - Tc)/d, где U -
соответствующий коэффициент теплопередачи, d - диаметр трубки и Тс -
температура теплообменной (например, охлаждающей) среды вне трубки.
Мы будем предполагать, что в реакторе имеется катализатор. Пусть
плотность жидкости равна pt, а ее теплоемкость на единицу объема при
постоянном давлении постоянна и равна Срг, плотность же катализатора
равна ps, а его теплоемкость на единицу объема также постоянна и равна
Cps. Тогда уравнения баланса массы и энергии можно представить в виде
= (-E/RT), (Р14-1)
(Г -Ус) + (1 -ерД-AHJk^cexpi-E/RT). (Р14-2>
z = 0 L =0
z - 1 I = L
Здесь т - время, eP - доля объема, занятая жидкостью, в v - (постоянная)
скорость жидкости. Начальные условия еле-
4.3. Задачи с распределенными параметрами
117"
дующие:
т = 0: c(t, 0) = c°(l), T(l, 0) = Т°(1). (Р14-3)
Пусть теперь на границе области выполняются граничные условия типа
Данквертса (рис. 4.10) (см. [4.8])
т > 0, 1 = 0: -Д. - = у (с0 - с),
-K-w = ve <САТе-Т)-
(Р14-4)-
т > 0, l = L: = = (PI4-5)
Здесь Со и То - концентрация и температура до входа в реактор. Введем
