ВЛИЯНИЕ ВИДА КИНЕТИЧЕСКИХ ФУНКЦИЙ НА ЗАКОНОМЕРНОСТИ ПРЕВРАЩЕНИЯ В ТОНКОМ СЛОЕ, РАСПОЛОЖЕННОМ НА ПОДЛОЖКЕ, В УСЛОВИЯХ УПРАВЛЕНИЯ ПОДВИЖНЫМ ИСТОЧНИКОМ ТЕПЛА
- Авторы: Караулова А.Ю1, Князева А.Г2
- Учреждения:
- Национальный исследовательский Томский государственный университет, Томск, Российская Федерация
- Институт физики прочности и материаловедения СО РАН, Томск, Российская Федерация
- Выпуск: № 6 (2024)
- Страницы: 25-37
- Раздел: Статьи
- URL: https://ered.pstu.ru/index.php/mechanics/article/view/4388
- DOI: https://doi.org/10.15593/perm.mech/2024.6.03
- Цитировать
Аннотация
Современные комбинированные лазерные технологии предполагают сочетание возможностей поверхностной обработки с синтезом материала. В этом случае предварительные оценки температуры только на основе теплофизических моделей представляются необоснованными. Требуется учет химических реакций, которые могут протекать как в твердой, так и в жидкой фазах. Для этого в работе сделан краткий обзор особенностей реакций с участием твердых веществ и видов кинетических функций, которые применяются для обработки эксперимента в химии и химической кинетике. Разнообразие возможных вариантов кинетических функций предоставляет широкие возможности для анализа экспериментов, что способствует более глубокому пониманию происходящих процессов. Далее на основе одномерной редуцированной модели с суммарной химической реакцией и с учетом накопления жидкой фазы показано, что вид кинетической функции, которая отражает механизм реакции на микроуровне, влияет как на максимальную температуру, так и на степень завершенности реакции. Редуцированные модели позволяют нам более эффективно исследовать каждый изучаемый фактор, что впоследствии может быть полезно для оптимизации технологических процессов. Для системы Al+TiO2 представлены зависимости характеристик квазистационарной стадии от параметров, отражающих роль условий синтеза. В этой системе превращения можно описать двумя суммарными реакциями: первая дает упрочняющую фазу, вторая приводит к изменению состава матрицы. Зависимость максимальной температуры от плотности мощности источника тепла оказываются линейными. После обработки покрытие содержит как продукты реакций, так и исходные вещества. Конкретный состав определяется условиями обработки.
Полный текст
При получении новых материалов в современных технологиях важную роль играют физико-химические превращения, следствием которых и являются новые состав, структура и свойства [1, 2]. Особое внимание в последние десятилетия уделяется процессам, которые происходят в твердом состоянии. Это объясняется тем, что в некоторых случаях именно твердофазные реакции являются наиболее эффективным способом для получения новых веществ и материалов. В некоторых случаях только процессы в твердом состоянии могут преобразовать вещество в форму, которая отвечает заданным требованиям, то есть позволяют создать материал с необходимыми свойствами, к примеру, изготовить кристаллы с определенной концентрацией нужных дефектов [3]. К современным и развивающимся классическим технологиям, в которых важную роль играют химические реакции, можно отнести спекание в твердой и жидкой фазах [4], управляемый синтез покрытий [5], синтез композитов [6], создание материалов и изделий в 3D-технологиях [7, 8] и др. В твердых телах механизмы и динамика физико-химических процессов отличаются от тех, которые наблюдаются в жидкостях или газах [2, 9], в частности, из-за значительно меньшей скорости диффузии, которая в твердых материалах на порядки ниже [10]. Это замедление приводит к тому, что концентрация веществ в системе не успевает выравниваться, что вызывает пространственную локализацию протекающих процессов [2, 11]. Локализация процессов влияет на кинетику реакции, где важную роль играют как удельная скорость процесса (или коэффициент диффузии), так и геометрия реакционной области. Эти характеристики, зависящие от геометрии, называются топохимическими [12, 13]. Кроме того, поскольку реакции происходят в определенных пространственных зонах, их скорость может зависеть как от процессов на границе раздела фаз (реакционный контроль), так и от скорости переноса реагентов к этой границе или удаления продуктов реакции (диффузионный контроль) [9, 11]. Любая твердофазная реакция совершается в виде существенно большего набора промежуточных состояний (которые энергетически мало отличаются друг от друга, по сравнению с реакциями в газовой и жидкой фазах) [14]. В простых системах, к которым применимы определенные модельные предположения, можно экспериментально выявить эти режимы по характеру зависимости степени превращения от времени. Первые попытки разработать достоверные модели кинетики твердого тела относятся к концу 1920-х годов (Tamman G., Yander G., Roginsky S и др.). Дальнейшие исследования привели к разработке более полных кинетических моделей [15], учитывающих различные механизмы реакции на микроуровне [16-19]. В кинетике химических реакций в твердых телах потеря массы твердого вещества-реагента выражается в так называемой степени превращения (Y), которая изменяется в диапазоне от 0 до 1 и определяется как [18] , (1) где m0 – начальный вес образца, mt – вес образца в момент времени t, а m∞ – конечный вес образца. Однако если вес образца в ходе реакции не изменяется, степень превращения можно ввести иным способом. Положим, и – текущие массы реагентов и продуктов реакции, так что и . Величину тоже можем назвать степенью превращения или долей продукта реакции. Заметим, что подобное понятие можно ввести и в рамках термодинамики необратимых процессов [20]. Для простой реакции вида (2) ее скорость может быть выражена через изменение числа частиц любого сорта –реагентов или продуктов : , где – координата реакции, степень полноты реакции или степень превращения; и – стехиометрические коэффициенты. Скорость твердофазной реакции в химии твердого тела представляет собой изменение во времени степени превращения вещества и пропорциональна произведению константы скорости k(T) на дифференциальную скорость реакции φ(Y): . (3) Константа скорости обычно подчиняется закону Аррениуса , (4) где k0 – предэкспонент, Е – энергия активации химической реакции, φ – кинетическая функция, R – газовая постоянная. Т.е., в неизотермических условиях это уже не константа, а функция температуры. Дифференциальная скорость реакции – кинетическая функция – φ(Y) является функцией степени превращения, и ее форма зависит от основного механизма твердофазной реакции на микроуровне. Полноправными стадиями реакций с участием твердых веществ являются зарождение и рост малой доли продукта в местах несовершенств (нуклеация); диффузия, механизмы которой могут быть разнообразными (объемная и поверхностная, по границам зерен и фаз, по механизмам внедрения и замещения) и зависят от наличия дефектов; геометрические особенности, приводящие к преимущественному росту продукта в разных направлениях и т.д. [12, 13]. Интегрируя уравнение (2) при T=const по времени, можно получить так называемую интегральную форму скорости реакции g(Y): . (5) В таблице 1 приведены наиболее часто используемые модели, учитывающие различные механизмы реакции. Первый столбец таблицы включает названия моделей. Таблица 1/ Table 1 Набор различных моделей твердофазных реакций и их скорости, представленные в дифференциальной и интегральной формах. Set of different models of solid-phase reactions and their rates represented in differential and integral forms. Модели Дифференциальная форма Интегральная форма Модели нуклеации Степенной закон (P2) 2Y1/2 Y1/2 Степенной закон (P3) 3Y2/3 Y1/3 Степенной закон (P4) 4Y3/4 Y1/4 Аварами-Ерофеева (A2) 2(1 - Y)[- ln(1 - Y)]1/2 [- ln(1 - Y)]1/2 Аварами-Ерофеева (A3) 3(1 - Y)[- ln(1 - Y)]2/3 [- ln(1 - Y)]1/3 Аварами-Ерофеева (A4) 4(1 - Y)[- ln(1 - Y)]3/4 [- ln(1 - Y)]1/4 Проут-Томпкинса (B1) Y(1 - Y) ln[Y/(1 - Y)]+cY Модели геометрического сужения Сокращающаяся область(R2) 2(1 - Y)1/2 1-(1 - Y)1/2 Сокращающийся объем(R3) 3(1 - Y)2/3 1-(1 - Y)1/3 Модели диффузии 1-D диффузия (D1) 1/(2Y) Y2 2-D диффузия (D2) -[1/ ln(1 - Y)] ((1 - Y)ln(1 - Y))+Y 3-D диффузия-Яндера (D3) [3(1 - Y)2/3]/[2(1-(1 - Y)1/3)] (1-(1 - Y)1/3)2 Гистлинга-Браунштейна(D4) 3/[2((1 - Y)-1/3-1)] 1-(2/3)Y-(1 - Y)2/3 Модели порядка реакции Нулевого порядка(F0/R1) 1 Y Первого порядка (F1) (1 - Y) -ln(1 - Y) Второго порядка (F2) (1 - Y)2 [1/(1 - Y)]-1 Третьего порядка (F3) (1 - Y)3 (1/2)[(1 - Y)-2-1] Все используемые здесь модели делятся на четыре семейства, a именно: модели зарождения, геометрического сжатия, диффузии и модели порядка реакции. Таким образом, в таблице мы можем найти каждое семейство (выделено курсивом) и под ними соответствующие модели. Названия моделей, приведенные в таблице, представлены в виде буквы, присоединенной к целому числу (например A2). Буква отражает тип модели, a целое число связано с механизмом реакции, рассматриваемым в каждом конкретном случае. B частности, символы P, A, R, D, F означают закон мощности, закон Аврами-Ерофеева, изменение объема смеси, диффузию и реакцию-порядок, соответственно. Подробные сведения обо всех используемых здесь моделях, механизмах реакций и выводах каждой модели, представленных в таблице, можно найти в литературе [12, 13, 21-25]. Кроме того, на втором и третьем столбце таблицы приведены формы так называемых дифференциальных и интегральных скоростей реакций соответственно для различных моделей. В [26] предложено обобщенное уравнение, которое в наших переменных имеет вид , где q, β и γ – характеристические параметры для каждого конкретного механизма реакции. В макрокинетике [27, 28] при описании синтеза новых материалов в режиме послойного горения или теплового взрыва, при синтезе в растворах [29] и в условиях экструзии [30] используется подобный подход к описанию кинетики реакций. Однако имеются и более «продвинутые» модели. Например, для учета торможения реакции слоем твердофазного продукта на основе решения специальных задач используются кинетические законы, включающие параметры торможения [31, 32]. Так, если продукт реакции и реагенты (или хотя бы один из реагентов) не плавятся, химическое взаимодействие осуществляется путем диффузии реагентов через слой образовавшегося продукта, который разделяет реагенты и препятствует дальнейшему протеканию реакции. Полного смешения реагентов не происходит, и реакция носит существенно гетерогенный характер. Вид кинетической функции в таких системах определяется механизмом переноса реагентов и структурой нарастающего слоя продуктов, что отражает кинетическая функция вида , (6) где – параметры торможения. Наиболее распространены следующие законы: – линейный, – параболический, – кубический, – экспоненциальный. Если механизм реакции нельзя свести к единственной суммарной реакции (2), то можно использовать и классические представления химической кинетики о параллельных, последовательных, конкурирующих реакциях и др. [33]. Например, реакцию, протекающую в две последовательные стадии можно представить так: , где B – промежуточный продукт первой реакции и реагент для второй реакции. В этом случае понятия концентрации и степени превращения не тождественны. Систему кинетических уравнений можно записать как относительно концентраций веществ, определенных аналогично предыдущему , так и относительно координат первой и второй стадий. Для первого варианта имеем систему кинетических уравнений , , а . Заметим, что классический формально-кинетический подход к описанию кинетики реакций, учитывающий параллельные и последовательные стадии, может быть вполне модифицирован за счет учета торможения реакций слоем продукта [34, 35]. Тогда для последовательной реакции с неплавящимся конечным продуктом P . Обзор двухуровневых моделей, в которых явно учитываются особенности диффузионно-контролируемых реакций, представлен в [36]. К сожалению, для большинства твердофазных реакций такой подход пока остается весьма условным. Цель настоящей работы состоит в исследовании влияния вида кинетического закона на термокинетические закономерности синтеза покрытия на подложке.Об авторах
А. Ю Караулова
Национальный исследовательский Томский государственный университет, Томск, Российская Федерация
А. Г Князева
Институт физики прочности и материаловедения СО РАН, Томск, Российская Федерация
Список литературы
- Spivak L. V., Shchepina N.E. Fiziko-khimicheskie osnovy protsessov mikro- i nanotekhnologii [Physicochemical foundations of micro- and nanotechnology processes]: ucheb. posobie: v 2 ch. – Perm. gos. nats. issled. un-t. – Perm', 2018. – Pt. 1. – 202 p.
- Butjagin P. Ju. Himicheskaja fizika tverdogo tela [Chemical physics of solids] : ucheb. dlja studentov, obuchajushhihsja po napravleniju 511700 "Himija, fizika i mehanika materialov" / Moskva : Izd-vo Mosk. un-ta (MGU), 2006. – 269 p
- Gorshkov V.S. Fizicheskaja himija silikatov i drugih tugoplavkih soedinenij [Physical chemistry of silicates and other refractory compounds], [Ucheb. dlja vuzov po spec. "Him. tehnologija tugoplavkih nemetal. i silikat. materialov"] / V.S. Gorshkov, V.G. Savel'ev, N.F. Fedorov. — M. : Vyssh. shkola, 1988. — 399 p
- Anciferov V.N., Gilev V.G. O roli ob#emnyh i massovyh jeffektov reakcij v processah reakcionnogo spekanija [On the role of volume and mass effects of reactions in reaction sintering processes]// Izv. vuzov. Poroshk. metallurgija i funkc. pokrytija. 2015. № 4. pp. 9–20 DOI: dx.doi.org/10.17073/1997-308X-2015-4-9-2
- Qilong Wu, Weimin Long, Lei Zhang, Hongwei Zhao, A review on ceramic coatings prepared by laser cladding technology, Optics Laser Technology, Volume 176, 2024, 110993, doi: 10.1016/j.optlastec.2024.11099
- Zigui Luo, Qi Jin, Rui Wang, Zhenjia Yang, Zhaoyang Wu, Guangqiang Li, Xi’an Fan High-performance FeSiAl/(Al 2 O3 -Ni) soft magnetic composites prepared by in situ synthesis method // Advanced Powder Technology 35 (2024) 104624 doi: 10.1016/j.apt.2024.10462
- Shishkovskij I.V., Jadroitcev I.A., Smurov I.Ju. Sozdanie ob#emnyh izdelij iz nikelida titana metodom poslojnogo lazernogo plavlenija [Creation of volumetric products from titanium nickelide by the method of layer-by-layer laser melting]// Pis'ma v ZhTF, 2013, tom 39, vyp. 24. pp.15-2
- Fomin V.M. [i dr.] In-situ sintez v processe prjamogo lazernogo vyrashhivanija funkcional'nogo gradientnogo materiala na osnove polikristallicheskogo bora i splava Ti64 [In-situ synthesis in the process of direct laser growth of a functional gradient material based on polycrystalline boron and Ti64 alloy]// Fizicheskaja mezomehanika, 2022, T.25, № 4. pp.34-43. doi: 10.55652/1683-805X_2022_25_4_3
- Tretyakov Yu. D. Tverdofaznye reaktsii [Solid phase reactions]. M.: Khimiya Publ., 1978, 360 p
- Mehrer H. Diffusion in Solids. Springer Series in Solid-State Sciences, 2007, vol. 155. 654 p. doi: 10.1007/978-3-540-71488-
- Braun. M. Reakcii tverdyh tel [Reactions of Solids] / M. Braun, D. Dollimor, A. Galvej ; red. V. V. Boldyrev . – M. : Mir, 1983 . – 359 p
- Rozovskij A. Ja. Kinetika topohimicheskih reakcij [Kinetics of topochemical reactions]. - M.: Himija, 1976. – 219 p
- Del'mon B. Kinetika geterogennyh reakcij [Kinetics of heterogeneous reactions]: Per. s franc.1972. 556 p
- Tretyakov Yu.D. Solid phase reactions //. Soros Educational Journal, 1999 – N4 – pp. 35-3
- Wyazovkin S., Wight C.A. Isothermal and non-isothermal kinetics of thermally stimulated reactions of solid // International Reviews in Physical Chemistry, 1998, Vol . 17, No . 3, 407- 43
- Vyazovkin, S.; Wight, C.A. Isothermal and Nonisothermal Reaction Kinetics in Solid: In Search of Ways toward Consensus. J. Phys. Chem. A 1997, 101, 8279–8284. doi: 10.1021/jp971889
- Roduit, B. Computational aspects of kinetic analysis. Part E: The ICTAC Kinetics Project-numerical techniques and kinetics of solid state processes // Thermochim. Acta 2000, 355, 171–180. doi: 10.1016/S0040-6031(00)00447-
- Khawam, A.; Flanagan, D.R. Solid-state kinetic models: Basics and mathematical fundamentals. J. Phys. Chem. B 2006, 110, 17315–17328. doi: 10.1021/jp062746a;
- Khawam, A.; Flanagan, D.R. Basics and Applications of Solid-State Kinetics: A Pharmaceutical Perspective. J. Pharm. Sci. 2006, 95, 472–498. doi: 10.1002/jps.20559 5-
- Prigozhin I., Kondepudi D. Sovremennaja termodinamika. Ot teplovyh dvigatelej do dissipativnyh struktur [Modern Thermodynamics: From Heat Engines to Dissipative Structures]/ 2002. 464 p.
- Yannis Vasilopoulos, Eliška Skoˇrepová, Miroslav Šoóš. COMF: Comprehensive Model-Fitting Method for Simulating Isothermal and Single-Step Solid-State Reactions // Crystals. – 2020. – 10 – 139. DOI:0.3390/cryst1002013
- Capart, R.; Khezami, L.; Burnham, A.K. Assessment of various kinetic models for the pyrolysis of a microgranular cellulose // Thermochim. Acta 2004, 417, 79–89. doi: 10.1016/j.tca.2004.01.02
- Zhou, D.; Grant, D.J.W. Model dependence of the activation energy derived from nonisothermal kinetic data // J. Phys. Chem. A 2004, 108, 4239–4246. doi: 10.1021/jp037917
- Vyazovkin, S.; Chrissafis, K.; Di Lorenzo, M.L.; Koga, N.; Pijolat, M.; Roduit, B.; Sbirrazzuoli, N.; Suñol, J.J. ICTAC Kinetics Committee recommendations for collecting experimental thermal analysis data for kinetic computations // Thermochim. Acta 2014, 590, 1–23. doi: 10.1016/j.tca.2014.05.03
- Dickinson, C.F.; Heal, G.R. A review of the ICTAC Kinetics Project, 2000: Part 1. Isothermal results // Thermochim. Acta 2009, 494, 1–14. doi: 10.1016/j.tca.2009.05.003
- Urbanovici E., Segal E. General kinetic equation for solid state reactions // Journal of thermal analysis and Calorimetry, 1999, V.55.P. 919-92
- Frank-Kameneckij D.A. Osnovy makrokinetiki. Diffuzija i teploperedacha v himicheskoj kinetike [Fundamentals of Macrokinetics. Diffusion and heat transfer in chemical kinetics]: uchebnik-monografija, 4th ed. – Dolgoprudnyj: Intellekt, 2008. 408 p.
- Rogachev A.S., Mukasjan A.S. Gorenie dlja sinteza materialov: vvedenie v strukturnuju makrokinetiku [Combustion for the synthesis of materials: introduction to structural macrokinetics] 2013. 400 p
- Materialoobrazujushhie vysokojekzotermicheskie processy: metallotermija i gorenie sistem termitnogo tipa [Material-forming highly exothermic processes: metallothermy and combustion of thermite-type systems] (Alymov M.I., Juhvid V.I., Levinskij Ju.V. - red) / M.: 2021, 376 p
- A. M. Stolin and P. M. Bazhin SHS Extrusion: An Overview // nternational Journal of SelfPropagating HighTemperature Synthesis, 2014, Vol. 23, No. 2, pp. 65–73. doi: 10.3103/S106138621402011
- A.P. Aldushin, A.G. Merzhanov, B.I. Khaikin Some features of combustion of condensed systems with high melting-point reaction products // Dokl. Akad. Nauk SSSR, 204 (1972), pp. 1139-114
- A.P. Aldushin [i dr.] Propagation of the front of an exothermic reaction in condensed mixtures with the interaction of the components through a layer of high-melting product // Comb. Explos. Shock Waves, 8 (1972), pp. 202-21
- Semiohin, I.A., Strakhov, B.V. and Osipov, A.I. (1995) Kinetics of Chemical Reactions. M. Publishing house of Moscow State University. 351 p
- Knjazeva A.G., Sorokova S.N. Modelirovanie formirovanija fazovoj struktury pokrytija v processe jelektronno–luchevoj obrabotki s ispol'zovaniem sinteza v tverdoj faze // Teoreticheskie osnovy himicheskoj tehnologij [Simulation of coating phase structure formation in solid phase synthesis assisted by electron-beam treatment // Theoretical Foundations of Chemical Engineering] 2008. T.42. №4. pp. 457– 46
- Sorokova S.N., Knjazeva A.G. Matematicheskoe modelirovanie ob’emnyh izmenenij v processe spekanija poroshkov sistemy T-Al // Fizicheskaja mezomehanika [Mathematical modeling of volumetric changes in the process of sintering powders of the T-Al system // Physical mesomechanics], 2008, T. 11, № 6. pp.95-10
- Knjazeva A.G. Two-Level Models of Composite Synthesis: History and Potential // Combustion Explosion and Shock Waves, 2024, Т.60, № 1. pp.48-62. doi: 10.15372/FGV2023.933
- Knyazeva A.G., Anisimova M.A. Two-level model of composite synthesis on a substrate from a powder mixture of TiO2 and Al. 1. Model description and simple examples // Nanoscience and Technology: An International Journal, 2023, V.14, No3. P.73-102. doi: 10.1615/NanoSciTechnolIntJ.202304704
- Knjazeva, A.G. Termomehanicheskaja svjazannaja model' sinteza pokrytija na podlozhke [Thermomechanical coupled model of coating synthesis on a substrate] / A.G. Knjazeva // Vestnik Permskogo nacional'nogo issledovatel'skogo politehnicheskogo universiteta. Mehanika. – 2024. – № 1. – pp. 58–74. doi: 10.15593/perm.mech/2024.1.0
- Fan R., Liu B., Zhang J., Bi J., Y. Yin Kinetic evaluation of combustion synthesis 3TiO2 + 7Al → 3TiAl + 2Al2O3 using non-isothermal DSC method // Materials Chemistry and Physics Volume 91, Issue 1, 15 May 2005, Pages 140-145 doi: 10.1016/j.matchemphys.2004.11.00
- Zhu H. G., Min J., Ai Y. L., Wu Q. Reaction Mechanism of an Al-TiO2 System // Advanced Materials Research. – 2010. – Vols. 97–101. – P. 1624–1627. doi: 10.4028/www.scientific.net/amr.97-101.16