Скорость образования дефектов кристаллической решетки рельсов электромагнитного ускорителя под воздействием плазменного образования

Язык труда и переводы:
УДК:
539.1.043
Дата публикации:
21 ноября 2022, 00:28
Категория:
Математическое моделирование физических процессов и технических систем
Авторы
Маслов Артемий Глебович
МГТУ им. Н.Э. Баумана
Горев Владимир Васильевич
НИЦ «Курчатовский институт» — ИРЕА
Аннотация:
Построена математическая модель воздействия потока ионов плазмы приэлектродной области электромагнитного ускорителя на кристаллическую структуру рельсов. Показано, что в локальной области размером несколько нанометров скорость разрушения кристаллической решётки может достигать до 90 % узлов за 1 мкс. Также в работе рассчитано изменение удельного сопротивления материала рельсов в области повреждения.
Ключевые слова:
рельсотрон, ионное легирование, радиационные дефекты, эрозия
Основной текст труда

Электромагнитный запуск объектов применяется для ускорения макроскопических тел [1] до скоростей, превышающих 2 км/с, ускорения плазмы [2, 3] и др. Однако для рельсотронов одной из серьезных проблем, ограничивающих скорость метания, является эрозия канала ускорителя [2, 3]. В качестве причин эрозии рассматривается тепловое воздействие протекающего тока и ускоряющей плазмы [2, 3] и механическое воздействие [4].

Тепловая эрозия в рельсотроне имеет свои особенности по сравнению с системами, в которых движения плазменной дуги не происходит. Поскольку в этом случае эрозионная масса уносится вдоль канала ускорителя, в зоне катодного пятна понижена концентрация и повышена длина свободного пробега частиц плазмы [5]. При этом в приэлектродном слое плазмы наблюдаются сильные электрические поля (до 109 В/м), что приводит к ускорению ионов до значительных энергий, способных вызывать образование радиационных дефектов по Шоттки и по Френкелю, следствием чего будет изменение физических и механических свойств материала в приповерхностном слое.

Схожий процесс известен как ионное легирование материалов. В рамках этого процесса движение одного высокоэнергетичного иона может приводить к выбиванию из узлов кристаллической решетки сотен атомов мишени прямыми и непрямыми столкновениями. Тем не менее зачастую междоузлия и вакансии, созданные каскадом от первично-выбитого атома, находятся в локальной области столкновения [6]. В этом случае энергия, выделенная в такой области, может быть использована для оценки количества дефектов, согласно теории образования радиационных дефектов Кинчина — Пиза:

{{N}_{def}}={}^{W}/{}_{2{{W}_{d}}},

где W – начальная энергия иона, Дж; Wd – энергия образования радиационного дефекта, Дж; Ndef — общее число дефектов.

В данной работе использована теория ионного торможения Линхарда — Шарфа — Шиотта (ЛШШ) [7], теория образования радиационных дефектов Кинчина — Пиза [6] и распределения ионов плазмы в приэлектродной области (fi(W)). Также для расчета изменения проводимости материала использованы соотношения, полученные Джонсом и Моттом [8].

 Для математического моделирования удобно выразить начальную энергию иона как обратную функции глубины его проникновения:

W=x{{\left(W\right)}^{-1}}={\frac {{{a}^{2}}{{\left({\text{LambertW}}\left(-1,-{{\text{e}}^{-{\frac {50x\ {{b}^{2}}+a}{a}}}}\right)+1\right)}^{2}}}{{b}^{2}}},

где x – глубина проникновения иона, м; LambertW — специальная функция Ламберта, а коэффициенты a, Дж/м, иb, Дж–1/2/м, заданы следующими соотношениями:

{\begin{matrix}a={\frac {{{2.7810}^{-22}}{{n}_{Me}}{{Z}_{i}}{{Z}_{Me}}{{M}_{i}}}{e{\sqrt {Z_{i}^{{}^{2}/{}_{3}}+Z_{Me}^{{}^{2}/{}_{3}}}}\left({{M}_{i}}+{{M}_{Me}}\right)}},\\b=0.0793\cdot {{10}^{-36}}{\frac {{{Z}_{i}}^{\frac {1}{6}}\cdot {{\left({{Z}_{i}}~{{Z}_{Me}}\right)}^{\frac {1}{2}}}\cdot {{\left({{M}_{i}}+{{M}_{Me}}\right)}^{\frac {3}{2}}}}{{{\left({{Z}_{i}}^{\frac {2}{3}}+{{Z}_{Me}}^{\frac {2}{3}}\right)}^{\frac {3}{4}}}~{{M}_{i}}^{\frac {3}{2}}~{{M}_{Me}}^{\frac {1}{2}}}}{\frac {{n}_{Me}}{e}},\\\end{matrix}}

индекс i относится к иону, а индекс Me – к металлу,e – заряд электрона, Кл;  Z – зарядовое число элемента, M – массовое число элемента, n – концентрация частиц, м-3. Также необходимо учесть площадь поверхности (S), подверженной влиянию плазмы, и дисперсию дефектов относительно выбранной глубины с помощью нормального распределения с дисперсией Δy. Окончательно зависимость концентрации дефектов от глубины проникновения иона задается следующей зависимостью:

{{n}_{de{{f}_{x}}}}\left(x\right)={\frac {1}{S\Delta y{\sqrt {2\pi }}}}\int \limits _{0}^{\infty }{{{f}_{i}}\left(W\right){\frac {W}{2{{W}_{d}}}}{{\left({\frac {dy\left(W\right)}{dW}}\right)}^{-1}}}{{e}^{-{\frac {{(x-y)}^{2}}{2\Delta {{y}^{2}}}}}}dy.

Для учета эффектов, связанных с изменением во времени напряженности электрического поля в канале ускорителя, концентрации ионов плазмы, состава плазмы и других характеристик удобнее перейти к скоростям образования дефектов. Для этого указанная величина отнесена ко времени, за которое все частицы в приэлектродной области провзаимодействуют с поверхностью:

\Delta t={\frac {2\lambda }{\sqrt {\frac {2\ {{10}^{3}}RT}{\pi {{M}_{i}}}}}},

где λ – длина свободного пробега иона, м; R – универсальная газовая постоянная; T – температура ионной компоненты, К. Типичный график зависимости скорости роста относительного количества дефектов приведен на рисунке, полученный при использовании распределения Максвелла при температуре 30 000 К, напряжённости электрического поля 5 кВ/мм, в процессе проникновения в медные электроды ионов азота из плазмы с концентрацией 1024 м-3, что соответствует разряду в атмосфере.

Распределение скорости роста относительного количества дефектов материале по глубине проникновения ионов

Согласно [8] изменение удельного сопротивления материала задается следующим соотношением:

\Delta {{\rho }_{e}}={\frac {n_{Me}^{1/3}{{\left(3{{\pi }^{2}}\right)}^{1/3}}\hbar }{{e}^{2}}}{\frac {{n}_{def}}{{n}_{Me}}}{{\sigma }_{eff}},

где \hbar   — рационализованная постоянная планка, σeff – сечение рассеяния на дефекте кристаллической решетки.

Исходя из представленных данных, авторы заключают, что ионная бомбардировка поверхности металла приводит к накоплению от 5 до 90 % дефектов в кристаллической структуре меди за микросекунду воздействия на глубине до 3 нм. В свою очередь дефектообразование в структуре приводит к увеличению удельного сопротивления материала на величину до 90 % относительно исходного.

Литература
  1. Schneider M., Vincent G., Hogan J.D., Spray J.G. The use of a railgun facility for dynamic fracture of brittle materials. IEEE Transactions on Plasma Science, May 2015. DOI: https://doi.org/10.1109/TPS.2015.2396081
  2. Резников Б.И., Бобашев C.В., Жуков Б.Г., Куракин Р.О., Поняев С.А., Розов С.И. Эффективный коэффициент эрозии и ограничение скорости плазмы в канале электромагнитного рельсового ускорителя. ЖТФ, 2014, т. 84, № 4, с. 31–34. DOI: https://doi.org/10.21883/jtf.2016.11.43812.1819
  3. Watt T., Stefani F., Crawford M., Mark H., Parker J. Investigation of damage to solid-armature railguns at startup. IEEE Transactions on Plasma Science, 2007, vol. 43, pp. 214–218. DOI: https://doi.org/10.1109/TMAG.2006.887432
  4. Meger R.A., Cooper K., Jones H., Neri J., Qadri S., Singer I.L., Sprague J., Wahl K.J. Analysis of rail surfaces from a multishot railgun. IEEE Transactions on Magnetics, 2005, vol. 41, pp. 211–213. DOI: https://doi.org/10.1109/ELT.2004.1398051
  5. Любимов Г.А., Раховский В.И. Катодное пятно вакуумной дуги. УФН, 1978, т. 125, № 4, с. 665–706.
  6. Wang G., Tian S., Morris M., Morris S., Obradovic B., Balamurugan G., Tasch A. A computationally efficient ion implantation model: modified Kinchin-Pease model in Microelectronic Device Technology, Austin, TX, USA: SPIE-Int Soc Opt Eng, 1997, pp. 324–333.
  7. Nastasi M., Mayer J.W. Ion Implantation and Synthesis of Materials. Springer, Berlin, Heidelberg, 2006.
  8. Mott N., Jones H. The theory of the properties of metals and alloys. New York, Oxford University Press – Dover Publications, 1958.
Ваш браузер устарел и не обеспечивает полноценную и безопасную работу с сайтом.
Установите актуальную версию вашего браузера или одну из современных альтернатив.