Degtyarenko_dlya_studentov_II_kursa_2013 / Детерминир. модели / Детерм. мод. - лаб. работы / Лаб. раб. 5 / Решенные лекционные задачи / Лаб. работа5.2

1.2. ПОСЛЕДОВАТЕЛЬНЫЕ РЕАКЦИИ

Случай 1.

Уравнение кинетической кривой для исходного реагента А

In[1]:= CA CA0 Exp k1 t

Out[1]= t k1 CA0

Аналитический способ нахождения кинетической кривой для промежуточного продукта B при помощи встроенной функции DSolve системы Mathematica

In[4]:= CP CA0 CA CB FullSimplify

CA0 1 t k2 k1 1 t k1 k2

Out[4]= k1 k2

Итак, получено аналитическое решение прямой кинетической задачи,

рассматриваемой в случае 1

2 Лаб. работа5.2.1.nb

Представим графически полученные для случая 1

результаты. Приведем графики кинетических кривых для всех участников реакции при определенных числовых значениях констант скоростей последовательных стадий и начальной концентрации исходного реагента А

In[5]:= CA0 1; k1 0.5; k2 1 3;

In[6]:= Needs "PlotLegends`"

In[7]:= Plot CA, CB, CP , t, 0, 10 , PlotRange 0, Full ,

PlotStyle Thickness 0.01 , Dashed, Thickness 0.02 , RGBColor 1, 0, 0 ,

DotDashed, Thickness 0.02 , RGBColor 0, 1, 0 , PlotLabel "Решение прямой кинетической задачи для последовательной реакции первого порядка

в случае 1", AxesLabel "t", "С" , PlotLegend "СA t ", "СB t ", "СP t " , LegendPosition 1.3, 0.5 , GridLines Automatic

Лаб. работа5.2.1.nb 3

Найдем в случае 1 численное решение исходной системы

ДУ с системой соответствующих начальных условий, если t 0;10 ,

при указанных выше числовых значениях констант скоростей последовательных стадий и начальных концентраций реагирующих веществ. Сравним графики численного решения и решения, полученного аналитическим путем.

In[8]:= CA .; CB .; CP .

In[9]:= sol NDSolve CA ' t k1 CA t , CB ' t k1 CA t k2 CB t ,

CP ' t k2 CB t , CA 0 CA0, CB 0 0, CP 0 0 , CA, CB, CP , t, 0, 10

Out[9]= CA InterpolatingFunction 0., 10. , ,

CB InterpolatingFunction 0., 10. , , CP InterpolatingFunction 0., 10. ,

In[10]:= Plot CA t . sol 1 , CB t . sol 1 , CP t . sol 1 ,

t, 0, 10 , PlotLabel "График численного решения

математической модели из случая 1", AxesLabel "t", "С" , PlotRange 0, Full , PlotStyle Thick, Thickness 0.01 , Dashed, Thickness 0.02 , RGBColor 1, 0, 0 ,

DotDashed, Thickness 0.02 , RGBColor 0, 1, 0 ,

PlotLegend "СA t ", "СB t ", "СP t " , LegendPosition 1.3, 0.5 ,

GridLines Automatic

Результаты аналитического и численного решения идентичны

Глядя на полученные графики, можно предположить, что для

кинетической кривой CB t характерно наличие точки локального максимума,

адля кинетической кривой CP t наличие точки перегиба графика этой функции. Выясним,

действительно ли это так и, если да, то что это за точки.

Исследуем сначала вопрос о существовании на

промежутке 0; точек локального экстремума функции СB t

4 Лаб. работа5.2.1.nb

In[11]:= CA0 .; k1 .; k2 .;

In[12]:= CB 3

t k1 1 t k1 k2 CA0 k1

Out[12]=

k1 k2

In[13]:= d1 D CB, t FullSimplify

t k1 CA0 k1 k1 t k1 k2 k2

Out[13]=

k1 k2

Критических точек на 0; функция СB t не имеет. Исследуем вопрос

о существовании стационарных точек этой функции на указанном промежутке

In[14]:= Reduce k1 t k1 k2 k2 0, k1 k2, k1 0, k2 0 , t, Reals

Log k1

k2

Out[14]= k1 0&& 0 k2 k1 k2 k1 &&t k1 k2

Log k1

k2

In[15]:= tst

k1 k2

Log k1

k2

Out[15]=

k1 k2

Cтационарная точка функции СB t на промежутке 0; существует и единственна,