В случае достаточно простой схемы реакции, когда возможно получение функции c(t) в явном виде, проверка адекватности сводится к вычислению теоретических значений функции при различных значениях времени и сравнении их с экспериментальными. К сожалению, часто системы дифференциальных уравнений, описывающие сложные химические реакции, не могут быть проинтегрированы в аналитическом виде, и приходится прибегать к численным методам интегрирования.

Для решения прямой кинетической задачи, т.е. по данному механизму, набору констант скорости k и начальных концентраций с0 рассчитать значения концентраций c(t) в любой конкретный момент времени, нами использовался метод Рунге-Кутта четвертого порядка (программа MODEL) — один из мощных одношаговых методов, легко обобщающийся на системы уравнений: yi+1=yi+([pic])(k1+2k2+2k3+k4), (2.3) где k1=f(xi,yi), k2=f(xi+h/2; yi+hk1/2), k3=f(xi+h/2; yi+hk2/2), k4=f(xi+h; yi+hk3).

Для получения решения с требуемой степенью точности использовалось последовательное дробление шага вычислений вдвое и алгоритм последовательного прохода экспериментальных точек.

Однако, как правило, значения констант скоростей отдельных стадий неизвестны, и для их нахождения снова приходится решать оптимизационную задачу:

Функционал (2.4) относится к числу наиболее сложных для оптимизации функций — так называемых «овражных». Очевидно, что константы скорости более медленных стадий будут являться существенными переменными, определяющими размерность «оврага» в окрестностях минимума.

Кроме того, функция f(k, c0, t) сама получается в результате оптимизационной процедуры, что делает ее значения «дорогими» в смысле вычислительных затрат, а значения частных производных (f/(ki недоступны в аналитическом виде и весьма труднодоступны для получения методом конечных разностей.

Поэтому для программы “КИНО”, решающей обратную задачу химической кинетики (задачу параметрической идентификации), был выбран комбинированный метод нулевого порядка (КМНП), не требующий вычисления частных производных и описанный в [42]. Его применение позволяет решать задачу параметрической идентификации, т.е. находить константы скорости отдельных стадий многостадийной химической реакции.

Для получения формальных кинетических уравнений использовался метод варьирования концентрации одного из реагентов при постоянстве всех остальных параметров. При этом определяемая из полулогарифмического представления кинетической кривой константа скорости псевдопервого порядка линейно зависит от варьируемой концентрации. Тангенс угла наклона этой прямой соответствует частному порядку реакции по выбранному реагенту. В случае, когда вся серия измерений проводилась в режиме измерения оптической плотности при одних и тех же параметрах АЦП, можно проводить обработку значений UT, так как они прямо пропорциональны значениям оптической плотности, а следовательно, концентрациям реагирующих веществ в каждый момент времени.

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