CC BY
20
6. Lisetskii F. N., Marinina O. A., Gadzhiev R. S., Vorobyeva E. Ya. Rationale for indicators of arable farming duration (based on research findings in the county of the antique polis of Kerkinitis) // The Social Sciences. - 2016. - V. 11. - Is. 13. -P. 3361-3365.
7. Lisetskii F. N., Smekalova T. N., Marinina O. A. Biogeochemical features of fallow lands in the steppe zone // Contemporary Problems of Ecology. - 2016. - V. 9. - No. 3. - P. 366-375. http://doi.org/10.1134/S1995425516030094
8. Lisetskii F., Stolba V. F., Marinina O. Indicators of agricultural soil genesis under varying conditions of land use, Steppe Crimea // Geoderma. - 2015. - V. 239-240. - P. 304-316. http://dx.doi.org/ 10.1016/j.geoderma.2014.11.006
9. Mertz W. The essential trace elements // Science. - 1981. - V. 213. - P. 1332-1338.
10. Motavalli P. P., Jintaridth B., Lehmann J., Goyne K. W., Gilles J. Assessing and managing soil quality for sustainable agricultural system / Presented at the SANREM CRSP Annual Meeting, Los Banos, The Philippines, 26-29 May 2008 // Sustainable Agriculture and Natural Resource Management (SANREM) Knowledgebase. 2016. http://hdl.handle.net/10919/68519
11. Pardo T., Clemente R., Epelde L., Garbisu C., Bernal M. P. Evaluation of the phytostabilisation efficiency in a trace elements contaminated soil using soil health indicators // Journal of hazardous materials. - 2014. - V. 268. - P. 68-76.
12. Speidel D. H., Agnew A. F. The natural geochemistry of our environment. Boulder, Colo. : Westview Press, 1982. -214 p.
13. Teng Y., Wu J., Lu S., Wang Y., Jiao X., Song L. Soil and soil environmental quality monitoring in China: a review // Environment international. - 2014. - V. 69. - P. 177-199.
ГЕОЛОГО-МИНЕРАЛОГИЧЕСКИЕ НАУКИ / GEOLOGY AND MINERALOGY
DOI: https://doi.org/10.23670/IRJ.2017.57.143 Пушкин А.А.1, Римкевич В.С.2
1ORCID: 0000-0002-4579-0507, Кандидат физико-математических наук, 2 ORCID: 0000-0002-4579-0507, Кандидат геолого-минералогических наук, Институт геологии и природопользования, Дальневосточное отделение Российской академии наук УСТАНОВЛЕНИЕ ЗОН ГЕТЕРОФАЗНЫХ РЕАКЦИЙ
Аннотация
Данная работа посвящена вопросу установления зон гетерофазных реакций. Определение зон реакций для каждой температуры производится нами по минимуму погрешностей аппроксимаций экспериментальных данных различными функциями. В качестве пробных функций выбраны четыре вида уравнений: степенная функция, экспоненциальный закон и уравнения Авраами и Праута-Томпкинса. В случае равенства погрешностей аппроксимаций для двух функций предпочтение отдается функции, погрешность константы скорости для которой меньше. В качестве примера рассматривается реакция фторирования анортозитов гидродифторидом аммония. Результат определения зон реакций для данного примера согласуется с результатом определения по традиционному критерию (по величине энергии активации).
Ключевые слова: зона реакции, погрешность аппроксимации, погрешность константы скорости, гидродифторид аммония, энергия активации.
Pushkin A.A.1, Rimkevich V.S.2
1ORCID: 0000-0002-4579-0507, PhD in Physics and Mathematics,
2ORCID: 0000-0002-4579-0507, PhD in Geology and Mineralogy, Institute of geology and nature management, Far Eastern Branch of RAS THE DETERMINATION OF HETEROPHASE REACTIONS ZONES
Abstract
This work is devoted to the question of heterophase reactions zones identification. The determination of the reactions zones is executed by us by minimum of errors of experimental data approximations by various functions. The four functions are selected as probes: the power function, exponential law and Avraami and Prout-Tompkins equations. Preference is given to the function with less rate constants errors in the case of equality of approximation errors. The reaction of fluorination of anorthosites by ammonia hydrodifluoride is given as an example. The result of reactions zones determination is agreed with result received by the traditional criterion (by value of activation energy).
Keywords: reaction zone, approximation error, rate constant error, ammonia hydrofluoride, activation energy.
В Институте геологии и природопользования ДВО РАН разрабатывается фторидная технология переработки силикатного сырья. С целью проведения термодинамических (apriori) [1, с. 108] и кинетических (aposteriori) [2, с. 37] расчетов гетерофазных реакций в ходе технологического процесса создана расчетная программа на языке программирования Visual Basic.
Известно, что скорость гомофазных реакций определяется столкновениями между частицами. Результаты столкновений зависят от энергий частиц, точнее, от разницы между энергией активированного комплекса и средней энергией частиц. От столкновений скорость реакции зависит и в случае гетерофазных реакций. Если скорость реакции не только зависит, но и определяется столкновениями, то говорят, что такая реакция протекает в кинетической зоне. В случае, когда ход реакции определяется диффузией, говорят, что реакция протекает в диффузионной зоне.
Таким образом, гетерофазные реакции протекают в две стадии: диффузия частиц через межфазную границу и химическое взаимодействие между частицами. Каждая из стадий оказывает сопротивление протеканию реакции. Сопротивлением реакции называется величина, обратная её скорости. Поэтому эффективное сопротивление гетерофазной реакции между веществами А и В складывается из сопротивления кинетической и диффузионной стадий и, следовательно, справедливо соотношение
1 _ 1 1
кфф = * ТР' (1)
где к - константа скорости химической реакции, / - отношение площади контакта к объему зерна, ( -коэффициент массопереноса [3, с. 31].
Если первое слагаемое в выражении (1) много больше второго, то сопротивление кинетической стадии много больше сопротивления стадии диффузии, следовательно, скорость много меньше, а продолжительность - много больше. Тогда продолжительность этой стадии определяет продолжительность двухстадийной реакции и, следовательно, реакция протекает в кинетической зоне.
Если первое слагаемое значительно меньше второго, то сопротивление кинетической стадии много меньше, скорость много больше и продолжительность много меньше продолжительности стадии диффузии. Следовательно, продолжительность гетерофазной реакции в этом случае определяется продолжительностью стадии диффузии, т.е., реакция протекает в диффузионной зоне.
Если оба слагаемых сопоставимы по величине, то реакция протекает в переходной зоне (обе стадии оказывают одинаковые сопротивления ходу реакции и одинаковое влияние на ход реакции). В переходной зоне справедливо отношение равенства скоростей диффузии и химической реакции. Рассмотрим диффузию вещества А через контактную площадку S между зернами веществ А и В (см. рис. 1). Скорость диффузии в переходном слое из закона Фика (произведение коэффициента диффузии на градиент концентрации и на площадь контактной площадки) приравняем к скорости химической реакции
0^5* к С А СВУ' (2)
где О - коэффициент диффузии, 5 - площадь контакта зерен, V - объем зерна В, СА и СВ - концентрации фаз А и
В.
Полагая, что концентрация вещества А уменьшается до 0 в пределах толщины переходного слоя 8, представим градиент в виде простой дроби
оСА5ък С а С В-V' (3)
где 8 — толщина зоны контакта.
Сократив на концентрацию вещества А, разделив обе части на объем V и сгруппировав множители, производя замену отношения коэффициента диффузии Л к толщине зоны контакта 8 коэффициентом массоотдачи (, а отношение контактной площадки S к объему зоны контакта V - удельной поверхностью контакта /, получим в переходной области выражение
(Г* к Св. (4)
В кинетической зоне левая часть выражения (4) значительно больше, чем правая, а первое слагаемое в правой части (1) значительно больше, чем второе. Процессы, протекающие в кинетической стадии, подчиняются степенному закону: скорость реакции связана с концентрацией формулой IV = к С1, где ж - скорость реакции, С - концентрация вещества, т - порядок реакции, к - константа скорости.
В этом случае, вследствие высокой скорости диффузии частицы имеют возможность контакта между собой, постоянно сталкиваясь.
В диффузионной зоне первое слагаемое в формуле (1) много меньше, чем второе, а, следовательно, и правая часть формулы (4) много больше, чем левая.
Процессы диффузии протекают медленно и занимают большую часть времени реакции. Диффузионные процессы описываются одним из следующих топохимических уравнений, экспоненциальным ( а = 1 — ехр (—кэфф), Праута-
Томпкинса ( I) = — к,фф 1) или Авраами ( а = 1 — ехр(—кэфф ■ Р™)). Во всех этих уравнениях а — степень превращения вещества.
Следовательно, эффективная константа скорости гетерофазной реакции в диффузионной зоне зависит от диффузионных характеристик. Если в формуле (1) пренебречь первым слагаемым, то выражение для эффективной константы скорости с учетом формулы (3) приобретет вид
кэ ф ф = (/. (6)
Таким образом, гетерофазная реакции может определяться как столкновениями молекул, так и диффузией. Традиционным методом определения зоны реакции считается количественное сравнение энергий активации. Считается, что если энергия активации имеет величину меньше или порядка 20кДж/моль, то реакция протекает в диффузионной зоне. Если реакция имеет энергию активации порядка (50 - 200кДж/моль), то она протекает в кинетической зоне [4, с. 346].
В данной работе для решения вопроса о механизме реакции мы вычисляем погрешность аппроксимации экспериментальных данных теоретическими кривыми. Используются четыре типа кривых: степенная, Авраами, экспоненциальная и Праута-Томпкинса. Параметры кривых (константы скорости, порядки реакции) вычисляются путем обработки экспериментальных данных. Численно погрешность аппроксимации мы оцениваем как корень квадратный из суммы квадратов разности между экспериментальными значениями и теоретическими
значениями ( ^ ( )) для каждой из кривых, деленной на количество точек п:
Ргес(/(0) =
I п •
(7)
"VI
Расчеты проводились при помощи расчетной программы. Отбор функций, которые дают лучшую аппроксимацию, осуществляется автоматически по двум параметрам: первый параметр - погрешность аппроксимации, второй -погрешность константы скорости. В случае равенства погрешностей аппроксимации, программа предпочтение отдает той функции, погрешность константы скорости у которой меньше.
Выбранная функция определяет механизм реакции для данной температуры, а по механизму реакции -определяем зону реакции: при степенном законе - кинетическая, а при других трех уравнениях - диффузионная.
Данные о константах скоростей, порядках реакций, погрешностях аппроксимаций и констант скоростей для всех четырех кривых помещаются на вкладку Кинетика. Результаты отбора (по одной функции для каждой температуры) выносятся на вкладку Зона реакции (рис. 2).
На рис. 2 показано, что реакция спекания анортозитов с гидродифторидом аммония протекает при всех температурах по экспоненциальному закону, а, следовательно, зона реакции - диффузионная. Энергия активации также одна для всех температур 19,8кДж/моль. Константы скоростей монотонно возрастают от 0,003955 мин-1 при 1000С до 0,018159 мин-1 при 2000С.
По результатам данного эксперимента можно отметить, что, согласно традиционному критерию, при энергии активации 19,8 кДж/моль зона реакции также диффузионная.
В заключение отметим, что экспоненциальный механизм реакции означает медленное протекание диффузии. Медленное настолько, что зародыш, образовавшийся на поверхности зерна, успевает полностью превратиться в продукт реакции, прежде чем произойдет образование нового зародыша.
Математически экспоненциальный закон напоминает пуассоновское распределение, для которого вероятность появления 0 событий к моменту времени £ равна ехр (—кэфф () , где произведение константы скорости кэфф на время £ представляет собой параметр Пуассона [5, с. 175].
В нашем случае полезным исходом является появление частицы продукта. Тогда вероятность отсутствия полезного исхода равна ехр (—кэфф () . Вероятность же хотя бы одного полезного исхода равна 1 — ехр (—кэфф () , что и приводит нас к экспоненциальному закону протекания реакции.
При законе Авраами разность скоростей диффузии и скорости химической реакции не так велика. Уравнение Авраами представляет собой комбинацию экспоненциальной функции, характерной для диффузии, и степенной функции, характерной для химического взаимодействия частиц. Об этом говорит порядок реакции т, присутствующий в формуле в качестве характеристики химической реакции. Заметим, что порядок реакции т присутствует также в степенном законе, но отсутствует в экспоненциальном законе и в уравнении Праута-Томпкинса. Так что, возможно, уравнение Авраами характеризует гетерофазную реакцию не в диффузионной зоне, а в переходной.
Далее обратимся к уравнению Праута-Томпкинса. Оно представляет собой сумму геометрической прогрессии с первым членом ехр (—кэфф £) и знаменателем (—ехр (—кэфф £) ). Представим сумму в виде ряда
1+ехр ( ^фф^) = еХР (—кФФ£) ( 1 — еХР (—кФФ£) + еХР ( — 2 кФ Ф£) — еХР ( — 3 кФ Ф£) + ' ' ' ) ■ (8)
Первый, экспоненциальный, множитель этого выражения означает вероятность отсутствия благоприятного исхода в зародыше первого порядка, а первые два члена ряда описывают вероятность благоприятного исхода в зародыше второго порядка. Далее, если сгруппировать второй и третий члены ряда и вынести ( —ехр (— кэфф£) ) за скобки, то получится вероятность благоприятного исхода в зародыше третьего порядка и т.д. Поэтому процесс, описываемый уравнением Праута-Томпкинса, описывает появление частиц продукта в зародышах второго и более высоких порядков при отсутствии частицы продукта в первом порядке.
Таким образом, в данной работе мы проанализировали процессы, описываемые каждым из используемых нами четырех уравнений. Анализ показывает, что протекание процесса в соответствии со степенным законом происходит в кинетической зоне, а в соответствии с уравнениями экспоненциальным или Праута-Томпкинса - в диффузионной зоне. В случае если процесс описывается уравнением Авраами, то он, по -видимому, при данной температуре протекает в переходной зоне. Увеличение температуры сместит равновесие в ту или иную сторону.
Рис. 1 - Контакт двух зерен веществ А и В: 5- толщина переходного слоя в зоне контакта,
8 - площадь контакта
Рис. 2 - Вкладка Зона реакции программы Кинетика с результатами расчета реакции спекания анортозитов Котозерского массива Карелии с гидродифторидом аммония
Список литературы / References
1. Пушкин А.А., Римкевич В.С. Автоматизация расчетов изменений термодинамических потенциалов в ходе химических реакций// Международный научно-исследовательский журнал. 2015. №6 (37). Часть 1. С. 108-119.
2. Пушкин А.А., Римкевич В.С. Автоматизация расчетов кинетических параметров химических реакций// Инновационная наука. 2016. №2. В 5 частях. Часть 5. С.37-42.
3. В.В.Бочкарев. Теория химико-технологических процессов органического синтеза. Гетерофазные и гетерогенно-каталитические реакции. Учебное пособие. - Томск: Издательство ТПУ, 2005. - 118с.
4. А.Г.Стромберг, Д.П.Семченко. Физическая химия. - М.: Высшая школа, 1999. -528с.
5. В.Феллер. Введение в теорию вероятностей и её приложения. В 2-х томах. Т.1: Пер. с англ. - М.: Мир, 1984. -528с.
Список литературы на английском языке / References in English
1. Pushkin A., Rimkevich V.S. Avtomatizatsiya raschetov izmenenii termodinamicheskih potentsialov v hode himicheskih reaktsii [Automatization of calculations of thermodynamic potentials change during chemical reactions]//Mezhdunarodny nauchno-issledovatelskii Zhurnal [International research-scientifical journal]. 2015. №6 (37). Part 1. P. 108-119.
2. Pushkin A., Rimkevich V.S. Avtomatizatsiya raschetov kineticheskih parametrov himicheskih reaktsii [Automatization of calculations of kinetic parameters of chemical reactions]//Innovatsionnaya nauka [Innovation science]. 2016. №2. In 5 parts. Part 5. P. 37-42.
3. V.V.Bochkarev. Teoriya himiko-tehnologicheskogo sinteza. Heterofaznie I heterogenno-kataliticheskie reaktsii. Uchebnoe posobie [The theory of chemical technological synthesis. Heterophase and heterogeneous catalytic reactions. Tutorial]. - Tomsk: Izdatelstvo TPU [Publisher of TPU], 2005. 118p.
4. A.G.Stromberg, D.P.Semchenko. Fizicheskaya himiya [Physical chemistry]. - M.: Vysshaya shkola [High school], 1999. - 528p.
5. V.Feller. Vvedenie v teoriyu veroyatnostei i eio prilozheniya V 2-uh chastiyah [Introduction in theory of probability and her applications. In 2 parts]. V. 1: Perevod s angliiskogo [Translation from English]. - M.: Mir [Peace], 1984. - 184p.
