УДК 53.087, 681.2.08
Проблемы измерения потоков энергии излучения на объектах транспортной инфраструктуры с излучающей плазмой
В. Ф. Лапшин
Петербургский государственный университет путей сообщения Императора Александра I, Российская Федерация, 190031, Санкт-Петербург, Московский пр., 9
Для цитирования: Лапшин В. Ф. Проблемы измерения потоков энергии излучения на объектах транспортной инфраструктуры с излучающей плазмой // Известия Петербургского университета путей сообщения. — СПб.: ПГУПС, 2022. — Т. 19. — Вып. 4. — С. 685-692. DOI: 10.20295/1815-588Х-2022-4-685-692
Аннотация
Цель: Объекты транспортной инфраструктуры включают в себя многочисленные технические устройства, содержащие излучающую плазму. К их числу относятся источники видимого и спектрального излучения, коммутационные аппараты на тяговых подстанциях, высокотемпературные теплообменники и камеры сгорания. При экспериментальном исследовании таких устройств проводят измерения мощности и спектрального состава излучения, испускаемого плазмой. Для этого, как правило, используют фотодиоды небольших размеров, устанавливаемые на некотором расстоянии от плазменного образования. Цель настоящей работы состоит в установлении связи между мощностью излучения плазмы и величиной потока излучения, падающего на рабочую поверхность фотодиода. Методы: Для решения поставленной задачи в работе используется метод прямого интегрирования уравнения переноса излучения в предположении однородной структуры плазмы, но в отсутствии предположения о наличии локального термодинамического равновесия. Рассматриваются случаи отражающих и поглощающих поверхностей, ограничивающих плазму. Результаты: Найдены явные выражения для потока излучения, выходящего с поверхности плазменного образования, и потока, падающего на поверхность фотодиода. Численно исследована зависимость отношения величин этих потоков излучения от геометрических размеров и оптической толщины плазмы. Для случая отражающих поверхностей, ограничивающих плазму, найдено простое асимптотическое выражение для величины отношения потоков и определена область его применимости. Практическая значимость: Соотношения, устанавливающие связь между мощностью излучения плазмы и величиной мощности потока излучения, падающего на рабочую поверхность фотодиода, позволяют решить основную задачу экспериментального исследования плазменных образований — восстановления характеристик плазмы по результатам измерения фототока. Результаты, полученные в работе, могут быть использованы при экспериментальном исследовании технических устройств, содержащих излучающую плазму.
Ключевые слова: Измерение потока излучения, излучающая плазма, уравнение переноса излучения.
Введение
Развитие объектов транспортной инфраструктуры требует решения большого числа технических задач. Одной из них является задача обеспечения надежной работы технических устройств, содержащих излучающую плазму. К таким устройствам относятся газоразрядные источники интенсивного видимого и селективного излучения, используемые для освещения объектов железнодорожной инфраструктуры, для обеззараживания воды в системах водоснабжения, для отвердевания строительных смесей в транспортном и гражданском строительстве и других целей. Излучающая плазма используется в высокотемпературных теплообменниках и камерах сгорания. Интенсивно излучающие дуги возникают в сильноточных электрических цепях: контактной сети электрифицированных железных дорог, в выключателях высокого напряжения, в коммутационных аппаратах на тяговых подстанциях и других устройствах [1-4]. В соответствии с широким применением технических устройств, включающих в себя излучающую плазму, ведется непрерывная работа по исследованию возможности улучшения их характеристик [5-7]. Центральное место в этих исследованиях занимает экспериментальное изучение. Плазменные образования всегда излучают в окружающее пространство значительную часть вложенной в них энергии. Измерения интенсивности и спектрального состава этого излучения позволяют получать разнообразную информацию о процессах, протекающих в технических устройствах. Для измерения мощности и спектра излучаемых плазмой электромагнитных волн используются фотодетекторы [8, 9]. Основной проблемой здесь является установление связи между энергией, поглощаемой фотодиодом, и энергией, излучаемой плазменным образованием. Сложность решения этой задачи связана с тем, что плазменный источник
излучения является протяженным и объемным, имеет конечную оптическую толщину и, кроме того, ограничен поверхностью, которая может быть отражающей или поглощающей. Цель настоящей работы состоит в том, чтобы установить связь между мощностью излучения плазменного образования в техническом устройстве и мощностью излучения, поглощаемого фотодиодом.
Решение уравнения переноса излучения
В работе рассматривается характерная для технических устройств цилиндрическая геометрия плазменного образования. Кроме того, поскольку в целом ряде случаев [9] свечение плазмы является однородным по объему, то поставленная задача решается в рамках модели однородной плазмы. Геометрия задачи изображена на рис. 1. Ось 2 совмещена с осью плазменного цилиндра радиуса Я, основания которого имеют координаты г = -к/2 и г = к/2. Фотодиод с небольшой по площади поверхностью ^^ располагается перпендикулярно радиусу плазменного цилиндра на расстоянии Ь от его оси. При этом фотодиод расположен напротив плазменного образования так, что координата его геометрического центра находится в пределах -к/2 < < к/2. Размеры рабочей поверхности фотодиода обычно малы по сравнению с размерами излучающей плазмы. По этой причине поверхность 8й фотодиода можно считать освещенной равномерно. Мощность Рй потока излучения, падающего на фотодиод, задается в этом случае соотношением [10]:
Pd(X) = Sd J (Q-n)Ix(Q)
d Q,
(1)
ЛП
где ¡х — спектральная интенсивность излучения, падающего на поверхность фотодиода; АО — телесный угол, в пределах которого фотоны достигают поверхности фотодиода.
d
Рис. 1. Геометрия плазменного столба: ^ — вектор единичной длины в направлении полета фотона; ег — вектор единичной длины в радиальном направлении; п — нормаль к поверхности фотодиода
Отметим, что полярный и азимутальный углы 0d и yd вектора Q, задающего направление полета фотона на поверхности фотодиода Sd, связаны с углами 0 и у этого же вектора на поверхности Spl плазменного цилиндра (см. рис. 1) соотношениями:
Lsin0d = Rsin0 и yd = y. (2)
Как видно из рис. 1, телесный угол AQ определяется возможными значениями азимутального
и полярного углов 0d и У: -0d max < 0d < 0d max и
ymin < у < ymax. Предельные значения углов задаются при этом соотношениями (см. рис. 1):
sin 0„
■R / L.
(3)
tg У min
0,5h + z,
L cos 0d + R cos 0 '
tg У m
0,5h - z,
L cos 0d + R cos 0
(4)
(5)
Здесь верхний знак в знаменателе соответствует случаю, когда поверхности, ограничивающие плазму сверху и снизу, являются отражающими (коэффициент отражения а = 1), то есть когда фотоны могут прилетать на поверхность фотодиода ^^ из любой точки на боковой поверхности плазмы 8р1. Нижний знак соответствует случаю поглощающих поверхностей, ограничивающих плазму (коэффициент отражения а = 0). При этом основной вклад в поток излучения на диод вносят только те лучи, которые выходят из
точек на боковой поверхности цилиндра, проти- где Ь^, 0, у) = 27соз0 для а = 1. Для а = 0 воположной от диода.
Для определения спектральной интенсивности ¡х в работе используется стационарное уравнение переноса излучения в однородной плазме [11, 12]:
L(z, 0, у) = 2Rcos0 при < у < y2, L(z, 0, y) = = (0,5h + z)ctg |y| при -n/2 < y < -^j , L(z, 0, y) = = (0,5h - z)ctg |y| при y2 < y < п/2. Здесь tg^ = = -(0,5h + z)/2Rcos0 и tgy2 = (0,5h - z)/2Rcos0, -0,< z < 0,5^.
QW, (r, z, Q) = w,- k, I, (r, z, Q),
(6)
где ux — спектральная испускательная способность единицы объема плазмы в единичный телесный угол; r — радиальная переменная; к'х—коэффициент поглощения плазмы с учетом поправки на вынужденное излучение [13, 14].
Отметим здесь, что поглощением и излучением среды, окружающей плазменный объем, можно пренебречь. В этом случае, как следует из (6), величина спектральной интенсивности ¡х во внешней среде остается постоянной вдоль каждого направления: ¡х(0) = const. Соответственно, задача определения мощности излучения, падающего на фотодиод, сводится к определению спектральной интенсивности излучения ¡х, выходящего с поверхности плазменного образования.
Для определения ¡х на поверхности плазмы в работе решается уравнение переноса излучения (6). При постановке граничных условий на боковой поверхности Spl плазменного цилиндра считается, что излучение не входит в плазму через ее боковую поверхность: ¡x(Q)|Spl = 0 при Q ■ er < 0. При постановке граничных условий на верхней и нижней поверхностях, ограничивающих плазменный столб, рассматриваются два случая: случай полностью отражающих (а = 1) и полностью поглощающих (а = 0) поверхностей. В этих условиях решение (6) имеет вид:
I. (z. е. < = £
1-
exp
k,L(z, е, у) cos у
(7)
Вывод основных соотношений
Расчетную формулу для мощности Pd потока излучения, падающего на фотодиод, получим, подставляя в исходное выражение (1) решение уравнения переноса излучения (7):
.. 0d max Vmax
Pd(X) = Sd J dQd cos0d J dycos2 yx
1 - exp
k
0
' k,l(z, е, у)л
cos у
/-I
(8)
Как видно из (4) и (5), пределы интегрирования в (8) имеют величину порядка ~ утах ~ к/Ь. При дальнейшем рассмотрении учтем, что в экспериментах обычно располагают фотодиод на расстоянии существенно большем, чем размеры плазменного образования, т. е. что к/Ь << 1. По этой причине во внутреннем интеграле в (8) с точностью до членов порядка (к/Ь)2, можно положить соэу ~ 1. Тогда, переходя с помощью (2) к интегрированию по переменной 0, получаем:
(,)=п к,
Здесь
(9)
п/2
f%= 4п-1 (1 - b2)-1 J de cos 0[1 - exp(-TD cos 9)]
0
(l - b2 sin2 e)1/2 ± bcos 0
где td = 2Rk\ — оптическая толщина плазменного столба;
x
Ь = Я/Ь. Знак «±» в последнем выражении по-прежнему означает два случая граничных условий: знак «+» соответствует случаю а = 1, а знак «-» — случаю а = 0.
Теперь найдем мощность потерь энергии на излучение Рр1 через боковую поверхность плазменного столба 8рГ Величина Рр1 выражается через спектральную интенсивность излучения на этой поверхности:
P
S
рр1 (Я) = J ds J (й e) ix (z, й)
d а.
SP1 (2П)
Используя здесь для величины ¡х выражение (7), получаем:
PSpl = nSpi U' £Я: kl
(10)
где 8р1 = 2пЯк — площадь излучающей поверхности плазмы. Величина 8Х имеет смысл спектральной степени черноты и для случая отражающих поверхностей, ограничивающих плазму (а = 1), имеет вид
п/2
п/2
£я (td) = 1 — J d0 cos 0 J dу cos2 у exp x n 0
(1 1)
f n Л
td cos 0
cos у
— = Y (T D, R / L, 5) Ppl 4nL2^ D ' J
(13)
где у = /х/ех — безразмерная функция, зависящая от оптической толщины тв и геометрических параметров Я/Ь и 5 = к/4Я.
Теперь отметим, что в экспериментах по исследованию излучения плазменных образований фотодетекторы располагают, как правило, на значительных расстояниях от плазмы. Сама плазма, имея небольшие размеры, в большинстве случаев является оптически прозрачной. Это позволяет, используя малость величин Я/Ь и тп, разложить подынтегральные функции в (11) и (12) в ряды по этим параметрам. Соответствующие интегралы вычисляются и в случае а = 1 приводят к явному выражению для у:
Y
С0 C1T D + С2 T D
(1 -2td/3 + 0,125td(a-lntd))'
(14)
Здесь с0 = 1 + 8Ь/3п + 7Ь2/8, с1 = 4(1 + 9пЬ/32 + + 0,9Ь2), с2 = 0,125(1 + 128Ь/45 + 11Ь2/12), а = 1,25 + + 21п2 - СЕ = 2,059079, СЕ = 0,5772157 — постоянная Эйлера. Отметим, что при выводе выражения (13) пренебрегалось членами порядка (к/Ь)2, а при записи (14), кроме того, опущены величины порядка (Я/Ь)4 и (тп)3.
а для поглощающей поверхности (а = 0):
1 п/2 п/2
£я(td,5) = — Jdx J d0cos0 J dycos2 yx
П 0 0 -n/2
x[1 - exp (-td£(x, 0, y))]. (12)
Здесь 5 = к/4Я, х = г/(0,5й), g = 5(1 + х)/81п(|^|) при -п/2 < у < -у1 , g = cos0/cosy при -у1 < у < у2 и g = 5(1 - х)^ту при у2 < у < п/2.
Связь между мощностью Рр, излученной плазменным образованием, и мощностью Рй, поглощаемой фотодиодом, задается соотношением:
Результаты вычислений и обсуждение
Приступая к обсуждению результатов, прежде всего отметим, что величина $>й / 4пЬ2 в (13) представляет из себя зависимость отношения Рй / Рр1 от расстояния Ь для простейшего случая точечного изотропного источника. Соответственно, величина у в (13) характеризует отличие этой зависимости от простейшего случая, обусловленное конечностью геометрических размеров и оптической толщины плазмы, а также свойствами поверхностей, ограничивающих плазму. На рис. 2 приведены результаты расчета зави-
pi
3,5-1 3,02,52,01,51,0
1,8 1,6 1,4 1,2
0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 R / L
0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 R / L
Рис. 2. Зависимость величины у от отношения Я/Ь в случае отражающих (а) и поглощающих (б) оснований плазменного цилиндра при разных значениях оптической толщины столба плазмы т0:
1 — 0,1; 2 — 0,5; 3 — 1,0; 4 — 1,5; 5 — 2,5
симости величины у = ^ / вх от геометрического параметра Я/Ь для случаев, когда поверхности, ограничивающие плазменный цилиндр сверху и снизу, являются отражающими (рис. 2, а, а = 1) и поглощающими (рис. 2, б, а = 0, к = Я). Хорошо видно, что величина у слабо зависит от оптической толщины столба плазмы тд: при изменении тп в 25 раз (от 0,1 до 2,5) у меняется в пределах 20 %. Также отметим здесь, что характер зависимости от расстояния Ь для случаев а = 0 и а =1 существенно отличается. По мере удаления фотодиода от излучающей плазмы, при а = 1 (отражающие основания плазменного цилиндра) отношение Рй / Рр1 стремится к значению, соответствующему точечному источнику. При а = 0 (поглощающие основания), наоборот, при удалении от плазменного излучателя величина Рс1 / Рр1 все более отличается от значения для точечного источника.
Заключение
В настоящей работе рассмотрена задача, возникающая при измерении потерь энергии на излучение в технических устройствах, содержащих излучающую плазму. Решено уравнение переноса излучения в условиях, характерных для экспериментального исследования излучательных
свойств плазмы. Получены соотношения, устанавливающие связь между мощностью излучения плазмы и величиной мощности потока излучения, падающего на рабочую поверхность фотодиода. Выполнены расчеты зависимости отношения этих мощностей от геометрических размеров плазменного образования и его оптической толщины. Для случая отражающих поверхностей, ограничивающих плазму, получена простая асимптотическая формула. Результаты, полученные в работе, могут быть использованы при экспериментальном исследовании технических устройств, содержащих излучающую плазму.
Библиографический список
1. Райзер Ю. П. Физика газового разряда / Ю. П. Рай-зер. — Долгопрудный: Издательский Дом «Интеллект», 2009. — 736 с.
2. Lapshin V F. Eco-friendly sun lamp for railway facilities / V. F. Lapshin // J. Phys.: Conf. Ser. — 2021. — Vol. 2131(4). — P. 042092.
3. Zabello K. K. Anode Surface State and Anode Temperature Distribution after Current Zero for Different AMF-Contact Systems / K. K. Zabello, I. N. Poluyanova, A. A. Logachev et al. // IEEE Trans. Plasma Sci. — 2019. — Vol. 47. — № 8(1). — Pp. 3563-3571.
4. Logachev A. A. Cathode Surface State and Cathode Temperature Distribution after Current Zero of Different
AMF-Contacts / A. A. Logachev, I. N. Poluyanova, K. K. Zabello et al. // IEEE Trans. Plasma Sci. — 2019. — Vol. 47. — № 8(1). — Pp. 3516-3524.
5. Wang L. J. Simulation results of influence of constricted arc column on anode deformation and melting pool swirl in vacuum arcs with AMF contacts / L. J. Wang, X. Zhang, X. Huang et al. // Phys. Plasmas. — 2017. — Vol. 24(11). — Р. 113511
6. Li T. Simulation of Electric Arc Characteristics Based on MATLAB/Simulink / T. Li, H. Li // IOP Conf. Ser.: Mater. Sci. Eng. — 2018. — Vol. 452(4). — Р. 042080.
7. Litvinova V V Stochastic model of thermal processes in the contact network at arc discharges occurring at high speeds of movement / V. V. Litvinova, V. I. Moiseev, E. V Runev // CEUR Workshop Proceedings. — 2020. — Vol. 2803. — Pp. 84-91.
8. Poluyanova I. N. Measurements of Thermal Radiation Brightness of Anode Surface After Current Zero for a Range of Current Levels / I. N. Poluyanova, K. K. Zabello, A. A. Logatchev et al. // IEEE Trans. Plasma Sci. — 2017. — Vol. 45. — № 8(2). — Pp. 2119-2125.
9. Баринов Ю. А. Мощность излучения сильноточной вакуумной дуги, стабилизированной аксиальным магнитным полем, в видимой и ультрафиолетовой областях спектра / Ю. А. Баринов, К. К. Забелло, А. А. Логачев и
др. // Письма в ЖТФ. — 2021. — Т. 47. — Вып. 3. — С. 18-20
10. Lapshin V. F. Radiative heat transfer in plasma of pulsed high pressure caesium discharge / V. F. Lapshin // J. Phys.: Conf. Ser. — 2016. — Vol. 669. — P. 012035.
11. Зельдович Я. Б. Физика ударных волн и высокотемпературных гидродинамических явлений / Я. Б. Зельдович, Ю. П. Райзер. — М.: Физматлит, 2008. — 656 с.
12. Abaszadeh M. Analysis of radiative heat transfer in two-dimensional irregular geometries by developed immersed boundary-lattice Boltzmann method / M. Abaszadeh, A. Safavinejad, A.A. Delouei et al. // JQSRT. — 2022. — Vol. 280. — P. 108086.
13. Baksht F. G. Modeling of the Plasma Waveguide on the Basis of the Pulse-Periodic High-Pressure Cesium Discharge / F.G. Baksht, VF. Lapshin // Plasma Phys. Rep. — 2020. — Vol. 46. — Pp. 846-849.
14. Loginov A. V. Radiative Constants in the W VII Ion Spectrum / A. V Loginov, V I. Nikitchenko // Optics and Spectroscopy. — 2020. — Vol. 128(8). — Pp. 1078-1081.
Дата поступления: 21.06.2022 Решение о публикации: 08.08.2022
Контактная информация:
ЛАПШИН Владимир Федорович — д-р физ.-мат. наук, проф.; [email protected]
Problems of Radiation Energy Flux Measurement at Transport Infrastructure Facilities with Radiating Plasma
V. F. Lapshin
Emperor Alexander I Petersburg State Transport University, 9, Moskovsky pr., Saint Petersburg, 190031, Russian Federation
For citation: Lapshin V. F. Problems of Radiation Energy Flux Measurement at Transport Infrastructure Facilities with Radiating Plasma // Proceedings of Petersburg Transport University, 2022, vol. 19, iss. 4, pp. 685-692. (In Russian). DOI: 10.20295/1815-588X-2022-4-685-692
Summary
Purpose: Transport infrastructure objects comprise numerous technical devices containing radiating plasma. They include visible and spectral radiation sources, switching devices at traction substations, high-temperature heat exchangers and combustion chambers. During experimental study of such devices, the measurements of power and spectral composition of radiation emitted by plasma are carried out. For this, as a rule, photodiodes of small sizes are used, installed at a certain distance from plasma formation. The purpose of the present work is to establish relationship between plasma emission power and radiation flux amount incident on photodiode working surface. Methods: To solve the set task, the method of direct integration of radiation transfer equation
in homogeneous plasma structure assumption, but in the absence of local thermodynamic equilibrium presence assumption, is used. The cases of reflecting and absorbing surfaces, limiting a plasma, are considered. Results: Explicit expressions are found for radiation flux which exits plasma formation surface and for flux incident on photodiode surface. The dependence for the ratio of the values of these radiation fluxes to plasma geometric sizes and optical thickness is numerically studied. For the case of reflecting surfaces that bound plasma, simple asymptotic expression is found for ratio value for the fluxes and applicability field of the expression is determined. Practical significance: The ratios, which establish relationship between plasma radiation power and power value for radiation flux incident on photodiode working surface, make it possible to solve the main task of the experimental study of plasma formations — the restoration of plasma characteristics according to photocurrent measurement results. Obtained in the work results can be used at experimental study of technical devices containing emitting plasma.
Keywords: Measurement of radiation flux, radiating plasma, radiative transfer equation.
References
1. Rayzer Yu. P. Fizika gazovogo razryada [Physics of gas discharge]. Dolgoprudnyy: «Intellekt» Publ., 2009. 736 p. (In Russian)
2. Lapshin V F. Eco-friendly sun lamp for railway facilities. J. Phys.: Conf. Ser. 2021, vol. 2131(4), pp. 042-092.
3. Zabello K. K., Poluyanova I. N., Logachev A. A., Begal D. I., Shkol'nik S. M. Anode Surface State and Anode Temperature Distribution after Current Zero for Different AMF-Contact Systems. IEEE Trans. Plasma Sci. 2019, vol. 47, I. 8 (1), pp. 3563-3571.
4. Logachev A. A., Poluyanova I. N., Zabello K. K., Barinov Yu. A., Shkol'nik S. M. Cathode Surface State and Cathode Temperature Distribution after Current Zero of Different AMF-Contacts. IEEE Trans. Plasma Sci. 2019, vol. 47, I. 8 (1), pp. 3516-3524.
5. Wang L. J., Zhang X., Huang X., Jia S. L. Simulation results of influence of constricted arc column on anode deformation and melting pool swirl in vacuum arcs with AMF contacts. Phys. Plasmas. 2017, vol. 24 (11), pp. 113-511.
6. Li T., Li H. Simulation of Electric Arc Characteristics Based on MATLAB/Simulink. IOP Conf. Ser.: Mater. Sci. Eng. 2018, vol. 452 (4), p. 042080.
7. Litvinova V V, Moiseev V I., Runev E. V Stochastic model of thermal processes in the contact network at arc discharges occurring at high speeds of movement. CEUR Workshop Proceedings. 2020, vol. 2803, pp. 84-91.
8. Poluyanova I. N., Zabello K. K., Logatchev A. A., Yakovlev V V., Shkol'nik S. M. Measurements of Thermal Radiation Brightness of Anode Surface After Current Zero for a Range of Current Levels. IEEE Trans. Plasma Sci. 2017, vol. 45, I. 8(2), pp. 2119-2125.
9. Barinov Yu. A. Moshchnost' izlucheniya sil'notochnoy vakuumnoy dugi, stabilizirovannoy aksial'nym magnitnym polem, v vidimoy i ul'trafioletovoy oblastyakh spektra [Radiation power of a high-current vacuum arc stabilized by an axial magnetic field in the visible and ultraviolet regions of the spectrum]. Pis'ma v ZhTF [Letters to ZhTF]. 2021, vol. 47, I. 3, pp. 18-20. (In Russian)
10. Lapshin V F. Radiative heat transfer in plasma of pulsed high pressure caesium discharge. J. Phys.: Conf. Ser. 2016, vol. 669, pp. 012-035.
11. Zel'dovich Ya. B., Rayzer Yu. P. Fizika udarnykh voln i vysokotemperaturnykh gidrodinamicheskikh yavleniy [Physics of shock waves and high-temperature hydrody-namic phenomena]. Moscow: Fizmatlit Publ., 2008. 656 p. (In Russian)
12. Abaszadeh M., Safavinejad A., Delouei A. A., Amiri H. Analysis of radiative heat transfer in two-dimensional irregular geometries by developed immersed boundary-lattice Boltzmann method. JQSRT. 2022, vol. 280, pp. 108086.
13. Baksht F. G., Lapshin V F. Modeling of the Plasma Waveguide on the Basis of the Pulse-Periodic High-Pressure Cesium Discharge. Plasma Phys. Rep. 2020, vol. 46, pp. 846849.
14. Loginov A. V Radiative Constants in the W VII Ion Spectrum. Optics and Spectroscopy. 2020, vol. 128(8), pp. 1078-1081.
Received: June 21, 2022 Accepted: August 08, 2022
Author's information:
Vladimir F. LAPSHIN — Dr. Sci. in Physics and Mathematics, Professor; [email protected]