Субботина Ю.А.
Практическая работа №10
"Решение обратной задачи химической кинетики
многостадийной химической реакции."
Цель работы:
Изучить решение обратной задачи химической кинетики многостадийной химической
реакции с использованием численных методов решения дифференциальных уравнений и
подпрограммы «Поиск решения».

Постановка обратной задачи химической кинетики:
Исходными данными для обратной задачи химической кинетики являются: кинетическая
схема реакции, начальные концентрации всех веществ, концентрации одного или
нескольких веществ в различные моменты времени (принадлежащие заданному
временному интервалу).
Необходимо определить константы скорости всех элементарных реакций.

Методы решения обратной задачи химической кинетики:
Решение обратной задачи химической кинетики может быть сведено к многократному
решению прямой задачи:
1) задаются начальные значения констант скорости элементарных стадий реакции;
2) решается прямая задача химической кинетики с выбранными константами скорости;
3) вычисляется сумма квадратов отклонений расчетных концентраций от экспериментальных;
4) полученная сумма квадратов отклонений минимизируется за счет изменения констант
скорости.
Вывод:
В данной работе нашли обновлённые коэффициенты на основе эксперементальных данных.
Для выполнения данной работы были использованы формулы исправленного метода Эйлера:
Значение функций y(x), z(x) в следующей точке рассчитывается на основе предыдущих
значений этих функций и их производных, рассчитанных в начале и конце отрезка (точках
xi и xi+h):
yi+1 = yi + hФ1i
zi+1 = zi + hФ2i, где
K1= f(xi, yi, zi) = значение производной функции f в предыдущей точке
L1= g(xi, yi, zi) = значение производной функции g в предыдущей точке
yэйл = yi + h K1 = значение функции f, рассчитанное простым методом Эйлера
zэйл = zi + h L1 = значение функции g, рассчитанное простым методом Эйлера
К2= f(xi+h, yэйл,zэйл) = значение производной функции f на конце отрезка
L2= f(xi+h, zэйл,zэйл) = значение производной функции g на конце отрезка
Ф1i= 0,5 (K1+ K2) = среднее значение производной функции f
Ф2i= 0,5 (L1+ L2) = среднее значение производной функции g
This site was made on Tilda — a website builder that helps to create a website without any code
Create a website