Plasma and magnetic field dynamics in POS: pic simulations

2D3V particle in cell computer simulations for plasma, magnetic field and current density spatio-temporal dynamics in plasma opening switch (POS) are presented for various initial plasma densities. The current channel formation in the POS’s plasma, its passage through the plasma bridge, and its dest...

Повний опис

Збережено в:
Бібліографічні деталі
Дата:2022
Автори: Manuilenko, O.V., Onishchenko, I.N., Pashchenko, A.V., Pashchenko, I.A., Yuferov, V.B.
Формат: Стаття
Мова:English
Опубліковано: Національний науковий центр «Харківський фізико-технічний інститут» НАН України 2022
Назва видання:Problems of Atomic Science and Technology
Теми:
Онлайн доступ:http://dspace.nbuv.gov.ua/handle/123456789/195886
Теги: Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:Plasma and magnetic field dynamics in POS: pic simulations / O.V. Manuilenko, I.N. Onishchenko, A.V. Pashchenko, I.A. Pashchenko, V.B. Yuferov // Problems of Atomic Science and Technology. — 2022. — № 6. — С. 55-59. — Бібліогр.: 12 назв. — англ.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-195886
record_format dspace
spelling irk-123456789-1958862023-12-08T12:31:28Z Plasma and magnetic field dynamics in POS: pic simulations Manuilenko, O.V. Onishchenko, I.N. Pashchenko, A.V. Pashchenko, I.A. Yuferov, V.B. Plasma dynamics and plasma-wall interaction 2D3V particle in cell computer simulations for plasma, magnetic field and current density spatio-temporal dynamics in plasma opening switch (POS) are presented for various initial plasma densities. The current channel formation in the POS’s plasma, its passage through the plasma bridge, and its destruction at the boundary of the plasma jumper are studied. It is shown, that during the penetration of the magnetic field into the POS’s plasma and current channel formation, electron vortices are excited. The passage speed of the POS’s plasma by vortices decreases with increasing plasma density, which is qualitatively consistent with estimates based on classical electron magnetohydrodynamics (EMHD). Наведено результати 2D3V комп’ютерного моделювання методом макрочастинок просторово-часової динаміки плазми, магнітного поля та густини струму в плазмовому комутаторі струму (ПКС) для різних початкових густин плазми. Досліджено формування струмового каналу в плазмі ПКС, його проходження через плазму та руйнування на межі плазмової перемички. Показано, що при проникненні магнітного поля у плазму ПКС та формуванні струмового каналу збуджуються електронні вихори. Швидкість їх проходження через плазму ПКС знижується зі зростанням густини плазми, що якісно узгоджується з оцінками на основі класичної електронної магнітогідродинаміки (ЕМГД). 2022 Article Plasma and magnetic field dynamics in POS: pic simulations / O.V. Manuilenko, I.N. Onishchenko, A.V. Pashchenko, I.A. Pashchenko, V.B. Yuferov // Problems of Atomic Science and Technology. — 2022. — № 6. — С. 55-59. — Бібліогр.: 12 назв. — англ. 1562-6016 PACS: 52.75.-d, 52.75.Kq, 52.40.Hf, 52.65,+z, 52.70.Kz, 94.20 DOI: https://doi.org/10.46813/2022-142-055 http://dspace.nbuv.gov.ua/handle/123456789/195886 en Problems of Atomic Science and Technology Національний науковий центр «Харківський фізико-технічний інститут» НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language English
topic Plasma dynamics and plasma-wall interaction
Plasma dynamics and plasma-wall interaction
spellingShingle Plasma dynamics and plasma-wall interaction
Plasma dynamics and plasma-wall interaction
Manuilenko, O.V.
Onishchenko, I.N.
Pashchenko, A.V.
Pashchenko, I.A.
Yuferov, V.B.
Plasma and magnetic field dynamics in POS: pic simulations
Problems of Atomic Science and Technology
description 2D3V particle in cell computer simulations for plasma, magnetic field and current density spatio-temporal dynamics in plasma opening switch (POS) are presented for various initial plasma densities. The current channel formation in the POS’s plasma, its passage through the plasma bridge, and its destruction at the boundary of the plasma jumper are studied. It is shown, that during the penetration of the magnetic field into the POS’s plasma and current channel formation, electron vortices are excited. The passage speed of the POS’s plasma by vortices decreases with increasing plasma density, which is qualitatively consistent with estimates based on classical electron magnetohydrodynamics (EMHD).
format Article
author Manuilenko, O.V.
Onishchenko, I.N.
Pashchenko, A.V.
Pashchenko, I.A.
Yuferov, V.B.
author_facet Manuilenko, O.V.
Onishchenko, I.N.
Pashchenko, A.V.
Pashchenko, I.A.
Yuferov, V.B.
author_sort Manuilenko, O.V.
title Plasma and magnetic field dynamics in POS: pic simulations
title_short Plasma and magnetic field dynamics in POS: pic simulations
title_full Plasma and magnetic field dynamics in POS: pic simulations
title_fullStr Plasma and magnetic field dynamics in POS: pic simulations
title_full_unstemmed Plasma and magnetic field dynamics in POS: pic simulations
title_sort plasma and magnetic field dynamics in pos: pic simulations
publisher Національний науковий центр «Харківський фізико-технічний інститут» НАН України
publishDate 2022
topic_facet Plasma dynamics and plasma-wall interaction
url http://dspace.nbuv.gov.ua/handle/123456789/195886
citation_txt Plasma and magnetic field dynamics in POS: pic simulations / O.V. Manuilenko, I.N. Onishchenko, A.V. Pashchenko, I.A. Pashchenko, V.B. Yuferov // Problems of Atomic Science and Technology. — 2022. — № 6. — С. 55-59. — Бібліогр.: 12 назв. — англ.
series Problems of Atomic Science and Technology
work_keys_str_mv AT manuilenkoov plasmaandmagneticfielddynamicsinpospicsimulations
AT onishchenkoin plasmaandmagneticfielddynamicsinpospicsimulations
AT pashchenkoav plasmaandmagneticfielddynamicsinpospicsimulations
AT pashchenkoia plasmaandmagneticfielddynamicsinpospicsimulations
AT yuferovvb plasmaandmagneticfielddynamicsinpospicsimulations
first_indexed 2025-07-17T00:09:13Z
last_indexed 2025-07-17T00:09:13Z
_version_ 1837850642130927616
fulltext ISSN 1562-6016. Problems of Atomic Science and Technology. 2022. №6(142). Series: Plasma Physics (28), p. 55-59. 55 PLASMA DYNAMICS AND PLASMA-WALL INTERACTION https://doi.org/10.46813/2022-142-055 PLASMA AND MAGNETIC FIELD DYNAMICS IN POS: PIC SIMULATIONS O.V. Manuilenko1,2*, I.N. Onishchenko1, A.V. Pashchenko1, I.A. Pashchenko1, V.B. Yuferov1 1National Science Center ”Kharkov Institute of Physics and Technology”, Kharkiv, Ukraine 2V.N. Karazin Kharkiv National University, Kharkiv, Ukraine *E-mail: ovm@kipt.kharkov.ua 2D3V particle in cell computer simulations for plasma, magnetic field and current density spatio-temporal dynamics in plasma opening switch (POS) are presented for various initial plasma densities. The current channel formation in the POS’s plasma, its passage through the plasma bridge, and its destruction at the boundary of the plasma jumper are studied. It is shown, that during the penetration of the magnetic field into the POS's plasma and current channel formation, electron vortices are excited. The passage speed of the POS's plasma by vortices decreases with increasing plasma density, which is qualitatively consistent with estimates based on classical electron magnetohydrodynamics (EMHD). PACS: 52.75.-d, 52.75.Kq, 52.40.Hf, 52.65,+z, 52.70.Kz, 94.20 INTRODUCTION A POS is a plasma bridge between two electrodes, usually in vacuum coaxial line. POSes are widely used for voltage multiplication in pulsed-power science and technology including high-current pulsed electron accelerators [1-4]. The high-voltage pulse generator with POS is an electric RLC circuit closed through the plasma bridge. One end of the POS is connected to capacitive energy storage, the other ‒ to a load, usually high current vacuum diode. When energy storage is turned on, the increasing in time current is closed through the plasma bridge, the electrical energy of the capacitive energy storage CUo 2/2 is converted into magnetic energy LI2/2. After some time, the plasma bridge impedance increases sharply, and the current is quickly decreased, an overvoltage pulse appears UPOS ~ LdI/dt > Uo, the accumulated energy is transferred to the load. The POS opening time is much shorter than the current rise time. This shortening causes a compression of the energy and as a result ‒ an increase in the load power. The main element of the POS is a plasma bridge with a current. Its spatio-temporal dynamics determines the energy accumulation time (current rise time) and the POS opening time (current interruption time in the plasma bridge), which, in turn, determines the magnitude of the voltage multiplication UPOS/Uo. The plasma bridge spatio-temporal dynamics is determined by the initial conditions and the current flow through the plasma. The only way to control the processes in the plasma bridge is the proper choice of the initial conditions - the plasma density and composition, its initial spatial distribution, the dimensions of the cathode, anode, and the length of the plasma bridge. Despite many years of using POS, their experimental and theoretical studies [1, 2, 5-7], presently there is no complete theory, let alone computational models for predicting the POS parameters. Usually, the POS plasma density ne lies in the range of 1012...1014 cm-3, the cathode radii Rc are 1...5 cm, the anode radii R are 4…10 cm, i.e., the characteristic distances between the cathode and the anode a are 3...5 cm. The composition of the plasma can be different, but in most cases it is a plasma that contains carbon ions, charge Z = 1...3. In this parameters range mec2/e2 < nea2 < Mic2/Ze2, where me ‒ electron mass; Mi ‒ ion mass; c is the speed of light, i.e. the plasma and electromagnetic field dynamics in a POS is described by EMHD [8], i.e., the current is carried mainly by electrons, and the penetration of the electromagnetic field into the plasma is determined by the rapid penetration of the magnetic field B into the plasma in the form of a quasi-magnetostatic wave [8] with the speed uf = cBr2/8e/r(1/ner2) [5]. Computer simulations for magnetic field penetration into POS plasma and current loop formation in it in EMHD approximation for various initial plasma densities, currents, and POS geometries [9] shown that the current loop dynamics in the POS is determined by the fast magnetic field penetration in plasma due to the Hall effect. The strong dependence of the current loop longitudinal velocity on the transverse coordinate lead to the formation of the narrow S-shaped current loop even in homogeneous plasma. It was shown that the control parameters influencing the dynamics of the magnetic field and the current loop motion in the POS are the initial plasma density, current and POS geometry. Below the results of 2d3v PIC simulations of plasma dynamics in POS for different plasma densities ne = 1012...1014 cm-3 are presented. The fully relativistic electromagnetic particle-in-cell code XOOPIC [10] was used for simulations. SIMULATION RESULTS AND DISCUSSION Fig. 1 shows the POS geometry. The simulation domain is the space between cathode Rc and anode R, as well as between z = 0 and z = ZL, ZL = 11 cm. Initially electrons and plasma ions occupy the region {Rc, R}{Zpin, Zpf}, Zpin = 1 cm, Zpf = 10 cm. Anode is a mailto:ovm@kipt.kharkov.ua 56 ISSN 1562-6016. Problems of Atomic Science and Technology. 2022. №6(142) metal boundary, where Ez = 0, electrons and ions that came to the anode are absorbed. Cathode is also a metal boundary, where Ez = 0, the particles that came to the cathode are absorbed. On the segment {Zpin, Zpf}, where the plasma touches the cathode, if the self-consistent electric field Er > 2105 V/cm, a self-consistent explosive emission of electrons from the cathode occurs. On the left border, an electromagnetic field corresponding to the current I, which is supplied by the energy source and passes through the cathode - anode gap, is set. Impedance boundary conditions for free space are set on the right boundary. Particles that have reached the right boundary are removed. Fig. 1. POS is a coaxial line segment filled with plasma. z, r ‒ longitudinal and transverse coordinate For spatially homogeneous plasma we have for penetration velocity of the electromagnetic field into the plasma uf = 9.941017 I/ner2, where uf is in cm/s, r ‒ in cm, I ‒ in A, ne ‒ in cm-3. Estimates show that at ne = 1013 cm-3, I = 100 kA, the quasi ‒ magnetostatic wave passage time throw POS plasma bridge with a length of 10 cm is from 1 ns, for cathode radius Rc = 1 cm, to 6.25 ns, for Rc = 2.5 cm. At such time intervals, the ions can be considered immobile, and their influence on the processes of current transfer through the plasma can be neglected. Therefore, in the simulations below, the ions were immobile. Test simulations with mobile ions showed that their movement practically does not affect the process of formation and passage of the current channel through the plasma bridge. Figs. 2-6 show space-time dynamics of plasma, current channel and magnetic field for plasma densities ne = {1012, 1013, 1014} cm-3 at I = 100 kA. Initially plasma density is uniformly distributed in space. Rc = 1 cm, R = 7.5 cm, the plasma bridge length Lp = 9 cm. Figs. 2, 3 show the projection of the phase space {r, z, ; r, z; t} onto the configuration space {r, z} for plasma electrons (see Fig. 2) and electrons emitted from the cathode (see Fig. 2) as a result of explosive emission at different times. Figs. 4, 5 show the corresponding to Figs. 2, 3 electron current density = {jr(r, z, t), jz(r, z, t)} (current loop) at same time moments, and Fig. 6 shows corresponding magnetic field B(r,z,t). Current density and magnetic field are depicted as lines of the same level, which is displayed in color. This makes it possible to clearly compare the dynamics of charged particles with the current density and the field. Fig. 2. Plasma electrons in the configuration space {r, z}. Lines correspond to the times {0.11, 0.33, 0.55, 2.2, 5.5} ns. Columns correspond to initial density ne0 = 1012 cm-3 (left); ne0 = 1013 cm-3 (middle); ne0 = 1014 cm-3(right) Fig. 3. Electrons emitted from cathode in the configuration space {r, z}. Lines correspond to the times {0.11, 0.33, 0.55, 2.2, 5.5} ns. Columns correspond to initial density ne0 = 1012 cm-3 (left); ne0 = 1013 cm-3 (middle); ne0 = 1014 cm-3(right) Analysis Figs. 2-6 shows that initially a narrow current radial channel is formed on the left side of the ISSN 1562-6016. Problems of Atomic Science and Technology. 2022. №6(142). 57 plasma bridge, closer to the energy source. The current is carried by plasma electrons and electrons generated as a result of emission from the cathode. Current channel electrons are accelerated by the radial electric field Er, when passing from the cathode to the anode, they are displaced to the right by the force of r  B. Fig. 4. Spatio-temporal dynamics of electron current density, jr component, in the POS's plasma. Lines correspond to the times {0.11, 0.33, 0.55, 5.5} ns. Columns correspond to initial density ne0 = 1012 cm-3 (left); ne0 = 1013 cm-3 (middle); ne0 = 1014 cm-3 (right) Fig. 5. Spatio-temporal dynamics of electron current density, jz component, in the POS's plasma. Lines correspond to the times {0.11, 0.33, 0.55, 5.5} ns. Columns correspond to initial density ne0 = 1012 cm-3 (left); ne0 = 1013 cm-3 (middle); ne0 = 1014 cm-3(right) At the beginning of the current channel formation process, the electron trajectories are radial, with a small displacement in the axial direction. This leads to the formation in the current channel of a narrow layer of positive charge on the side of the energy source and a narrow layer of negative charge on the side of the load, which, in turn, leads to the appearance of a longitudinal electric field Ez in the current channel and to the radial drift of electrons in crossed Ez  B fields. Fig. 6. Spatio-temporal dynamics of magnetic field penetration into the POS's plasma. B (r,z,t) component is presented. Lines correspond to the times {0.33, 0.55, 5.5} ns. Columns correspond to initial density ne0 = 1012 cm-3 (left); ne0 = 1013 cm-3 (middle); ne0 = 1014 cm-3(right) At the cathode and anode Ez = 0. Near the cathode, in the plasma, the radial component of the electric field Er is increased, compared to the vacuum field. This is due to the fact that, firstly, as a result of a significant difference in masses, even at the same temperatures, electrons quickly go to the cathode, forming a positively charged layer near the cathode along the entire length, where the cathode touches the plasma, and, secondly, in the current channel, a narrow layer of positive charge is formed on the side of the energy source as a result of the movement of current-carrying electrons, which at the same time are shifted to the side of the load. This positive volume charge enhances Er, which leads to the explosive emission of electrons from the cathode and their acceleration into the current channel (see Fig. 3). Near the anode, in the plasma, Er is suppressed by the positive volume charge. Electron emission from the cathode “ties” the left border of the current channel near the cathode to the left border of the plasma bridge (see Fig. 3). Over time, the length of the cathode region emitting electrons increases, its right border shifts in the direction of the load (see Fig. 3). The current channel near the anode moves in the longitudinal direction due to the drift of electrons in crossed Er  B fields (see Figs. 4, 5). When the current channel reaches the right border of the plasma bridge, it consists mainly of the electrons emitted by the cathode (see Figs. 2, 3). If this flow of electrons falls into the mode of magnetic self-isolation, then the current between the anode and the cathode is completely interrupted, the energy of the magnetic field is transferred to the load. In most cases, there is not a complete breakdown of the current between the cathode and the anode, part of the electrons of the current channel reaches the anode. With the help of an external magnetic field, it is possible to cut off this part of the 58 ISSN 1562-6016. Problems of Atomic Science and Technology. 2022. №6(142) electron current from the anode and achieve a complete interruption of the current. As can be seen from Figs. 2-6, during the penetration of the magnetic field into the plasma at its left boundary, electron vortices are excited, in which electrons rotate against the hour hand. In the center, they have a positive charge - the density of the electronic component of the plasma is lower than in undisturbed plasma. The speed of vortices movement in the load direction is (0.2....5)uf. As the plasma density increases, the time of the magnetic field penetration into the plasma increases, which is qualitatively consistent with estimates based on classical EMHD – uf  1/ne. The characteristic size of vortices decreases with increasing plasma density. As can be seen from Figs. 4, 5, the current channel passes along the boundary of vortices – undisturbed plasma. From Fig. 6, it can be seen that the magnetic field is penetrated into the POS plasma, filling the volume of the vortices. Electrons captured by the vortex are excluded from the current transfer between the cathode and the anode, which can lead to an increase in the resistance of the POS at the energy accumulation stage, if this is not compensated by electrons emitted from the cathode due to explosive emission. In addition, after POS opening, the electrons captured by the vortices can quickly short- circuit the cathode – anode gap of the POS. Such vortices can be described [11, 12] by the complete Maxwell's equations and by the relativistic plasma hydrodynamics equations, from which the equations for the electronic fluid vorticity can be obtained: eee vvdivv t   )()()( +−=+   , where )/( cAevmrot eee  −=  ; ArotB  = ; e ‒ relati- vistic factor; ev  ‒ speed of the plasma electronic component. As can be seen from (1),   is transferred along the electron flow, and the terms on the right-hand side of (1) play the role of a source or sink of vorticity. Eq. can be rewritten as ee vndtd  )()/(/ = , where )(// +=  evtdtd . If 0)( = ev  , then enI /=  does not change along the electron flow. This is the case for flat geometry, if the electrons move in the {x, y} plane, and   has only z-component. Only in this case, the structure of the electron vortex was investigated and it was shown that the electron density in its center is significantly lower than the ion density, i.e., there is a violation of quasi-neutrality, the magnetic field in the center of the vortex is strengthened, compared to the peripheral magnetic field [11]. We see a similar structure of vortices in numerical simulations. Analysis Figs. 2-6 shows that four stages of POS operation can be distinguished in our case: (1) the initial stage, when a current channel is formed near the left border of the plasma bridge, (2) drift movement of a current channel part located near the anode in the load direction with simultaneous penetration of the magnetic field (and the current channel, since ( ) Hrottrj  ~, ) into the POS plasma in the form of quasi-magnetostatic wave and vortices on its front, and by pressing the part of the current channel located near the cathode to the cathode, (3) the exit of the current channel to the right border of the plasma jumper, (4) exit of the current channel behind the plasma bridge and complete or partial magnetic self- isolation of electrons, which leads to a complete or partial interruption of the current in the cathode-anode gap. In the third and fourth stages, the current channel consists mainly of electrons emitted from the cathode. CONCLUSIONS Computer simulations of spatio-temporal dynamics of plasma and electromagnetic fields in the POS for a wide range of plasma densities is presented. Numerical simulations were performed using the particle-in-cell method with solution complete Maxwell's equations and relativistic equations for charged particles motion. It is shown that theoretical models based on EMHD describe quite well the movement of the current channel through the POS plasma. This allows someone to use a simple expression for the velosity of magnetic field penetration into the plasma in the form of a quasi-magnetostatic wave for estimating the energy storage time in inductive energy storage devices with POS. The general picture of the current channel dynamic in the POS plasma, which follows from the numerical simulations, is as follows. Initially a narrow radial current channel is formed on the left side of the plasma bridge, closer to the energy source. Electrons in the current channel are accelerated by the radial electric field. When they passing from the cathode to the anode, they are shifted in the load direction by the force r  B.. This leads to the formation in the current channel of a narrow positively charged layer on the energy source side and negatively charged layer on the load side, which leads to the appearance of a longitudinal electric field in the current channel and a radial drift of electrons in crossed fields Ez  B. Near the cathode, in the plasma, the radial component of the electric field is increased, compared to the coaxial line vacuum field. This is due to the forming a positively charged layer near the cathode (sheath). This positive volume charge strengthens the radial component of the electric field, which leads to the explosive emission of electrons from the cathode. Over time, the length of the cathode region emitting electrons increases, its right border shifts in the load direction, and the current channel is pressed against the cathode (magnetic self- isolation). The current channel near the anode moves in the longitudinal direction due to the electrons drift in crossed Er  B fields. In the process of penetration of the magnetic field into the plasma (and passing by the current channel through the plasma jumper), electron vortices are excited. The vortices motion velosity in the load direction is (0.2...0.5)uf. As the plasma density increases, the time of passage of the current channel through the plasma bridge increases, which is consistent with estimates based on classical electron magnetohydrodynamics: uf  1/ne. When the current channel reaches the right border of the plasma bridge, it consists mainly of electrons emitted by the cathode. If this flow of electrons falls into the mode of magnetic ISSN 1562-6016. Problems of Atomic Science and Technology. 2022. №6(142). 59 self-isolation, then the current between the anode and the cathode is interrupted. In most cases, there is an incomplete current breaking. Part of the current channel electrons reaches the anode. Four stages of POS operation can be distinguished: the initial stage, when a current channel is formed, penetration of the magnetic field and the current channel into the POS plasma in the form of quasi-magnetostatic wave with vortices on its front, the exit of the current channel to the border of the plasma bridge, and exit of the current channel behind the plasma jumper, complete or partial magnetic self-isolation of electrons, which leads to interruption of the current in the cathode-anode gap. REFERENCES 1. C.W. Mendel, S.A. Goldshtein. A fast opening switch for use in REB diode experiments // Journal of Applied Physics. 1977, v. 48, p. 1004-1007. 2. R.A. Meger, R.J. Commisso, G. Cooperstein, S.A. Goldstein. Vacuum inductive store/pulse compression experiments on a high power accelerator using plasma opening switches // Applied Physics Letters. 1983, v. 42, p. 943-949. 3. E.I. Skibenko, V.B. Yuferov. Small-size direct-action electron accelerator with a high-efficiency nanosecond plasma-current switch // Problems of Atomic Science and Technology. Series “Plasma Electronics and New Methods of Acceleration” (122). 2019, № 4, p. 10-14. 4. V.B. Yuferov, E.I. Skibenko, I.N. Onishchenko, V.G. Artyuch, O.S. Druy. Investigation of high-current plasma opening switch at low gas pressure // Problems of Atomic Science and Technology. Series “Nuclear Physics Investigations”. 2000, № 2, p. 100-102. 5. A.V. Gordeev, A.S. Kingsep, L.I. Rudakov. Electron magnetohydrodynamics // Physics Reports. 1994, v. 243, p. 215-315. 6. G.I. Dolgachev, E.D. Kazakov, Yu.G. Kalinin, et al. RS-20MR high-current relativistic electron beam generator based on a plasma opening switch and its applications // Plasma Physics Reports. 2019, v. 45, p. 315-324. 7. S.V. Loginov. Experimental investigation of voltage scaling in a plasma switch // Russian Physics Journal. 2020, v. 62, p. 1976-1981. 8. A.S. Kingsepp, Yu.V. Mokhov, K.V. Chukbar. Non- linear skin effects in a plasma // Soviet Journal of Plasma Physics. 1984, v. 10, p. 495-499. 9. O.V. Manuilenko, I.N. Onishchenko, A.V. Pash- chenko, I.A. Pashchenko, V.A. Soshenko, V.G. Svi- chensky, V.B. Yuferov, B.V. Zajtsev. Magnetic field dynamics in plasma opening switch // Problems of Atomic Science and Technology. Series “Nuclear Physics Investigations” (136). 2021, № 6, p. 61-66. 10. J.P. Verboncoeur, A.B. Langdon, N.T. Gladd. An object-oriented electromagnetic PIC code // Computer Physics Communications. 1995, v. 87, № 1, 2, p. 199-211. 11. A.V. Gordeev, S.V. Levchenko. Is the Abrikosov model applicable for describing electronic vortices in plasmas? // JETP Letters. 1998, v. 67, № 7, p. 482-488. 12. A.V. Gordeev. «Mixing» of Lagrange invariants and a new class of integrable equations for the Hall quasi- equilibria // Plasma Physics Reports. 1997, v. 23, № 2, p. 92-101. Article received 13.10.2022 ДИНАМІКА ПЛАЗМИ І МАГНІТНОГО ПОЛЯ У ПКC: МОДЕЛЮВАННЯ МЕТОДОМ МАКРОЧАСТИНОК О.В. Мануйленко, І.М. Оніщенко, А.В. Пащенко, І.А. Пащенко, В.Б. Юферов Наведено результати 2D3V комп'ютерного моделювання методом макрочастинок просторово-часової динаміки плазми, магнітного поля та густини струму в плазмовому комутаторі струму (ПКС) для різних початкових густин плазми. Досліджено формування струмового каналу в плазмі ПКС, його проходження через плазму та руйнування на межі плазмової перемички. Показано, що при проникненні магнітного поля у плазму ПКС та формуванні струмового каналу збуджуються електронні вихори. Швидкість їх проходження через плазму ПКС знижується зі зростанням густини плазми, що якісно узгоджується з оцінками на основі класичної електронної магнітогідродинаміки (ЕМГД).