Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki
[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2021, vol. 25, no. 1, pp. 193-202 ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtu1831
Short Communications
MSC: 74F15, 74F05, 74C20
Healing of cracks in plates by strong electromagnetic field
© K. V. Kukudzhanov1, A. L. Levitin1, U. Kh. Ugurchiev2
1 Ishlinsky Institite for Problems in Mechanics, Russian Academy of Sciences, 101—1, pr. Vernadskogo, Moscow, 119526, Russian Federation.
2 Mechanical Engineering Research Institute of the Russian Academy of Sciences, 4, M. Khariton'evskii per., Moscow, 101990, Russian Federation.
Abstract
The problem of a pulsed high-energy electromagnetic field action on an edge crack in a thin plate, reproducing the pioneering experiment of Soviet scientists on the destruction of the crack tip by a strong electromagnetic field, is considered. The numerical simulation is based on the proposed electrothermomechanical model of a short-pulse high-energy electromagnetic field (HEMF) action on a material with a crack. The model takes the phase transformations (melting and evaporation) of the material occurring in the vicinity of defects and the corresponding changes in the rheology of the material in the areas of these transformations into account, as well as the possibility of electric current flowing between the free surfaces of the crack (breakdown due to electron emission). All physical and mechanical properties of the material are considered temperature-dependent. The model equations are coupled and solved together on a moving finite element grid using the arbitrary Euler-Lagrangian method. The processes of localization of the current density and temperature fields, phase transformations (melting and evaporation) at the crack tip, autoelectronic and thermoelectronic emissions between free crack surfaces, and the effect of these processes on crack healing are investigated. The simulation results are compared with the available experimental data on the pulse field action on the edge crack in the plate. The average metal heating rate, temperature gradients and time forming of the crater obtained in the vicinity of the crack tip are in good quantitative agreement with the experimental data. Away from the crack, as well as on the crack sides away from the tip, the temperature rises slightly. The process of modeling the electromagnetic field action, similar to the experiment,
Review Article
Q ©® The content is published under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/) Please cite this article in press as:
Kukudzhanov K. V., Levitin A. L., Ugurchiev U. Kh. Healing of cracks in plates by strong electromagnetic field, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2021, vol. 25, no. 1, pp. 193-202. https://doi.org/10.14498/vsgtu1831 (In Russian). Authors' Details:
Konstantin V. Kukudzhanov & https://orcid.org/0000-0001-9060-2838 Cand. Phys. & Math. Sci.; Senior Researcher; e-mail:kconstantin@mail.ru Alexander L. Levitin © https://orcid.org/0000-0003-2077-8808 Leading Programmer; e-mail: alex_lev@ipmnet.ru Umar Kh. Ugurchiev © https://orcid.org/0000-0003-2072-6354 Researcher; e-mail: umar77@bk.ru
© Samara State Technical University
193
was accompanied by melting at the crack tip, as well as metal evaporation. Thus, under the considered current action, a crater is formed at the crack tip, which prevents the further spread of the crack, leading to its healing. It was not possible to obtain similar results using the previously proposed models.
Keywords: pulsed electromagnetic action, defect healing, crack arrest, high-energy field.
Received: 11th October, 2020 / Revised: 15th February, 2021 / Accepted: 10th March, 2021 / First online: 23rd March, 2021
Introduction. An external electromagnetic field is applied to the conduc-
tive sample, which induces a current in the material with a current density of
108 to 1011 A/m2 during short time interval. A pulsed electromagnetic field is called a high-energy electromagnetic field (HEMF) if electromagnetic action during 10-4 — 10-3 s leads to amount of the scattering specific electromagnetic energy in the material within the range 108 ^ q ^ 1010 J/m3.
The presence of defects in the conductive solid leads to the localization (concentration) of the electric field in their vicinity. This causes a locally inhomoge-neous temperature distribution in the material, due to the increased dissipation of electromagnetic energy in the area of defects, or rather in their sharp tips. The use of this circumstance to explain the healing by short pulses of high-density current of cracks propagating in siliceous iron under the influence of pulsed HEMF was carried out for the first time in [1].
Experiments [1,2] show that HEMF action leads to a sharp temperature increase at the microcracks tips. As a result, melting, evaporation, microexplosion, and crater formation occur at the macrocracks tips.
Thus, the resulting crater prevents the further crack propagation and leads to crack healing.
In particular, these experiments showed that when short pulses of high-density current induced by the HEMF were applied to the edge crack in a thin plate in the vicinity of its tip, a sharp temperature increase was observed: the measured heating rate was 107 °C/s, the temperature gradients were 106-107 °C/m, and the temperature away from the crack did not exceed 10 °C. Such rates and localities of heat release led to melting and evaporation of the material at the tip of the crack, accompanied by a microexplosion, which ejected the explosion products from the crack tip in a direction perpendicular to the plate plane. At the same time, a crater formed at the crack tip [1].
A number of researchers have attempted to model the electrical, thermal, and mechanical processes that occur during the above-described pulsed action of the HEMF on a crack in a plate.
These processes were considered both in quasistatics [3-9] and in dynamics [11]. The plate material was assumed to be thermoelastic [7,9], thermo-rigid-plastic [10], thermoelastoplastic [11], and the deformations were small. The process of thermal conductivity is described by Fourier's law. In the works [3-5], the magnetic effects arising under the HEMF pulse action were also taken into account.
The simulation showed that an inhomogeneous temperature field occurs near the macrocrack and temperature localization reaches the melting point in the vicinity of the macrocrack tip. Compressive stresses occurring in the vicinity of the crack tip led not only to the crack arrest, but also to the convergence of its sides (partial closure of the defect). However, these models allowed to describe the process of heating and deformation only qualitatively: the heating rate, temperature gradients, and the time of the beginning of melting in the vicinity of the crack tip, obtained as a result of modeling, were significantly different from those observed in experiments [1,2]. Within the considered above models the temperature at the crack tip could not reach the evaporation temperature, therefore these models failed to describe the processes of explosion and the formation of a crater in the tip.
Due to mathematical difficulties in the analytical solution of electrothermoe-lastic (rigid-plastic) and numerical solution of electro-thermoelastic problems, none of the models [3-10] allowed (in the context of the assumption of small deformations) to trace the effects of the stress-strain state in the vicinity of the defect tip on the electric and temperature fields consistently. In addition, none of these models took into account the phase transformations in the metal and the resulting changes in rheology, including the dependence of the physical and mechanical properties of the metal on temperature after changes in the aggregate state of the substance during deformation. Therefore, the models under consideration did not allow coupling study of the determining processes occurring at the crack tip under the current pulse action.
Meanwhile, as shown in [1], the absolute temperature values (and its gradients) obtained in experiments are very large, they obviously exceed the melting point and reach the evaporation temperature of the metal. Under such conditions, it is necessary to take into account the changes in the aggregate state of the substance, the dependence of the properties of the metal on temperature and their influence on the stress-strain state.
To eliminate these lacks, we propose a quasi-stationary model of pulsed HEMF action on cracks in metal and a numerical method for solving the resulting system of equations.
On the basis of the developed model, the problem of the short-term HEMF action on the edge crack (with a rounded tip) in the plate is solved, which numerically reproduces the experiment performed in [1]. Changes in the electromagnetic and temperature fields, as well as in the aggregate states in the vicinity of the crack, are studied.
The study of these processes will allow us to understand the mechanism of healing of microcracks in plates under the HEMF action better and get closer to explaining the experimentally observed processes of healing fatigue cracks by welding their sides.
1. Electrothermomechanical model. Consider the dynamic response of a conducting plate to a HEMF pulse that induces a high-density electric current in the plate material for 10-4 s. There is an edge crack in the plate, as shown in Fig. 1. Our purpose is to reproduce the experiment carried out in the work [1] and to quantify the observed effects.
We formulate the assumptions and equations of the analytical model for describing the electrical, thermal and mechanical processes occurring in the plate
Figure 1. The sketch of the specimen 15x50x0.3 mm with the edge crack (the crack length is 7.5 mm, including the rounding at the tip with a radius of 10 pm, the crack has parallel sides and thickness is 20 pm); only the upper half of the sample is simulated in FEM due to the symmetry; jo(i) are electrical current density vector applied to plate sides
material. The characteristic time of each of these processes is approximately inversely proportional to the propagation velocity of the corresponding perturbations. If we neglect the electromagnetic induction and viscosity during the propagation of mechanical disturbances, then the time required to establish the electromagnetic and mechanical processes considered in this model is 10-13-10-12 s and 10-8-10-7 s, respectively. This is significantly less than the time of the external action of the HEMF source (10-4 s).
Therefore, to obtain the electric potential <p in a conducting material, we use the law of conservation of charge (assuming that the current in the sample is steady). Since we search for a numerical solution by the finite element method, we use it in the variational formulation (1) [12]. At the same time, the laws of Ohm and Joule-Lenz are considered valid for the material.
The displacement field u is determined from the equilibrium equations written in the form of the virtual work principle (2) [12]. The deformations are assumed to be finite. In this case, the additivity of the rates of elastic, plastic, and temperature deformations is assumed (2). For the velocities of elastic and plastic deformations, Hooke's law for an isotropic body (3) and the associated flow law with the Mises plasticity (4) condition are assumed to be valid, respectively. The rates of temperature deformations are assumed to be linearly related to the temperature derivative (5).
The analytical estimates show that during pulsed electromagnetic field action, the current density at the crack tip is such that the relative temperature change at the distances of the order of the free path of phonons (and electrons) is greater than one. It is also evidenced by the measured temperature gradients in the vicinity of the crack vertices, which turn out to be very large [1] It is not correct to apply the Fourier law of thermal conductivity in this case. In addition, the time of electromagnetic impact on the material is small, 10-4 s. Therefore, thermal conductivity should be neglected and the process of HEMF exposure should be considered adiabatic (6).
The temperature field T is determined from the law of conservation of energy. In this case, we take into account the heat released in the unit volume in the current configuration of the body per unit time (6) due to the flow of electric current in accordance with the Joule-Lenz law (7), the heat released during plastic
deformation (8), as well as the latent heat absorbed in the processes of melting and evaporation. Therefore, the resulting evolutionary equation for temperature should be extended with equations (9) and (10), which represent the conditions for the absorption of latent heat during the transition of a substance from one aggregate state to another: from solid to molten during melting and from liquid to gaseous during evaporation.
As noted above, during the current action on the crack in the plate, the metal temperature changed from room temperature to the evaporation temperature [1]. Therefore, all the physical and mechanical characteristics of the model (density, specific heat, electrical conductivity, coefficient of thermal expansion, elastic modulus, yield strength, etc.) are assumed to depend on the temperature over the entire range until the evaporation temperature is reached.
At the points where condition (9) is realized, the material is considered molten. In this case, there is a sharp change in all the physical properties of the material: electrical conductivity, heat capacity, density, coefficient of linear expansion and all other mechanical characteristics of the material. Such a change in the properties of the material corresponds to the available experimental data given in [14-16], which shows the dependence of the properties of various metals on temperature. In this case, there is a decrease in some (elastic modulus, yield strength, linear expansion coefficient, etc.) and an increase in others (density, etc.) of the physical characteristics of the material [15,16].
Thus, in the model under consideration, when the melting temperature was reached the material does not lose the ability to conduct an electric current [15,16] and after the metal completely passes into the molten state, the melt is heated further. As the result, elastic modulus decreases and the yield limit falls to zero which allows us to describe the behavior of the molten material by the same constitutive equations (3)-(5), which degenerate for the melt into thermoviscoelastic with nonlinear viscosity.
At the points where condition (10) is realized the material is considered to have completely evaporated. At all subsequent times, the current density j = 0, the stress tensor a = 0 and the temperature are assumed to be constant T = Tevap. Therefore, in the model under consideration, from the time when the material completely evaporates, the metal loses its ability to conduct an electric current, its further heating does not occur, it loses the properties of a viscous liquid and is considered as a rarefied gas.
The complete system of equations of the considered electrothermomechanical model has the form [12,13]:
j = aE(T) E = -aE(T)^, (1) -, é= ¿el + èpl + ¿th, (2)
& = X(T) èel : I + 2 ) èel
(3)
epl =À ^-=À s, <r = aY(T),
deth = a(T) I dT, p(T) c(T) T = rE + rpl + rmelt + reval
(5)
(6)
rE = rjE jE = r]E Vip-aE-Vip, (7)
rpl = a:£pl, (8)
T = Tmelt, ¿sol < t < iiiq, / ^ + rpl)di = P Amelt, (9)
^ sol /* ^ evap
T = Tevap, ieliq < t < ievap, / (rE + dt = p Aevap. (10)
^eliq
The notation in this equations is given in [12,13].
The boundary and contact conditions for the electrical potential displacements u and stresses a at the boundaries of the integration domain used in this model, as well as the initial conditions, coincide with those considered in [12,13].
2. Results of numerical simulation. Equations (1)-(10) of the electrother-momechanical problem coupled with boundary, contact, and initial conditions were solved by the finite element method for a plate with an edge crack in a 3D formulation. The calculations were done using 20-node quadratic brick finite elements.
Simulation was done for zinc samples, the physical and mechanical properties of which and their dependence on temperature were taken in accordance with [15,16].
The electric potential difference and its change from t used in the calculations were selected so that the induced current density coincided with the one measured in the experiment [1]. The pulse duration was 160 ps.
Fig. 2 a shows the stress field a33 in the vicinity of the crack tip at t = 29.5 ps. Fig. 2 a shows that compressive stresses occur at the crack tip, which prevent further crack propagation. This result corresponds to the results obtained using models [3-6,10], which also predicted the appearance of stresses compressing the crack sides at the crack tip. At the same time, calculations show that the residual plastic deformations in the vicinity of the crack tip reach tens of percent. This gives reason to believe that the restoration of the distance between the banks of the crack to the original size will not occur and after the relaxation of the temperature the sides will remain close together.
A sharp increase in temperature in the vicinity of the crack tip under HEMF pulse action was observed in the experiments presented in [1]. The measured
Figure 2. a) normal stress field <733 along vertical axis in the vicinity of the crack tip at t = 29.5 ps; b) temperature field in the vicinity of the crack tip at time t = 38 ps immediately after the forming of the crater (the gray area is the crater zone)
heating rate in the vicinity of the crack tip was on the order of 107 °C/s, the local temperature gradients were 6.0 ■ 106 - 5.0 ■ 107 C/m, and the temperature change did not exceed 10 C away from the crack. At the same time, melting and evaporation of the material was observed at the tip of the crack, accompanied by a microexplosion with the release of products from the crack tip in the direction perpendicular to the plate plane and the formation of a crater.
The results obtained by the proposed model are in good quantitative agreement with the experiment: the calculated average heating rate (in the vicinity of the crack tip with a radius of 10 pm) was 6.3 ■ 107 C/s, and the local temperature gradients were 107-108 C/m. Away from the crack, the temperature did not rise more than 10 C (at times exceeding 100 ps), nor did heating occur on the sides of the crack away from the tip.
The process was accompanied by melting at the crack tip, as well as evaporation of the metal. Melting begins at the crack tip at time t = 6.5 ps, and the melt moves in a direction perpendicular to the plate plane. Evaporation begins at the time t = 31 ps. The time of crater formation can be associated with an interval of 31-38 ps. In the experiment, the crater was formed at 38-45 ps.
Thus, under the considered current action, a crater is formed at the crack tip, which prevents the further spread of the crack, leading to its healing.
Fig. 2 b shows the temperature field in the vicinity of the crack tip at time t = 38 ps (immediately after the forming of the crater).
Calculations based on the proposed model showed that the autoelectronic (tunnel) emission between the free surfaces (sides) of the crack does not occur due to the insufficient potential difference for this under the considered HEMF action. While the thermoelectronic emission current between the crack banks can be neglected, since its density is significantly less than the current density induced by the HEMF in the vicinity of the crack tip.
3. Conclusions. Thus, there is reason to believe that the model reproduces the main features of electrothermomechanical processes in the vicinity of microde-fect correctly, describing the process of healing a crack and forming a crater.
Competing interests. We declare that we have no conflicts of interest in the authorship and publication of this article.
Authors' contributions and responsibilities. K.V. Kukudzhanov: Idea of study; Formulation of research goals and aims; Theoretical analysis; Analysis of calculation results; Supervision and consulting; Writing — original draft and review & editing. A.L. Levitin: Performing calculations, their analysis and verification; Writing — original draft and review & editing. U.Kh. Ugurchiev: Literature review; Analysis of experiments; Analysis of calculation results; Writing — part of the original draft. The authors are absolutely responsible for submitting the final manuscript in print. Each author has approved the final version of manuscript.
Funding. The study of K.V. Kukudzhanov is supported by the Russian Science Foundation (project no. 19-19-00616).
References
1. Finkel V. M., Golovin Yu. I., Sletkov A. A. Disintegration of a crack tip with a strong electromagnetic field, Sov. Phys. Dokl., 1977, vol.22, pp. 683-685.
2. Finkel V. M., Golovin Yu. I., Sletkov A. A. Possibility of braking rapid cracks by pulses of current, Sov. Phys. Dokl., 1976, vol.21, no. 4, pp. 216-218.
3. Kudryavtsev B. A., Parton V. Z., Rubinskii B. D. Electromagnetic and thermoelastic fields in a conducting plate with a cut of finite length, Mech. Solids, 1982, vol. 17, no. 1, pp. 110-118.
4. Parton V. Z., Kudryavtsev B. A., Rubinskii B. D. Crack propagation under the action of an electromagnetic field, Dokl. Akad. Nauk SSSR, 1980, vol. 250, no. 5, pp. 1096-1100 (In Russian).
5. Cai G. X., Yuan F. G. Stresses around the crack tip due to electric current and self-induced magnetic field, Adv. Eng. Software, 1998, vol.29, no. 3-6, pp. 297-306. https://doi.org/ 10.1016/S0965-9978(97)00078-1.
6. Cai G. X., Yuan F. G. Electric current-induced stresses at the crack tip in conductors, Int. J. Fract., 1999, vol.96, no. 3, pp. 279-301. https://doi.org/10.1023/A: 1018670829822.
7. Liu T. J. C. Effects of temperature-dependent material properties on stress and temperature in cracked metal plate under electric current load, World Academy of Science, Engineering and Technology. Int. J. Mechanical and Mechatronics Eng., 2010, vol.4, no. 5, pp. 474-479. https://doi.org/10.5281/zenodo.1072134.
8. Yu J., Zhang H., Deng D., Hao S., Iqbal A. Numerical calculation and experimental research on crack arrest by detour effect and joule heating of high pulsed current in remanufacturing, Chin. J. Mech. Eng., 2014, vol.27, no. 4, pp. 745-753. https://doi.org/10.3901/CJME. 2014.0414.075.
9. Gallo F., Satapathy S., Ravi-Chandar K. Melting and crack growth in electrical conductors subjected to short-duration current pulses, Int. J. Fract., 2011, vol. 167, no. 2, pp. 183-193. https://doi.org/10.1007/s10704-010-9543-0.
10. Ovchinnikov I.V. Influence of electric current on the plasticity of metals, PhD Thesis. Moscow, 1989, 123 pp. (In Russian)
11. Kukudzhanov K. V., Levitin A. L. Deformation processes of elastoplastic material with defects under electrodynamic loading, PNRPU Mechanics Bulletin, 2015, no. 1, pp. 106-120 (In Russian). https://doi.org/10.15593/perm.mech/2015.1.07.
12. Kukudzhanov K.V. Modeling the treatment of high-energy pulsed electromagnetic field of the micro-cracks in a polycrystalline metal, PNRPU Mechanics Bulletin, 2015, no. 4, pp. 138-158 (In Russian). https://doi.org/10.15593/perm.mech/2015.4.09.
13. Kukudzhanov K. V., Levitin A. L. Modeling the healing of microcracks in metal stimulated by a pulsed high-energy electromagnetic field. Part I, Int. J. Nanomech. Sci. Tech., 2015, vol.6, no. 3, pp. 233-249. https://doi.org/10.1615/NanomechanicsSciTechnolIntJ.v6. i3.60.
14. Sivukhin D. V. Termodinamika i molekulyarnaya fizika [Thermodynamics and Molecular Physics]. Moscow, Fizmatlit, 2014, 544 pp. (In Russian)
15. Pikunov M. V. Plavka metallov. Kristallizatsiia splavov. Zatverdevanie otlivok [Metal Smelting. Alloy Crystallization. The Solidification of Castings]. Moscow, MISiS, 1997, 374 pp. (In Russian)
16. Pikunov M. V. Metallurgiia rasplavov. Kurs lektsii [Metallurgy of the Melts. A Course of Lectures]. Moscow, MISiS, 2005, 286 pp. (In Russian)
Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2021. Т. 25, № 1. С. 193-202
ISSN: 2310-7081 (online), 1991-8615 (print)
d https://doi.org/10.14498/vsgtu1831
УДК 517.958:531-133
Залечивание трещин в пластинах сильным электромагнитным полем
© К. В. Кукуджанов1, А. Л. Левитин1, У. Х. Угурчиев2
1 Институт проблем механики им. А. Ю. Ишлинского РАН, Россия, 119526, Москва, просп. Вернадского, 101, корп. 1.
2 Институт машиноведения им. А. А. Благонравова РАН, Россия, 101990, Москва, М. Харитоньевский пер., 4.
Аннотация
Рассматривается задача о воздействии импульсным высокоэнергетическим электромагнитным полем на краевую трещину в тонкой пластине, воспроизводящая пионерский эксперимент советских ученых по разрушению вершины трещины сильным электромагнитным полем. Численное моделирование осуществляется на основе предложенной электромеханической модели воздействия короткоимпульсным высокоэнергетическим электромагнитным полем на материал с трещиной. Модель учитывает фазовые превращения (плавление и испарение) материала, происходящие в окрестности дефектов, и соответствующие изменения реологии материала в областях этих трансформаций, а также возможность протекания электрического тока между свободными поверхностями трещины (пробоя за счет эмиссии электронов). Все физико-механические характеристики материала считаются зависящими от температуры. Уравнения модели связаны и решаются совместно на подвижной конечно-элементной сетке с применением смешанного метода Эйлера-Лагран-жа. Исследуются процессы локализации полей плотности тока и температуры, фазовых превращений (плавления и испарения) в вершине трещины, автоэлектронной и термоэлектронной эмиссии между свободными поверхностями трещины и влияние этих процессов на залечивание трещины. Проводится сравнение результатов моделирования с имеющимися экспериментальными данными по воздействию импульсного поля на краевую трещину в пластине. Полученные в окрестности вершины трещины средняя скорость нагрева металла и градиенты температуры неплохо количественно согласуются с экспериментальными данными.
Обзор
3 ©® Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru) Образец для цитирования
Кукуджанов К. В., Левитин А. Л., Угурчиев У. Х. Залечивание трещин в пластинах сильным электромагнитным полем // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2021. Т. 25, № 1. С. 193-202. https://doi.org/10.14498/vsgtu1831. Сведения об авторах
Константин Владимирович Кукуджанов А https://orcid.org/0000-0001-9060-2838 кандидат физико-математических наук; старший научный сотрудник; e-mail: kconstantin@mail. ru
Александр Леонидович Левитин © https://orcid.org/0000-0003-2077-8808 ведущий программист; e-mail: alex_lev@ipmnet.ru
Умар Хажбикарович Угурчиев © https://orcid.org/0000-0003-2072-6354 научный сотрудник; e-mail: umar77@bk.ru
© Самарский государственный технический университет
201
Вдали от трещины, а также на берегах трещины вдали от вершины температура поднималась незначительно. Процесс моделирования воздействия электромагнитным полем, аналогично эксперименту, сопровождается плавлением в вершине трещины, а также испарением металла. Таким образом, при рассматриваемом воздействии током в вершине трещины формируется кратер, который препятствует дальнейшему распространению трещины, приводя к ее залечиванию. Получить аналогичные результаты с помощью ранее предложенных моделей не удавалось.
Ключевые слова: импульсное электромагнитное воздействие, залечивание дефектов, торможение трещины, высокоэнергетическое поле.
Получение: 11 октября 2020 г. / Исправление: 15 февраля 2021 г. / Принятие: 10 марта 2021 г. / Публикация онлайн: 23 марта 2021 г.
Конкурирующие интересы. Заявляем, что в отношении авторства и публикации этой статьи конфликта интересов не имеем.
Авторский вклад и ответственность. К.В. Кукуджанов — идея исследования, формулировка целей и задач исследования, теоретический анализ, анализ результатов расчетов, руководство и консультирование, черновик и чистовик рукописи. А.Л. Левитин — выполнение расчетов, анализ, визуализация и верификация результатов расчетов, черновик и чистовик рукописи. У.Х. Угурчиев — анализ литературы, анализ экспериментов, анализ результатов расчетов, черновик части рукописи. Авторы несут полную ответственность за предоставление окончательной рукописи в печать. Окончательная версия рукописи была одобрена всеми авторами. Финансирование. Исследование К.В. Кукуджанова выполнено при поддержке Российского научного фонда (проект № 19-19-00616).