ON THE CHEMICAL AFFINITY TENSOR IN CHEMO-MECHANICS OF DEFORMABLE SOLIDS

Abstract


В работе выводится выражение тензора химического сродства для химических реакций в деформируемом теле. Рассматривается реакция между диффундирующим и деформируемым твердым компонентами. Продукт реакции также является деформируемым твердым материалом. Реакция сопровождается собственной деформацией превращения. Для открытой системы «деформируемые и диффундирующая компоненты», в которой протекает химическая реакция, локализованная на фронте реакции, для случая конечных деформаций записываются балансы массы, импульса, энергии и второй закон термодинамики в виде неравенства Клаузиуса – Дюгема. Вводится параметр, характеризующий степень участия в реакции диффундирующего компонента, расположенного по разные стороны фронта реакции. В результате выводится выражение для производства энтропии вследствие химической реакции на распространяющемся фронте реакции. Соответствующее выражение для диссипации энергии имеет вид произведения нормальной компоненты тензора, имеющего смысл тензора химического сродства, на скорость реакции на ориентированной площадке фронта. Выражение следует из фундаментальных законов без использования определяющих соотношений компонентов реакции. Тензор сродства равен линейной комбинации тензоров химического потенциала твердых компонентов реакции, равных тензорам энергии-импульса Эшелби, деленным на плотности материала, и химического потенциала диффундирующего компонента. Тензорность сродства объясняется протеканием реакции не в точке деформируемого тела, а на ориентированной площадке, проходящей через точку, подобно тому, как усилия в негидростатически напряженном теле, в отличие от давления в точке, рассматривают на ориентированных площадках и определяют тензором напряжений. С точки зрения механики конфигурационных сил тензор сродства определяет конфигурационную силу, управляющую распространением фронта реакции. Непротиворечивость результата подтверждается совпадением с классическим выражением скалярного сродства в случае гидростатически напряженного тела. Выведенное выражение тензора сродства дает возможность сформулировать кинетическое уравнение для реакции, в котором напряженно-деформированное состояние влияет на реакцию через сродство. Это уравнение является замыкающим для связанной задачи «химия–диффузия–механика». Важность результата обусловлена также возможностью его использования при моделировании объемных химических реакций, которые на мезоуровне протекают на распространяющихся фронтах реакции.

Full Text

Хемомеханика изучает взаимосвязи химических реакций и напряженно-деформируемого состояния. В данной работе рассматривается реакция между деформируемым твердым и диффундирующим компонентами. Продукт реакции также является деформируемым твердым материалом. Реакция сопровождается собственной деформацией превращения, которая порождает напряжения, влияющие в свою очередь на скорость реакции, и изменением механических параметров твердых компонентов, таких как модули упругости, и, вообще говоря, изменением реологических свойств материала, например, упругий материал может становиться вязкоупругим. К таким процессам можно отнести возникновение напряжений в оксидных пленках [1], при окислении кремния при создании и эксплуатации микроэлектронных устройств (см., например, [2–5]). Важными примерами взаимосвязи химической реакции, деформирования и разрушения являются окисление микронных элементов МЭМС [6–8], формирование интерметаллидов в бессвинцовых припоях [9–11], коррозии металлов под напряжением [12–16]. Особую актуальность эти исследования приобретают в связи с разработкой литий-ионных батарей (см., например, [17 – 20]). Одним из наиболее перспективных материалов для анода литий-ионных батарей является кремний (см., например, [21, 22]). Но реакция литизации кремния сопровождается деформацией превращения до 300% [23], которая порождает напряжения, ведущие к замедлению реакции и деградации деформационно-прочностных свойств материала анода в циклах заряда-разряда. При этом на первых циклах реакция локализована на распространяющемся фронте превращения, а на последующих циклах протекает в объеме [24]. Влияние напряжений на химическую реакцию может моделироваться через влияние на диффузию диффундирующего компонента реакции (см., например, [25 – 30]). Но многие реакции контролируются непосредственно скоростью реакции, если скорость диффузионного подвода достаточно велика (см., например, [31, 32]). В этом случае требуется теория, описывающая влияние напряжений именно на скорость реакции. В классической физической химии химические реакции описывают с помощью параметра, называемого степенью превращения (степенью завершенности реакции) [33], при этом скорость реакции – скорость изменения этого параметра – определяется скалярным химическим сродством, которое появляется как множитель в выражении производства энтропии вследствие химической реакции. В работах [34 – 36] был развит подход к описанию локализованных реакций на основе полученного выражения тензора химического сродства (см. также [37, 38]). Было показано, что диссипация энергии вследствие распространения фронта химической реакции равна произведению нормальной компоненты тензора сродства на скорость реакции на ориентированной площадке фронта аналогично произведению скалярного сродства на скорость реакции в точке при классических подходах. Тензорность сродства в случае реакции в деформируемом твердом теле объясняется рассмотрением протекания реакции не в точке тела, а на ориентированной площадке, проходящей через точку, подобно тому, как усилия в негидростатически напряженном теле, в отличие от давления в точке, рассматривают на ориентированных площадках и определяют тензором напряжений. Было показано, что тензор химического сродства определяется тензорами химического потенциала твердых компонентов, равных тензорам энергии-импульса Эшелби, деленным на плотность материала (тензорами химических потенциалов твердых компонентов реакции), и химическим потенциалом диффундирующего компонента, доставленного к фронту реакции диффузией через превращенный материал. Тогда скорость реакции на фронте и, следовательно, скорость фронта определяются нормальной компонентой тензора сродства. Напряжения при этом влияют на скорость реакции через сродство. На основании развитого подхода был решен ряд связанных задач «химия – диффузия – механика», продемонстрировавших, что механические напряжения могут замедлять и блокировать распространение фронта реакции (см., например [39 – 41]). При выводе формулы тензора сродства рассматривалась локализованная реакция между твердым и диффундирующим компонентами, которая происходит в результате доставки диффундирующего компонента к фронту реакции через превращенный материала. Диффузия в непревращенном материале не происходила. Это предположение согласуется с экспериментальными данными (см., напр., [24]), но оставило открытым вопросы о влиянии диффузии в непревращенном материале, а также роли тензора сродства в случае объемных химических реакций. Эти вопросы взаимосвязаны. Возможность диффузии впереди локализованного фронта превращения может привести к «размыванию» фронта реакции и переходу в режим объемной реакции. Но объемная химическая реакция может быть рассмотрена на двух масштабных уровнях, называемых далее макро- и мезоуровнями. Если с самого начала рассматривать объемную реакцию как изменение степени превращения как макроскопического (согласно определениям [42] – макрофеноменологического) параметра состояния, то тензорность химического сродства может оказаться скрытой в макрофеноменологическом определяющем уравнении в виде зависимости скорости превращения от скалярного сродства. В работе [43] для объемных реакций в деформируемом теле при ряде дополнительных предположений было получено выражение для скалярного химического сродства как множителя при скорости реакции в выражении для производства энтропии. Рассмотрение велось на макроуровне, когда реакция в точке тела была представлена изменением степени превращения – объемной доли превращенного материала, а напряжения и деформации в точке вычислялись по правилу смеси. Было показано, что это сродство формально может быть вычислено как нормальная компонента полученного ранее тензора сродства, если, в духе подходов метода фазового поля определить нормаль как нормированный градиент степени превращения [44, 45]. При этом остались скрытыми физико-химические процессы, протекающие в представительном объеме на следующем масштабном уровне. Представительный объем, соответствующий материальной точке макроуровня, на мезоуровне представляет собой смесь «чистых» исходного и полностью превращенного материала, через которую проходит диффундирующий компонент. Твердые компоненты разделены границами – локальными фронтами реакции. Скорость реакции на макроуровне определяется скоростью фронтов реакции на мезоуровне, движение которых меняет соотношение объемов, занятых превращенным и непревращенным материалами в представительном объеме, меняя степень превращения на макроуровне. При этом необходимо, чтобы на мезоуровне диффузия проходила как через превращенный, так и непревращенный материалы. Поэтому учет влияния диффузии в непревращенном материале является важным не только для моделирования реакций, локализованных на макроуровне, когда тело в целом разделено фронтом реакции на области превращенного и непревращенного материалов, но и для объемных реакций, если принимать во внимание физико-химические процессы, протекающие на масштабе мезоуровня, определение которого требует отдельного обсуждения (см., напр., монографию [42]), и ведущие также к изменению степени локализации реакции. В данной работе выводится выражение тензора сродства в случае локализованной реакции с учетом возможности диффузии через весь материал и участия в реакции диффундирующего компонента, находящегося по обе стороны фронта реакции. Рассмотрение ведется на макроуровне, когда изначально предполагается, что реакция локализована на фронте превращения. Это также открывает возможность для разработки модели объемных химических реакций, в которой скорость превращения на макроуровне определяется кинетикой распространения фронтов превращения в представительном объеме, а макроскопическим величинам соответствуют средние значения соответствующих величин на мезоуровне.

About the authors

A. B Freidin

Institute for Problems in Mechanical Engineering of the Russian Academy of Sciences

E. V Mikheev

Institute for Problems in Mechanical Engineering of the Russian Academy of Sciences

References

  1. Stringer J. Stress generation and relief in growing oxide films // Corrosion Science. – 1970. – pp. 513-543. https://doi.org/10.1016/S0010-938X(70)80036-
  2. Kao D., McVitie J., Nix W., Saraswat K. Two-dimensional thermal oxidation of silicon-I. Experiments // IEEE Trans. Electron Devices. – Vol. ED-34. – pp. 1008-1017. https://doi.org/10.1109/T-ED.1987.2303
  3. Kao D., McVitie J., Nix W., Saraswat K. Two-dimensional thermal oxidation of silicon-II. Modeling stress effect in wet oxides. // IEEE Trans. Electron Devices. – Vol. ED-35. – pp. 25-37. http://dx.doi.org/10.1109/16.241
  4. Mihalyi A., Jaccodine R.J., Delph T.J. Stress effects in the oxidation of planar silicon substrates // Applied Physics Letters. – Vol. 74. – pp. 1981-1983. https://doi.org/10.1063/1.12372
  5. Yen J.Y., Hwu J.G. Enhancement of silicon oxidation rate due to tensile mechanical stress // Applied Physics Letter. – 2000. – Vol. 76. – pp. 1834-1835. http://dx.doi.org/10.1063/1.12618
  6. Muhlstein C., Brown S., Ritchie R. High-cycle fatigue and durability of polycrystalline silicon thin films in ambient air // Sensors and Actuators. – 2001. – Vol. A94. – pp. 177-188. http://dx.doi.org/10.1557/PROC-657-EE5.
  7. Muhlstein C., Ritchie R. High-cycle fatigue of micron-scale polycrystalline silicon films: fracture mechanics analyses of the role of the silica/silicon interface // Int. J. of Fract. – 2003. – Vol. 119/120. – pp. 449-474. http://dx.doi.org/10.1023/A:102498803139
  8. Muhlstein C., Stach E., Ritchie R. A reaction-layer mechanism for the delayed failure of micron-scale polycrystalline silicon structural films subjected to high-cycle fatigue loading // Acta Materialia. – 2002. – Vol. 50. – pp. 3579-3595. https://doi.org/10.1016/S1359-6454(02)00158-
  9. Böhme T., Müller W.H., Weinberg K. Numerical modeling of diffusion induced phase transformations in mechanically stressed lead-free alloys // Computational Materials Science. – 2009. – Vol. 45. – pp. 837-844. https://doi.org/10.1016/j.commatsci.2008.08.01
  10. Wang Y., Wang Y., Ma L., Han L., Guo F. Effect of Sn grain c-axis on Cu atomic motion in Cu reinforced composite solder joints under electromigration // Journal of Electronic Materials. – 2020. – Vol. 49. – pp. 2159-2163. http://dx.doi.org/10.1007/s11664-019-07897-
  11. Morozov A., Freidin A.B., Kinkov V.A., Semencha A.V., Muller W.H., Hauck T. Experimental and theoretical studies of Cu-Sn intermetallic phase growth during high-temperature storage of eutectic SnAg interconnects // Journal of Electronic Materials. – 2020. – Vol. 49. – pp. 7194-7210. http://dx.doi.org/10.1007/s11664-020-08433-
  12. Gutman E., Bergman R., Levitsky S. Influence of internal uniform corrosion on stability loss of a thin-walled spherical shell subjected to external pressure // Corrosion Science. – 2016. – Vol. 111. – pp. 212-215. https://doi.org/10.1016/j.corsci.2016.04.01
  13. Pronina Y. An analytical solution for the mechanochemical growth of an elliptical hole in an elastic plane under a uniform remote load // European Journal of Mechanics — A/Solids. – 2017. – Vol. 61. – pp. 357-363. https://doi.org/10.1016/j.euromechsol.2016.10.00
  14. Pronina Y., Khryashchev S.M. Mechanochemical growth of an elliptical hole under normal pressure. // Materials Physics and Mechanics. – 2017, – Vol. 31. – pp. 52-55. https://doi.org/10.1016/j.euromechsol.2016.10.00
  15. Sedova O., Pronina Y. The thermoelasticity problem for pressure vessels with protective coatings, operating under conditions of mechanochemical corrosion // International Journal of Engineering Science. – 2022. – Vol. 170. https://doi.org/10.1016/j.ijengsci.2021.10358
  16. Evstafeva I., Pronina Y. On the mechanochemical dissolution of shells and its temperature dependence: Discussion of dfferent models // International Journal of Engineering Science. – 2023. – Vol. 190. https://doi.org/10.1016/j.ijengsci.2023.10388
  17. Sethuraman V., Van Winkle N., Abraham D., Bower A., Guduru P. Real-time stress measurements in lithium-ion battery negative-electrodes // Journal Power Sources. – 2012. – Vol. 206. – pp. 334-342. https://doi.org/10.1016/j.jpowsour.2012.01.03
  18. McDowell M.T., Lee S.W., Nix W.D., Cui Y. 25th anniversary article: Understanding the lithiation of silicon and other alloying anodes for lithium-ion batteries // Advanced Materials. – 2013b. – Vol. 25. – pp. 4966-4985. https://doi.org/10.1002/adma.20130179
  19. Tian J., Chen Z., Zhao Y. Review on modeling for chemo-mechanical behavior at interfaces of all-solid-state lithium-ion batteries and beyond // ACS Omega. – 2022. – Vol. 7. – pp. 6455-6462. http://dx.doi.org/10.1021/acsomega.1c0679
  20. McDowell M., Cortes F.J.Q., Thenuwara A.C., Lewis J.A. Toward high-capacity battery anode materials: Chemistry and mechanics intertwined // Chemistry of Materials. – 2020. – Vol. 28. – pp. 8755-8771. http://dx.doi.org/10.1021/acs.chemmater.0c0298
  21. Wu H., Cui Y. Designing nanostructured si anodes for high energy lithium-ion batteries // Nano Today. – 2012. – Vol. 7. – pp. 414-429. https://doi.org/10.1016/j.nantod.2012.08.00
  22. Liang B., Liu Y., Xu Y. Silicon-based materials as high-capacity anodes for next generation lithium-ion batteries // Journal of Power Sources. – 2014. – Vol. 267. – pp. 469-490. https://doi.org/10.1016/j.jpowsour.2014.05.09
  23. McDowell M., Ryu I., Lee S., Wang C., Nix W., Cui Y. Studying the kinetics of crystalline silicon nanoparticle lithiation with in situ transmission electron microscopy // Adv. Mater. – 2012. – Vol. 24. – pp. 6034-6041. https://doi.org/10.1002/adma.20120274
  24. McDowell M., Lee S.W., Harris J.T., Korgel B.A., Wang C.M., Nix W.D., Cui Y. In situ TEM of two-phase lithiation of amorphous silicon nanospheres // Nano Letters. – 2013a. – Vol. 13. – pp. 758-764. https://doi.org/10.1021/nl304450
  25. Knyazeva A. Cross effects in solid media with diffusion // Journal of Applied Mechanics and Technical Physics. –2003. – Vol. 44. – pp. 373-384. https://doi.org/10.1023/A:102348522403
  26. Knyazeva A. Application of irreversible thermodynamics to diffusion in solids with internal surfaces. // Journal of Non-Equilibrium Thermodynamics. – 2020. – Vol. 45. – pp. 401-418. http://dx.doi.org/10.1515/jnet-2020-002
  27. Knyazeva A. Reaction rates in a deformable solid in terms of the thermodynamics of irreversible processes // Russian Physics Journal. – 2022. – Vol. 64. – pp. 1907-1914. http://dx.doi.org/10.1007/s11182-022-02540-
  28. Levitas V., Attariani H. Anisotropic compositional expansion and chemical potential for amorphous lithiated silicon under stress tensor // Sci. Rep. – 2013. – Vol. 3. https://doi.org/10.1038/srep0161
  29. Brassart L., Suo Z. Reactive ow in solids // J. Mech. Phys. Solids. – 2013. – Vol. 61. – pp. 61-77. https://doi.org/10.1016/j.jmps.2012.09.00
  30. Zhang X., Zhong Z. A coupled theory for chemically active and deformable solids with mass diffusion and heat conduction // Journal of the Mechanics and Physics of Solids. – 2017, – Vol. 107. – pp. 49-75. https://doi.org/10.1016/j.jmps.2017.06.01
  31. Liu X., Zheng H., Zhong L., Huang S., Karki K., Zhang L., Liu Y., Kushima A., Liang W., Wang J., et al. Anisotropic swelling and fracture of silicon nanowires during lithiation // Nano Lett. – 2011. – Vol. 11. – pp. 3312-3318. http://dx.doi.org/10.1021/nl201684
  32. Jia Z., Li T. Stress-modulated driving force for lithiation reaction in hollow nano-anodes // Journal of Power Sources. – 2015. – Vol. 275. – pp. 866-876. https://doi.org/10.1016/j.jpowsour.2014.11.08
  33. Prigogine I., Defay R. Chemical thermodynamics. London: Longmans. – 1988
  34. Freidin A.B. Chemical affinity tensor and stress-assist chemical reactions front propagation in solids // ASME International Mechanical Engineering Congress and Exposition, Volume 9: Mechanics of Solids, Structures and Fluids, p. V009T10A102. – San Diego, California, USA. – 2013. http://dx.doi.org/10.1115/IMECE2013-6495
  35. Freidin A.B., Vilchevskaya E.N., Korolev I.K. Stress-assist chemical reactions front propagation in deformable solids // International Journal of Engineering Science. – 2014. – pp. 57-75. https://doi.org/10.1016/j.ijengsci.2014.03.00
  36. Фрейдин А.Б. О тензоре химического сродства при химических реакциях в деформируемых материалах // Механика твердого тела. – 2015. – Т. 3. – С. 35-68
  37. Freidin A.B., Vilchevskaya E.N. Chemical affnity tensor in coupled problems of mechanochemistry // In: Encyclopedia of Continuum Mechanics. Berlin: Springer Berlin Heidelberg. – 2020. – pp. 264-280
  38. Фрейдин А.Б. О конфигурационных силах в механике фазовых и химических превращений // Прикладная математика и механика. – 2022. – Т. 86, № 4. – С. 571-583
  39. Poluektov M., Freidin A.B., Figiel Ł. Modelling stress-affected chemical reactions in non-linear viscoelastic solids with application to lithiation reaction in spherical Si particles // International Journal of Engineering Science. – 2018. – Vol. 128. – pp. 44-62. https://doi.org/10.1016/j.ijengsci.2018.03.00
  40. Morozov A., Freidin A.B., Müller W.H. On stress-affected propagation and stability of chemical reaction fronts in solids // International Journal of Engineering Science. – 2023. https://doi.org/10.1016/j.ijengsci.2023.10387
  41. Freidin A.B., Rublev I.A., Korolev I.K. On the initial boundary value problem for the propagating chemical reaction front in an elastic solid // Acta Mechnica. – 2023. – Vol. 235. – pp. 4289-4310. http://dx.doi.org/10.1007/s00707-024-03945-
  42. Трусов П.В., Швейкин А.И. Многоуровневые модели моно- и поликристаллических материалов: теория алгоритмы, примеры применения. – М.: Издательство сибирского отделения российской академии наук, 2019. – 605 с
  43. Poluektov M., Freidin A.B. Localisation of stress-affected chemical reactions in solids described by coupled mechanics-diffusion-reaction models // International Journal of Engineering Science. – 2024. – Vol. 196. https://doi.org/10.1016/j.ijengsci.2023.10400
  44. Schneider D. On the stress calculation within phase-field approaches: a model for finite deformations // Computational Mechanics. – 2016. – Vol. 60. – pp.203-217. https://doi.org/10.1007/s00466-017-1401-
  45. Beckermann C. Modeling Melt Convection in Phase-Field Simulations of Solidification // Joirnal of Computational Physics. – 1999. – Vol. 154. – pp.468-496. https://doi.org/10.1006/jcph.1999.6323
  46. Wilmanski K. Thermomechanics of continua. Springer-Verlag. Berlin, Heidelberg, New York. 1998. 273

Statistics

Views

Abstract - 7

PDF (Russian) - 5

Cited-By


PlumX


Copyright (c) 2025 Freidin A.B., Mikheev E.V.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies