Научная статья на тему 'НЕЛИНЕЙНЫЕ КИНЕТИЧЕСКИЕ ЗАКОНЫ СОХРАНЕНИЯ В ЗАКРЫТОМ БЕЗГРАДИЕНТНОМ РЕАКТОРЕ'

НЕЛИНЕЙНЫЕ КИНЕТИЧЕСКИЕ ЗАКОНЫ СОХРАНЕНИЯ В ЗАКРЫТОМ БЕЗГРАДИЕНТНОМ РЕАКТОРЕ Текст научной статьи по специальности «Химические науки»

CC BY
34
5
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ХИМИЧЕСКАЯ КИНЕТИКА / НЕЛИНЕЙНЫЕ РЕАКЦИИ / ЗАКРЫТЫЙ ИЗОТЕРМИЧЕСКИЙ РЕАКТОР ИДЕАЛЬНОГО СМЕШЕНИЯ / НЕЛИНЕЙНЫЕ КИНЕТИЧЕСКИЕ ЗАКОНЫ СОХРАНЕНИЯ / ОБРАТНАЯ ЗАДАЧА

Аннотация научной статьи по химическим наукам, автор научной работы — Кольцов Николай Иванович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

NONLINEAR KINETIC CONSERVATION LAWS IN CLOSED GRADIENTLESS REACTOR

A simple method has been developed that allows one to a priori determine the type, calculate the conditions of existence, and discover new nonlinear kinetic conservation laws for chemical reactions proceeding according to linear and nonlinear mechanisms within the framework of the law of acting masses in a closed isothermal ideal mixing reactor. Such nonlinear kinetic conservation laws are a kind of multiexperimental autonomous kinetic invariants that can be observed on the basis of data measured in one nonequilibrium experiment with any initial conditions given in advance. The method is applicable to reactions involving free reagents, i.e. reagents not interacting with other reagents. Such reagents, as a rule, are the starting materials involved in the irreversible stages. Stoichiometric conditions (criteria) for the existence of reactions involving free reagents are obtained. It is shown that this criterion is always fulfilled for reactions whose mechanisms necessarily include irreversible stages of expenditure of the starting materials, but do not contain the stages of their formation. Examples of possible reaction mechanisms involving free reagents are given. An effective algorithm has been developed to obtain an exact solution to the inverse problem of chemical kinetics for reactions involving free reagents without the use of complex optimization methods. Examples of the implementation of the developed method for specific model oxidation reactions of carbon, carbon monoxide and hydrogen are presented, illustrating the possibility of establishing new nonlinear kinetic conservation laws and solving the inverse problem of reconstructing the rate constants of the stages of linear and nonlinear reactions in a closed isothermal ideal mixing reactor.

Текст научной работы на тему «НЕЛИНЕЙНЫЕ КИНЕТИЧЕСКИЕ ЗАКОНЫ СОХРАНЕНИЯ В ЗАКРЫТОМ БЕЗГРАДИЕНТНОМ РЕАКТОРЕ»

Т 65 (1)

ИЗВЕСТИЯ ВЫСШИХ УЧЕБНЫХ ЗАВЕДЕНИЙ. Серия «ХИМИЯ И ХИМИЧЕСКАЯ ТЕХНОЛОГИЯ»

2022

V 65 (1) ChemChemTech 2022

DOI: 10.6060/ivkkt.20226501.6268 УДК: 541.12

НЕЛИНЕЙНЫЕ КИНЕТИЧЕСКИЕ ЗАКОНЫ СОХРАНЕНИЯ В ЗАКРЫТОМ БЕЗГРАДИЕНТНОМ РЕАКТОРЕ

Н.И. Кольцов

Николай Иванович Кольцов (ORCID 0000-0003-2264-1370)*

Кафедра физической химии и высокомолекулярных соединений, Чувашский государственный университет им. И.Н. Ульянова, Московский пр., 15, Чебоксары, Российская Федерация, 428015 E-mail: koltsovni@mail.ru*

Разработан простой метод, позволяющий априори определять вид, рассчитывать условия существования и обнаруживать новые нелинейные кинетические законы сохранения для химических реакций, протекающих по линейным и нелинейным механизмам в рамках закона действующих масс в закрытом изотермическом реакторе идеального смешения. Такие нелинейные кинетические законы сохранения представляют собой разновидность муль-тиэкспериментных автономных кинетических инвариантов, которые можно наблюдать на основе данных, измеренных в одном неравновесном эксперименте с любыми, наперед заданными начальными условиями. Метод применим к реакциям, протекающим с участием свободных реагентов, т.е. реагентов, не взаимодействующих с другими реагентами. Такими реагентами, как правило, являются исходные вещества, участвующие в необратимых стадиях. Получены стехиометрические условия (критерий) существования реакций с участием свободных реагентов. Показано, что этот критерий всегда выполняется для реакций, механизмы которых обязательно включают необратимые стадии расходования исходных веществ, но не содержат стадии их образования. Приведены примеры возможных механизмов реакций с участием свободных реагентов. Разработан эффективный алгоритм получения точного решения обратной задачи химической кинетики для реакций с участием свободных реагентов без применения сложных оптимизационных методов. Приведены примеры реализации разработанного метода для конкретных модельных реакций окисления углерода, монооксида углерода и водорода, иллюстрирующие возможность установления новых нелинейных кинетических законов сохранения и решения обратной задачи восстановления констант скоростей стадий линейных и нелинейных реакций в закрытом изотермическом реакторе идеального смешения.

Ключевые слова: химическая кинетика, нелинейные реакции, закрытый изотермический реактор идеального смешения, нелинейные кинетические законы сохранения, обратная задача

NONLINEAR KINETIC CONSERVATION LAWS IN CLOSED GRADIENTLESS REACTOR

N.I. Kol'tsov

Nikolay I. Kol'tsov (ORCID 0000-0003-2264-1370)*

Department of Physical Chemistry and Macromolecular Compounds, Chuvash State University named after I.N. Ulyanov, Moskovskiy ave., 15, Cheboksary, 428015, Russia E-mail: koltsovni@mail.ru*

A simple method has been developed that allows one to a priori determine the type, calculate the conditions of existence, and discover new nonlinear kinetic conservation laws for chemical reactions proceeding according to linear and nonlinear mechanisms within the framework of the

law of acting masses in a closed isothermal ideal mixing reactor. Such nonlinear kinetic conservation laws are a kind of multiexperimental autonomous kinetic invariants that can be observed on the basis of data measured in one nonequilibrium experiment with any initial conditions given in advance. The method is applicable to reactions involving free reagents, i.e. reagents not interacting with other reagents. Such reagents, as a rule, are the starting materials involved in the irreversible stages. Stoichiometric conditions (criteria) for the existence of reactions involving free reagents are obtained. It is shown that this criterion is always fulfilled for reactions whose mechanisms necessarily include irreversible stages of expenditure of the starting materials, but do not contain the stages of their formation. Examples of possible reaction mechanisms involving free reagents are given. An effective algorithm has been developed to obtain an exact solution to the inverse problem of chemical kinetics for reactions involving free reagents without the use of complex optimization methods. Examples of the implementation of the developed method for specific model oxidation reactions of carbon, carbon monoxide and hydrogen are presented, illustrating the possibility of establishing new nonlinear kinetic conservation laws and solving the inverse problem of reconstructing the rate constants of the stages of linear and nonlinear reactions in a closed isothermal ideal mixing reactor.

Key words: chemical kinetics, nonlinear reactions, closed isothermal reactor of ideal mixing, nonlinear kinetic conservation laws, inverse problem

Для цитирования:

Кольцов Н.И. Нелинейные кинетические законы сохранения в закрытом безградиентном реакторе. Изв. вузов. Химия и хим. технология. 2022. Т. 65. Вып. 1. С. 23-29

For citation:

Kol'tsov N.I. Nonlinear kinetic conservation laws in closed gradientless reactor. ChemChemTech [Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol.]. 2022. V. 65. N 1. P. 23-29

ВВЕДЕНИЕ

Кинетические законы сохранения (КЗС) можно разделить на два принципиально различных типа - линейные (ЛКЗС) и нелинейные (НКЗС) [1-5]. Соответственно, они представляют собой не зависящие от времени (автономные) комбинации неравновесных значений концентраций реагентов и кинетических параметров реакции в виде линейных и нелинейных алгебраических уравнений. Проведенные в работах М.Д. Корзухина, Г.С. Яблонского и др. [1-3] исследования позволили установить все возможные ЛКЗС для линейных и нелинейных реакций, протекающих в изотермическом безградиентном закрытом реакторе идеального смешения (РИС). НКЗС исследованы меньше и пока обнаружены только для линейных и отдельных специфических по форме (редко встречающихся на практике) нелинейных реакций, поэтому точное их число в закрытых системах не известно [4-5]. В связи с этим представляет интерес разработать метод, позволяющий устанавливать НКЗС для более широких классов распространенных на практике механизмов реакций. Ниже приведен метод, с помощью которого можно находить НКЗС в сложных многостадийных нелинейных химических реакциях, протекающих в закрытом изотермическом

безградиентном РИС. Метод базируется на оригинальных идеях использования в химической практике пар экспериментов с разными начальными условиями, впервые предложенных Г.С. Яблонским с соавторами [6-10]. Обобщение этих идей для более широкого круга экспериментов изложено в [11-18]. Отличие предложенного метода состоит в использовании свойств реагентов, не взаимодействующих с другими реагентами. Наличие в механизме реакции пар таких реагентов позволяет находить новые НКЗС. Следует отметить, что все законы сохранения, ЛКЗС и НКЗС, могут использоваться для решения обратной задачи (ОЗ) химической кинетики без применения оптимизационных методов, что ранее в литературе не отмечалась, см. например [19-29]. Это естественным образом следует из структуры законов сохранения, связывающих кинетические параметры с концентрациями реагентов химических реакций. Возможности использования НКЗС для решения ОЗ приведены ниже.

ТЕОРЕТИЧЕСКАЯ ЧАСТЬ

Рассмотрим механизмы сложных химических реакций с участием любого числа реагентов А] и протекающих через нелинейные стадии вида

ЪaijАj = Ъa-ijА], i=1,...,s,]=1,...,п, (1)

N.1. КО^ОУ

где щ, а-у — стехиометрические коэффициенты реагентов в левых и правых частях стадий. Выделим среди этих механизмов такие, в которых участвуют как минимум два «свободных» (назовем их так) реагента А^Ли и В=Ай, не взаимодействующие с другими реагентами (однобуквенные обозначения реагентов будем использовать для краткости и далее). Тогда механизм реакции должен включать хотя бы две стадии «без взаимодействий» разных исходных веществ вида А^2С, 2B^D+E и т.д., причем свободные реагенты не должны взаимодействовать с другими реагентами и в других стадиях. Правые части этих стадий могут быть любыми, т.к. не влияют на скорость образования исходных веществ. Динамика таких реакций в закрытом изотермическом РИС в рамках закона действующих масс описывается системой обыкновенных дифференциальных уравнений (ОДУ) [2, 5, 30]:

Л '= Е(а— —аЩ Л/——П, Л—, 7=1,... ,п, (2)

где А] — концентрации реагентов А;, мол. доли; Ло7 — начальные условия (н.у.); кг, к— — константы скоростей стадий в прямом и обратном направлениях, 1/с. Выразим критерий существования свободных реагентов в механизме реакции (1) через стехиометрию и кинетические параметры

а,*к Ф 0, а*к а* = 0, к—* = 0, (3)

где 1*=11,12,... — номера стадий с участием свободных реагентов. При выполнении критерия (3) скорость образования свободных реагентов зависит только от их концентрации Лк и констант скоростей стадий с их участием. Соответствующие уравнения (2) для двух свободных реагентов примут вид Л '= Zг*(a-г■*k-aг■*k)kг■*П^■kЛaг к = =/л(а,*к,а—*к ,Л,Ло,кг*), (4)

В ' = 2г*(а—г*к—*к)кг *П./кВЖ к = =/в(аг*к,а—г*к ,В,ВоМ, (5)

где /л,/в — полиномы, степени которых определяются максимальной молекулярностью стадий, как правило, не превышающей трех. Уравнения (4)-(5) не зависят от концентраций других реагентов. Поэтому для конкретных реакций они могут быть решены точно, но в общем случае их решения можно записать только в интегральном виде

Л(0 = / /л Ж = Ел^Щ+Сл^ъа—г'к ,Ло,к*), (6) В(0 = / /в Ж = Ев((,кг*)+Св(аг*к,а—г*к ,Во,кг*), (7) где I - время, с; Рл, Рв - первообразная; Сл, Св — константы интегрирования. Выразим из любого из уравнений (6)-(7) время

^ = т(а1*к,а—*к ,Ло,Во,кг*). (8)

Подставим его в другое, т.е. исключим время из уравнений (6)-(7), и получим автономный НКЗС для пары свободных реагентов А, В:

Елв(А],аг*к,а—г*к ,Во,кг*) = 0. (9)

Перенесем в (9) все слагаемые, содержащие переменные величины (концентрации реагентов), влево, а остальные слагаемые (константы, н.у. и др.) вправо. Получим удобную эквивалентную (симметричную) запись НКЗС

Клв = Рлв(А],кг*) = Р°лв(Ло,Во,кг*). (10) Аналогичные НКЗС можно получить для всех других пар свободных реагентов, если их несколько. Каждое из соотношений (9)-(10) зависит от концентраций реагентов, кинетических параметров реакции и н.у., но не зависит от времени (автономно). Следовательно, эти зависимости могут наблюдаться экспериментально при физичных значениях входящих в них параметров для любых реакций с участием двух и более свободных реагентов. Число таких НКЗС, очевидно, не меньше числа комбинаций по два из всех реагентов. Например, механизм реакции A+C=B+D+E с участием двух свободных реагентов А и С может иметь вид Л^2В, С^2Б, В+Б=Е, для которого можно определить НКЗС с помощью соотношений (10).

Опишем кратко, как для стадий с участием свободных реагентов на основе неравновесных данных можно решить и обратную задачу (ОЗ) без использования оптимизационных методов. Так, если подставить начальные концентрации свободных реагентов в (9) или (10) и решить полученные уравнения, то определяются комплексы констант соответствующих стадий. Чтобы найти каждую константу в отдельности, необходимо подставить те же начальные концентрации свободных реагентов в равенства (6)-(7) и найти все решения полученных пар систем нелинейных (но, что существенно, независимых) уравнений, методы решения которых хорошо разработаны, см. например, [31]. Для получения точных физичных (положительных) решений достаточно выбрать любой момент времени (кроме начального и конечного-равновесного) и «измерить» соответствующие концентрации исходных веществ. Для оценки погрешности измерений нужно поварьировать «измеренные» концентрации в пределах погрешности измерительных приборов (5-20%) и вновь, но уже приближенно, решить ОЗ.

Пример 1. Рассмотрим брутто-реакцию Л+С= Е протекающую по схеме

1) Л = В, 2) 2С = 2Б, 3) В + Б = Е. (1.1) Например, реакция С+О2 = СО2 может протекать через стадии 1) С=С*, 2) О2=2О*,

3) С*+2О*=СО2. Схема вида (1.1) описывается системой ОДУ

Л '= —кЛ + к—В, (1.2)

В' = Ы -к-В -кзВВ + к-зЕ, (1.3)

С' = -2к2С2 + 2к-2В?, (1.4)

В' = 2к2С2 - 2к-2&-кзВВ + к-Е, (1.5)

Е'= кзВВ -к-зЕ. (1.6)

В этой реакции участвуют пять реагентов, три из которых независимы. Следовательно, согласно [3, 5], данная реакция имеет два линейных стехиометрических закона сохранения (ЛСЗС) - их нетрудно выписать, но не имеет ЛКЗС. Покажем, что при этом существуют НКЗС и найдем их вид. Из (1.2)-(1.6) следует, что если в реакции (1.1) стадии 1 и 2 необратимы к- = к- = 0, то выполняются условия (3), и скорость расходования исходных веществ А и С не зависит от концентраций остальных реагентов. Тогда исходные реагенты становятся свободными, а соответствующие им ОДУ (1.2) и (1.4) принимают вид (4), т.е. становятся независимыми и интегрируются

А = Лоо+ (Ао -Лоэ)ехр(-к^), (1.7)

С = Со+ 11(2к21+1/Со), (1.8)

где А о = Со = 0 - равновесные концентрации реагентов А и С. Из (1.8) следует

Х = (1/С-1/Со)/(2к2). (1.9)

Подставим (1.9) в (1.7) и найдем точный НКЗС А = Аоехр[(-к1(1/С-1/Со)/(2к2)], (1.10) который удобно переписать в симметричном виде

К = 1пА+к1/(2к2С) = ¡пЛо+Ы^кСо), (1.11) Отсюда, например, при к1 = 1, к2 = 1/2, к- = к-2 = 0 и н.у. Ао = Со = 1/2, Во = Во = Ео = 0 получим К = 1пА+1/С = 1п(1/2)+2 и т.д.

Решим обратную задачу. Для этого подставим соответствующие н.у. Ао = Со = 1/2 в (1.10) или (1.11) и получим одно уравнение 1пА+к;/(2к2С) = 1пАо+к;/к2, откуда находим комплекс констант к1/к2 = 1п(А/Ао)/(1-1/2С). Чтобы найти каждую константу в отдельности, подставим н.у. в (1.7) -(1.8) и получим систему из двух уравнений А = Аоехр(-ЬХ), С = 1/(2к2Х+1/Со), откуда находим к: = 1п(А</А)/Х, к2 = (1/С-1/Со)/(2Х). (1.12) Эти константы принимают положительные значения при Ао > А, Со > С. Выберем любой момент времени, отвечающий этим условиям, например Х* = 4, и найдем соответствующие концентрации исходных веществ (точные) А(Х*) « 0,0092, С(1;*) « 0,1667. Подставим эти значения в (1.12) и получим точные значения к1 = 1,00, к2 = 0,50. Оценим влияние погрешности измерений, считая, что измеренные значения найдены с 10%-ошибкой А(Х*) » 0,01, С(Х*) » 0,17. Подставим их в (1.12) и получим приближенные значения к1 « 0,9780, к2 ~ 0,4853. Увеличим погрешность измерений до 10%, считая А(Х*) » 0,011, С(Х*) » 0,2. Подставим их

в (1.12) и получим более грубые значения к1» 0,9542, к2 » 0,375 и т.д.

Пример 2. Рассмотрим реакцию А+Б = С+Б, протекающую по схеме

А = В + С, 2) D = 2Е, 3) В + 2E=F. (2.1) Примером служит реакция 2С0+02=2С02 протекающая по схеме 1) СО = С* + О*, 2) О2 = 20*, 3) С* + 20* = СО2. Такие реакции описываются ОДУ

А'= -к ¡А + к-гВС, (2.2)

В'= кА - к-ВС - кзВЕ2 + к-зР, (2.3) С'= кА -к-ВС, В' = -к2Б + к-Е2, (2.4) Е' = 2к2В - 2к-Е2 - 2кзВЕ2 + 2к-зР, ¥'= кзВЕ2 -к-Р. (2.5)

В этой реакции шесть реагентов, три из которых независимые. Следовательно, в ней есть три ЛСЗС, но нет ЛКЗС. Найдем НКЗС. Из (2.2)-(2.5) следует, что при к- = к-2 = 0 выполняются условия (3) и появляются свободные реагенты А и D, а соответствующие им ОДУ интегрируются

А = Ао+ (Ао - АО)ехр(-к1Х), (2.6)

В = Во+ (Во -Во)ехр(-к2Х), (2.7)

где Ао= Во= 0. Выразим из (2.7) время, поставим в

(2.6) и найдем НКЗС

к21п(А/Ао) = кЛп(В/Во), (2.8)

который удобно переписать в симметричном виде К = Ак2/Вк1 = Аок2/Вок1 . (2.9)

Отсюда при к1 = 1, к2 = 1/2, к- = к-2 = 0 и н.у. Ао=Во = 1/2 получим К = А1/2/В = Ао1/2/Во = 21/2 и т.д.

Решим обратную задачу. Подставим н.у. А0 = Во = 1/2 в (2.8) и получим уравнение к21п(А/Ао) = кЛп(В/Во), откуда находим комплекс констант к1/к2 = 1п(А/Ао)/1п(В/Во). Подставим н.у. в (2.6)-

(2.7) и получим систему из двух уравнений А = Аоехр(-к;Х), В = Воехр(-к2Х). Отсюда

кг = 1п(Ао/А)/Х, к2 = 1п(Во/В)/Х. (2.10)

Эти константы принимают положительные значения при Ао > А, Во > В. Выберем любой момент времени, например Х* = 4, и найдем соответствующие концентрации исходных веществ А(Х*) « 0,0061, В(Х*) « 0,0902. Подставим эти значения в (2.10) и получим практически точные значения к1 ~ 1,00, к2 ~ 0,50. Оценим влияние погрешности измерений, считая, что измеренные значения найдены с 10%-ошибкой А(Х*) » 0,01, В(Х*) » 0,10. Подставим их в (2.10) и получим приближенные значения к1« 0,8766, к2 « 0,4743 с ошибкой 6% и т.д.

Пример 3. Реакция Н2 + 1/202 = Н2О протекающая по схеме 1) Н2 = 2Н*, 2) 02 = 20*, 3) 2Н* + 0* = Ш0 (3.1) описывается системой ОДУ А' = -аА + к-1В2, В' = -к2В+к-2Е2, (3.2)

C' = k3D2E -к-зС, D' = = 2kiA - 2k-iD2 - 2k3D2E + 2к-зС, (3.3) E' = 2k2B - 2k-2E2 - k3D2E + к-зС, (3.4) где A = [H2], B = [O2], С = [H2O], D = [H*], E = [O*] -концентрации реагентов. В этой реакции пять реагентов, три из которых независимые. Следовательно, в этой реакции есть два ЛСЗС, но нет ЛКЗС. Найдем НКЗС. Из (3.2)-(3.4) следует, что при k-1 = k-2 = 0 выполняются условия (3) и появляются свободные реагенты H2 и O2, а соответствующие им ОДУ (3.2) интегрируются

A = A«+ (Ao -A«)exp(-kit), (3.5)

B = B«+ (Bo -Ba)exp(-k2t), (3.6)

где A «=6^=0. Выразим из (3.6) время, поставим в (3.5) и найдем НКЗС

k2ln(A/Ao) = kiln(B/Bo), (3.7)

который удобно переписать в симметричном виде K =Ak2/Bki = Aok2/Bok1. (3.8)

Отсюда при ki = 1, k2 = 1/2, k-i = k-2 = 0 и н.у. Ao = Bo = 1/2 получим K=Al12/B = Ao1/2/Bo = 21/2 и т.д.

Решим обратную задачу. Подставим н.у. Ao = Bo = 1/2 в (3.7) и получим уравнение k2ln(A/Ao) = kiln(B/Bo), откуда находим комплекс констант ki/k2 = ln(A/4o)/ln(B/Bo). Подставим н.у. в (3.5)-(3.6) и получим систему из двух уравнений A = Aoexp(-kit), B = Boexp(-k2t). Отсюда

ki = ln(Ao/A)/t, k2 = ln(Bo/B)/t. (3.9)

Эти константы принимают положительные значения при Ao > A, Bo > B. Выберем любой момент времени и соответствующие концентрации исходных веществ t* = 4, A(t*) « 0,0092, B(t*) « 0,0677. Подставим эти значения в (3.9) и получим ki« 0,9999, k2 ~ 0,5000. Оценим влияние погрешности измерений, считая что измеренные значения найдены с 10%-ошибкой A(t*) » 0,01, B(t*) » 0,07. Подставим

ЛИТЕРАТУРА

1. Корзухин М.Д. Линейные законы сохранения в химической кинетике. Журн. физ. хим. 1972. Т. 46. № 7. С. 1845-1847.

2. Горбань А.Н., Быков В.И., Яблонский Г.С. Очерки о химической релаксации. Новосибирск: Наука. 1986. 320 с.

3. Алексеев Б.В., Кольцов Н.И., Федотов В.Х Линейные инварианты в формальной кинетике химических реакций. Журн. физ. хим. 1992. Т. 66. № 12. С. 3219-3224.

4. Алексеев Б.В., Кольцов Н.И., Федотов В.Х. О нелинейных законах сохранения в каталитических реакциях. Журн. физ. хим. 1988. Т. 62. № 11. С. 3069-3072.

5. Кольцов Н.И. Математическое моделирование каталитических реакций. Чебоксары: Изд-во Чувашского ун-та. 2007. 294 с.

6. Yablonsky G.S. Decoding Complexity of Chemical Reactions. Theor. Found. Chem. Eng. 2014. V. 48. N 5. P. 608-613. DOI: 10.1134/S004057951405025X.

их в (3.9) и получим приближенные значения ki «

0.9780. k2 « 0,4915 и т.д. Выберем другой момент времени (вблизи с.с.) и соответствующие концентрации t* = 10, A(t*)» 6,1806e-5, B(t*) » 0,0056. Подставим эти значения в (3.9) и вновь получим ki« 0,9997, k2 « 0,5000.

ВЫВОДЫ

Установлены новые виды точных нелинейных кинетических законов сохранения (НКЗС) для нелинейных реакций, протекающих в закрытом изотермическом реакторе идеального смешения. Механизмы таких реакций включают стадии с реагентами, скорость образования которых не зависит от концентраций других реагентов. Эти механизмы характеризуются наличием необратимых стадий превращения исходных веществ и отсутствием стадий их последующего синтеза. Такие вещества существенно отличаются от «обычных» и названы нами свободными реагентами. Для всех стадий реакций с участием таких реагентов существует возможность достаточно точного решения обратных задач химической кинетики без использования сложных оптимизационных алгоритмов, что показано на примерах конкретных реакций.

БЛАГОДАРНОСТИ

Автор выражает благодарность Федотову В.Х. за обсуждение работы.

Авторы заявляют об отсутствии конфликта интересов, требующего раскрытия в данной статье.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

The author is grateful to V.Kh. Fedotov for discussing the work.

The authors declare the absence a conflict of interest warranting disclosure in this article.

REFERENCES

1. Korzukhin M.D. Linear conservation laws in chemical kinetics. Zhurn. Fiz. Khim. 1972. V. 46. N 7. P. 1845-1847 (in Russian).

2. Gorban A.N., Bykov V I., Yablonsky G.S. Essays on chemical relaxation. Novosibirsk: Nauka. 1986. 320 p. (in Russian).

3. Alekseev B.V., Kol'tsov N.L, Fedotov V.Kh. Linear invariants in the formal kinetics of chemical reactions. Zhurn. Fiz. Khim. 1992. V. 66. N 12. P. 3219-3224 (in Russian).

4. Alekseev B.V., Kol'tsov N.L, Fedotov V.Kh. On nonlinear conservation laws in catalytic reactions. Zhurn. Fiz. Khim. 1988. V. 62. N 11. P. 3069-3072 (in Russian).

5. Kol'tsov N.L Mathematical modeling of catalytic reactions. Cheboksary: Chuvash State University. 2007. 294 p. (in Russian).

6. Yablonsky G.S. Decoding Complexity of Chemical Reactions. Theor. Found. Chem. Eng. 2014. V. 48. N 5. P. 608-613. DOI: 10.1134/S004057951405025X.

7. Branco P.D., Yablonsky G.S., Marin G.B., Constales D.

Mapping the kinetic events in a linear two-step irreversible reversible reaction mechanism. Chem. Eng. Sci. 2017. V. 158. P. 370-380. DOI: 10.1016/j.ces.2016.10.032.

8. Peng B., Yablonsky G.S., Constales D., Marin G.B., Mueh-ler M. Experimental confirmation of a new invariant for a nonlinear chemical reaction. Chem. Eng. Sci. 2018. V. 191. P. 262-267. DOI: 10.1016/j.ces.2018.06.065.

9. Branco P.D., Yablonsky G.S., Marin G.B., Constales D. New invariances of chemical reactions from Scaled Incremental Conversion (SIC). Chem. Eng. Sci. 2018. V. 184. P. 25-32. DOI: 10.1016/j.ces.2018.02.020.

10. Yablonsky G.S., Branco P.D., Marin G.B., Constales D. Conservatively perturbed equilibrium in chemical kinetics. Chem. Eng. Sci. 2019. V. 196. P. 384-390. DOI: 10.1016/j.ces.2018.11.010.

11. Федотов В.Х., Кольцов Н.И. Метод определения квазиинвариантов в кинетике химических реакций. Изв. вузов. Химия и хим. технология. 2016. Т. 59. Вып. 5. С. 72-76.

12. Федотов В.Х., Кольцов Н.И. Обобщенный метод расчета нелинейных двумерных автономных квазиинвариантов динамических систем. Вестн. технол. ун-та. 2019. Т. 22. № 1. С. 119-124.

13. Федотов В.Х., Кольцов Н.И. Квазиинварианты химических реакций c неидеальной кинетикой. Хим. физика. 2019. Т. 38. № 4. С. 23-27. DOI: 10.1134/S0207401X19040046.

14. Федотов В.Х., Кольцов Н.И. Кинетические квазиинварианты химических реакций в закрытых системах. Изв. вузов. Химия и хим. технология. 2019. Т. 62. Вып. 6. С. 47-52. DOI: 10.6060/ivkkt.20196206.5881.

15. Федотов В.Х., Кольцов Н.И. Кинетические квазиинварианты химических реакций в открытых системах. Изв. вузов. Химия и хим. технология. 2019. Т. 62. Вып. 8. С. 76-80. DOI: 10.6060/ivkkt.20196208.5891.

16. Федотов В.Х., Кольцов Н.И. Автономные кинетические инварианты линейных химических реакций. Кинетика и катализ. 2019. Т. 60. № 6. С. 756-762. DOI: 10.1134/S0453 881119060042.

17. Федотов В.Х., Кольцов Н.И., Косьянов П.М. Точные автономные инварианты химических реакций с участием двух реагентов. Хим. физика. 2020. Т. 39. № 3. С. 48-52.

18. Кольцов Н.И. Точные автономные инварианты нелинейных химических реакций. Кинетика и катализ. 2020. Т. 61. № 3. С. 812-817. DOI: 10.1134/S0023158420030143.

19. Исмагилова А.С., Спивак С.И. Обратные задачи химической кинетики, Saarbrucken: Lap Lambert Acad. Publ. 2013. 118 с.

20. Денисов A.M., Дмитриев В.И. Обратные и некорректно поставленные задачи. Вестн. Московск. ун-та. Сер. 15. Вычислит. Матем. и кибернетика. 2005. С. 23-30.

21. Tarantola A. Inverse Problem Theory and Methods for Model Parameter Estimation. SIAM: Philadelphia. 2005. 344 p. DOI: 10.1137/1.9780898717921.

22. Kaipio J., Somersalo E. Statistical and computational inverse problems. New York: Springer. 2010. 339 p.

23. Chavent Guy. Nonlinear Least Squares for Inverse problems. New York: Springer. 2010. 292 p. DOI: 10.1007/97890-481-2785-6.

24. Aster R.C., Borchers B., Thurber C.H. Parameter estimation and inverse problems. New York: Elsevier. 2013. 360 p. DOI: 10.1016/B978-0-12-385048-5.00010-0.

25. Леонов А.С. Решение некорректно поставленных обратных задач: очерк теории, практические алгоритмы и демонстрации в МАТЛАБ. М.: Либроком. 2015. 336 с.

7. Branco P.D., Yablonsky G.S., Marin G.B., Constales D.

Mapping the kinetic events in a linear two-step irreversible reversible reaction mechanism. Chem. Eng. Sci. 2017. V. 158. P. 370-380. DOI: 10.1016/j.ces.2016.10.032.

8. Peng B., Yablonsky G.S., Constales D., Marin G.B., Muehler M. Experimental confirmation of a new invariant for a non-linear chemical reaction. Chem. Eng. Sci. 2018. V. 191. P. 262-267. DOI: 10.1016/j.ces.2018.06.065.

9. Branco P.D., Yablonsky G.S., Marin G.B., Constales D. New invariances of chemical reactions from Scaled Incremental Conversion (SIC). Chem. Eng. Sci. 2018. V. 184. P. 25-32. DOI: 10.1016/j.ces.2018.02.020.

10. Yablonsky G.S., Branco P.D., Marin G.B., Constales D. Conservatively perturbed equilibrium in chemical kinetics. Chem. Eng. Sci. 2019. V. 196. P. 384-390. DOI: 10.1016/j.ces.2018.11.010.

11. Fedotov V.Kh., Kol'tsov N.I Method for the determination of quasi-invariants in the kinetics of chemical reactions. ChemChemTech [Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol.]. 2016. V. 59. N 5. P. 72-76 (in Russian).

12. Fedotov V.Kh., Kol'tsov N.I Generalized method for calculating nonlinear two-dimensional autonomous quasi-invariants of dynamical systems. Vestn. Tehnol. Un-ta. 2019. V. 22. N 1. P. 119-124 (in Russian).

13. Fedotov V.Kh., Kol'tsov N.I. Quasi-Invariants of Chemical Reactions with Nonideal Kinetics. Rus. Journ. Phys. Chem. B. 2019. V. 13. N 2. P. 262-266. DOI: 10.1134/S1990793119020155.

14. Fedotov V.Kh., Kol'tsov N.I. Kinetic quasi-invariants of chemical reactions in closed systems. ChemChemTech [Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol.]. 2019. V. 62. N 6. P. 47-52 (in Russian).

15. Fedotov V.Kh., Kol'tsov N.I. Kinetic quasi-invariants of chemical reactions in open systems. ChemChemTech [Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol.]. 2019. V. 62. N 8. P. 76-80 (in Russian). DOI: 10.6060/ivkkt. 20196208.5891.

16. Fedotov V.Kh., Kol'tsov N.I. Autonomous Kinetic Invariants of Linear Chemical Reactions. Kinet. Catal. 2019. V. 60. N 6. P. 776-782. DOI: 10.1134/S002315841906003X.

17. Fedotov V.Kh., Kol'tsov N.I., Kosianov P.M. Exact autonomous invariants of chemical reactions involving two reagents. Rus. J. Phys. Chem. B. 2020. V. 14. N 2. P. 284-289. DOI: 10.1134/S 1990793120020049.

18. Kol'tsov N.I. Exact autonomous invariants of nonlinear chemical reactions. Kinet. Сatal. 2020. V. 61. N 3. P. 812-817 (in Russian). DOI: 10.1134/S0023158420030143.

19. Ismagilova A.S., Spivak S.I. Inverse problems of chemical kinetics, Saarbrucken: Lap Lambert Acad. Publ. 2013. 118 p.

20. Denisov A.M., Dmitriev V.I. Inverse and incorrectly posed problems. Vestn. Moscovsk. Univ. Ser. 15. Vychislit. Matem. Kibernetika. 2005. N 5. P. 23-30 (in Russian).

21. Tarantola A. Inverse Problem Theory and Methods for Model Parameter Estimation. SIAM: Philadelphia. 2005. 344 p. DOI: 10.1137/1.9780898717921.

22. Kaipio J., Somersalo E. Statistical and computational inverse problems. New York: Springer. 2010. 339 p.

23. Chavent Guy. Nonlinear Least Squares for Inverse problems. New York: Springer. 2010. 292 p. DOI: 10.1007/97890-481-2785-6.

24. Aster R.C., Borchers B., Thurber C.H. Parameter estimation and inverse problems. New York: Elsevier. 2013. 360 p. DOI: 10.1016/B978-0-12-385048-5.00010-0.

26. Ягола AX., Янфей В., Степанова И.Э., Титаренко

В.Н Обратные задачи и методы их решения. М.: Бином. Лаборатория знаний. 2014. 216 с.

27. Федотов В.Х., Кольцов Н.И Разработка метода решения обратной задачи химической кинетики для каталитических реакций с участием основных веществ в каждой стадии.Хим. физика. 2016. Т. 35. № 10. С. 9-15. DOI: 10.7868/S0207401X1610006X.

28. Кольцов Н.И. Исследование адсорбции CO2 на хромок-сидном катализаторе по нелинейным временам релаксации. Изв. вузов. Химия и хим. технология. 2018. Т. 61. Вып. 2. С. 46-52.

29. Федотов В.Х., Кольцов Н.И. Исследование адсорбции CO2 на хромоксидном катализаторе по нестационарным концентрациям. Изв. вузов. Химия и хим. технология. 2018. Т. 61. Вып. 7. С. 37-43. DOI: 10.6060/ ivkkt.20186107.5714.

30. Bykov V.I., Tsybenova S.B., Yablonsky G.S. Chemical complexity via simple models. Berlin. New York: Germany. De Gruyter. 2018. 364 p. DOI: 10.1515/9783110464948.

31. Бабенко К.И. Основы численного анализа. М.-Ижевск: НИЦ «Регулярная и хаотическая динамика». 2002. 848 с.

25. Leonov A.S. Solution of ill-posed inverse problems: an essay on the theory, practical algorithms, and demonstrations in MATLAB. M.: Librokom. 2015. 336 p. (in Russian).

26. Yagola A.G., Yunfey V., Stepanova I.E., Titarenko V.N. Inverse problems and methods of their solution. M.: Binom. Laboratoriya znaniy. 2014. 216 p. (in Russian).

27. Fedotov V.Kh., Kol'tsov N.I. Method of Solving the Inverse Problem of Chemical Kinetics for Catalytic Reactions in Which Each Step Involves Main Reactants. Rus. J. Phys. Chem. B. 2016.V. 10. N 5. P. 753-759. DOI: 10.1134/S1990793116050195.

28. Kol'tsov N.I. Investigation of CO2 adsorption on a chromoxide catalyst by nonlinear relaxation times. ChemChemTech [Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol.]. 2018. V. 61. N 2. P. 46-52 (in Russian).

29. Fedotov V.Kh., Kol'tsov N.L Investigation of CO2 adsorption on a chromoxide catalyst by nonstationary concentrations. ChemChemTech [Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol.]. 2018. V. 61. N 7. P. 37-43 (in Russian). DOI: 10.6060/ivkkt.20186107.5714.

30. Bykov V.I., Tsybenova S.B., Yablonsky G.S. Chemical complexity via simple models. Berlin. New York: Germany. De Gruyter. 2018. 364 p. DOI: 10.1515/9783110464948.

31. Babenko K.I. Fundamentals of numerical analysis. M.Izhevsk: NITs «Regularnaya i khaoticheskaya dinamika». 2002. 848 p. (in Russian).

Поступила в редакцию 26.05.2021 Принята к опубликованию 15.11.2021

Received 26.05.2021 Accepted 15.11.2021

i Надоели баннеры? Вы всегда можете отключить рекламу.