CC BY
22
KINETIC SIMULATION OF VACUUM PLASMA EXPANSION BEYOND THE "PLASMA APPROXIMATION"
Vasily Y. Kozhevnikova, Andrey V. Kozyrevb, Aleksandr O. Kokovinc, Natalia S. Semeniukd
Institute of High Current Electronics, Laboratory of Theoretical Physics, Tomsk, Russian Federation
a e-mail: [email protected], corresponding author, ORCID iD: https://orcid.org/0000-0001-7499-0578
b e-mail: [email protected],
ORCID iD: https://orcid.org/0000-0002-7078-7991 c e-mail: [email protected],
ORCID iD: https://orcid.org/0000-0003-2068-7674 d e-mail: [email protected], ORCID iD: https://orcid.org/0000-0002-5972-2839
DOI: 10.5937/vojtehg70-37337; https://doi.org/ 10.5937/vojtehg70-37337
FIELD: Plasma physics
ARTICLE TYPE: Original scientific paper
Summary:
Introduction/purpose: One of the key approaches to solving an entire class of modern plasma physics problems is the so-called "plasma approximation". The most general definition of the "plasma approximation" is a theoretical approach to the electric field calculation of a system of charges under the electric quasi-neutrality condition. The purpose of this paper is to compare the results of the numerical simulation of the kinetic processes of the quasi-neutral plasma bunch expansion to the analytical solution of a similar kinetic model but in the "plasma approximation".
Methods: The given results are obtained by the methods of deterministic modeling based on the numerical solution of the system of Vlasov-Poisson equations.
ACKNOWLEDGMENT:
This paper is dedicated to the blessed memory of our good friend, permanent Chair of the TELFOR Telecommunications Forum Steering Committee Prof. Dr. Dorde Paunovic (19382021). During his long and rich career, he was a professor, head, mentor, colleague, friend, and, in every sense, much more than that.
The work was carried out within the framework of the state assignment of the Ministry of Science and Higher Education of the Russian Federation on the topics FWRM-2021-0007, FWRM-2021-0014.
Results: The provided comparison of the analytical expressions for the solution of kinetic equations in the "plasma approximation" and the numerical solutions of the Vlasov-Poisson equations system convincingly show the limitations of the "plasma approximation" in some important cases of the considered problem of plasma formation decay.
Conclusion: The theoretical results of this work are of great importance for understanding the shortcomings of the "plasma approximation", which can manifest themselves in practical applications of computational plasma physics.
Key words: physical kinetics, vacuum plasma, plasma expansion.
Introduction
The "plasma approximation" is known to be applicable to a low-frequency and steady-state phenomenon, popular among plasma scientists in various fields (Chen, 1984). Sometimes more appropriable term is used here "the plasma condition" (Nishikawa & Wakatani, 1990), i.e. the number of electrons in a Debye sphere is large enough to effect charge shielding. But utilizing it leads to inconsistencies in the equation of motion and prevents a proper, field-theoretic treatment of a condensed matter in the plasma state. This circumstance takes place due to the fact that if plasma reaches a quasi-neutral state ne« n, its space charge is approximately equal to zero p « 0. According to Poisson equation, this leads to E = 0. But the "plasma approximation" states that E * 0 and the electric field can be found elsewise (Chen, 1984). Such a separation of the initially consistent solution of plasma and field equations in the most cases leads to an ambiguous multivalued interpretation of the electric field definition. That seems to be the main methodological drawback of the "plasma approximation".
As only formal definitions equate with formal mathematics, so the "quasineutral" term correlates to "the plasma condition" ne« n. It refers to the profound tendency of plasma electrons to change their positions as a response to the electrostatic potential of ions to exponentially attenuate the Coulomb field, and is often taken as the definition of the "plasma approximation." We know that the the traditional term "neutral" already embraces the implications of quasineutrality, as no discrete medium remains neutral on characteristic scales sufficiently smaller to resolve isolated charges. In quasineutral media, the microscopic field fluctuates strongly, but on the particle scale it averages out as the differential volume element grows. The author of the monograph (Chen, 1984) claims that "the plasma approximation is almost the same as the condition of
CO
<u
o
O >
CM
of
UJ
a:
ZD O o
_J
<
o
X
o
LU
I—
>-
a: <
i—
< -j
CD >o
X LU I—
o
o >
0
quasineutrality discussed earlier but has a more exact meaning ... is a mathematical shortcut that one can use even for wave motions. it is usually possible to assume ne = ni and div E * 0 at the same time".
The physical kinetics of plasma provides a different point of view on a "quasineutrality" concept given from more fundamental positions. Indeed, if we consider simple two component plasma that consists of least of electrons and single-charged ions, then the ensemble of each type of particles in terms of physical kinetics is characterized by its distribution function, here fe and f\, respectively. Hence the number density of each particle type is a special case of the distribution function zero-moment
where (r, p) - phase-space coordinates, t - time variable.
Such zero-moments, like e.g. (1), of the particle distribution function, do not characterize the microscopic state of the ensemble of particles. They just represent particular macroscopic characteristics of certain plasma components. That is why the approximation of ne « n is the equality in a "weak form", i.e. the identity fe = f does not follow from the quasineutrality condition. Frequently, depending on a particular physical problem based, one has to introduce extra conditions for two or more additional (higher) moments of the distribution function in order to satisfy the "plasma approximation". If we assume that the "plasma approximation" is a convenient computational approach to a number of physical problems, then one has to determine the limits of its use.
This paper is aimed to clarify the details of the two-component vacuum plasma bunch expansion into free space by numerically solving the system of Poisson-Vlasov equations. Its main purpose is to show the features of this process without using the "plasma approximation". For simplicity, but without loss of generality, we solve the problem of plasma expansion in a one-dimensional Cartesian spatial configuration. The calculation results are compared with the exact self-similar solutions of the Vlasov equations with the "plasma approximation" (Dorozhkina & Semenov, 1998) pointing out possible shortcomings of the "plasma approximation".
Vacuum plasma expansion General terms
Let us consider a one-dimensional planar plasma bunch, consisting only of electrons and singly charged ions, located around the point x0, on
the x-axis. The bunch has a localization region of the order of Xc (spatial distribution half-width). We assume that both ions and electrons inside the bunch have Maxwellian velocity distributions with slightly different temperatures Te for electrons and T for ions. Assuming that the intial plasma is quasi-neutral, the particle distribution functions can be written in the following form
/ o \
fe,t ( X P )
N
SxyJ 8k
exp
m„
2m .kT
x
x exp
e, t
2
X
2
x J
c
(2)
where N0 - full number of particles in the bunch, me and mi - electron and ion rest mass, respectively, (x, p) - one-dimensional phase-space coordinates, S - bunch transversal cross section, Xc - characteristic spatial scale of a plasma bunch, and Te,i - initial ion and electron temperatures, respectively. If the initial electron distribution is assumed to be a non-Maxwellian, then it evolves to a Maxwellian one due to electron-electron elastic collisions. The Maxwellian distribution conserves, since collisions no longer have any influence on the electron distribution function, because the relevant term for elastic collisions is zero for a Maxwellian distribution. As the effects of collisions of the other kind are much more negligible with regard to electron-electrons, so the restriction to the Maxwellian initial distribution form is justified enough.
The electron and ion distribution function comply with the collisionless Boltzmann (Vlasov) equations without a magnetic field (Vlasov, 1968)
m
dt f
dt
dx
p^df± dx
qE
m
+ qE
f dp
ddf_ dp
0,
0,
(3)
where E - the electric field vector component along the x-axis, q - the electron charge.
During the expansion process, the assumption of a collisionless plasma is valid if the electron-electron collision time Tee >> r* = L/Cs, where L is the characteristic size and Cs = (qTe/m)1/2 is the ion sound velocity. As
Tee = 7e3/2/(5-10"6neA), where Coulomb log is A « 10 and A is the atomic number, we can estimate the condition for a number density ne << n* = Te2/(5- 10-12A12AL). For an antimony plasma (A = 51) at characteristic lengths of L = 1 cm and Te = 5 eV, the value of n* ~ 7-1010 cm-3 (Baitin & Kuzanyan, 1998).
In most cases where the domination of space-charge effects is significant and global plasma quasi-neutrality condition is not met, the system of equations (3) is complemented by the Poisson's equation in order to account the electric field in a self-consistent way
02( q ( x „ 0(
TT = —(ni - ne ) , E = -~T , (4)
ox SQ ox
where ( - electrostatic (electric) potential, S0 - vacuum dielectric permittivity, ne and ni are electron and ion number densities, respectively, that have to be found from (1).
"Plasma approximation"
Regarding the problem of interest, the first paper where the "plasma approximation" and the associated theoretical approach have been introduced is Gurevich's paper (Gurevich et al, 1966). The case of a halfinfinite plasma with a sharp boundary has been considered by using the Boltzmann electron distribution in order to find the electric field under the quasi-neutrality assumption
f ( +W
ne ( x
( x, t) = n0 exp
q(( x, t)
kT
(5)
where n0 - is the constant characteristic number density of ions. This assumption is not appropriate, since the electrostatic potential is non-stationary and the electron distribution function is different from a Boltzmann for a collisionless plasma. So, the total electron energy changes in time and the electron thermal energy becomes a source of the ion component acceleration in the expanding plasma. In a bounded collisionless plasma, the electrons are trapped and oscillate in a potential well, which is formed in order to satisfy a quasi-neutrality condition. Since the parameters of the potential well change during plasma expansion, so the variation of the electron energy becomes significant.
A more complicated approach to the implementation of the "plasma approximation" is based directly on the physical kinetics principles (Dorozhkina & Semenov, 1998). From the plasma quasi-neutrality
condition ne « n, the left and right parts of equations in (3) are multiplied by p and then integrated by the moment variable over the phase space. Considering the number densities definition and the boundary conditions for the distribution functions at |p| ^ <», the following equations can be obtained
j" PfedP + — J P2fedP + qEne = 0 o tJ m oxJ
a id (6)
— f pfidp +——J p2 fidP - qEn =
/ir * m ri-v *
In order to obtain the electric field strength from expression (6), the authors of the approach (Dorozhkina & Semenov, 1998) impose an additional approximation, namely, they assume that the plasma bunch is "currentless", so
q\ PfedP PfidP (7)
in this case by subtracting the second equation from the first one in (6), the following expression for the electric field is obtained
— f P2 fedP--- f P2 f,dP
m 5xJ m 5xJ (8) E = —e---.
2qne
We can now compare (8) to the electric field obtained from Gurevich's formula (5)
77 T d 1
E =--—In n. (9)
q ox
As it could be estimated, the given formulas (8) and (9) represent significantly different electric field values for the same plasma parameters. Namely, formula (9) gives a stronger electric field, which is an order of magnitude higher than the similar value obtained by formula (8). These arguments demonstrate the inconsistency of the "plasma approximation" concept. The whole point is the ambiguity of the definition of the electric field, and there are other approaches to obtaining the electric field in "plasma approximation" leading to different electric field estimations which are all different (Baitin & Kuzanyan, 1998).
CO
<u
o
O >
CM
of
UJ
a:
ZD O o
_J
<
o
X
o
LU
I—
>-
a: <
i—
< -j
CD >o
X LU I—
o
o >
0
The presented discussion highlights the main disadvantages of the so-called "plasma approximation". First, the variety of the electric field representations in the "plasma approximation" is determined by different theoretical approaches (liquid, kinetic or particle-in-cell) in use, obviously depriving the unambiguity of such approaches. Secondly, the calculation of the electric field requires some additional approximations that come far beyond the basic "quasi-neutrality" condition which is find to be insuffucient. In the scientific literature, the use of such approximations is given without sufficient justification (Dorozhkina & Semenov, 1998). Finally, there is no unambiguous way to identify physical situations where the electric field in the "plasma approximation" smoothly transforms into the field determined from the Poisson's equation or Maxwell's system of equations in the transition regions.
Numerical simulation
The obvious difficulties in choosing the correct formulation of the "plasma approximation" lead to the fact that the most accurate plasma dynamics has to be explained in terms of the complete Vlasov-Poisson system solution (3)-(4) where the electric field is determined in a self-consistent way. Here we choose a one-dimensional formulation of the problem in the Cartesian coordinates. Its advantages are obvious: the obtained solution results can be directly compared to the analytical formulas in the "plasma approximation" from (Dorozhkina & Semenov, 1998).
The direct numerical integration of (4) by using the trapezoidal or Simpson methods leads to significant inaccuracies associated with the accumulation of errors. To accurately determine the electric field and potential in this work, we used the advanced fourth order method (Knorr et al, 1980). As the computational phase space is restricted to the finite spatial interval xe[xmin, Xmax], so we apply E = 0 boundary conditions at the both sides of it.
In this paper, the system of partial differential equations (3) was solved numerically on a rectangular uniform phase-space grid (x, p) having 5000 per 2001 grid points for electrons and ions. The Vlasov equations have been solved by using the high-order Cheng-Knorr semi-Lagrangian method similar to that previously used (Zubarev et al, 2020, Kozhevnikov et al, 2021). The numerical solution algorithm was implemented in Mathworks MATLAB exploiting the embedded high-performance CPU capabilities. The results of numerical calculations have been validated by the computational grid and the computational time step value refining.
As an example, here we consider a two-component metallic plasma consisting of electrons and single-charged antimony ions Sb+. This plasma components are typical for vacuum discharge in diodes with antimony cathodes (Anders, 1997). For this plasma composition, a number of numerical calculation series have been carried out. The simulations were processed for a wide range of total number of particles No = 107 - 1013. In each calculation, it was assumed that the plasma bunch had a characteristic scale of Xc = 100 ^m, while the plasma was assumed to be nonequilibrium, i.e. Te = 1 eV, 7] = 5 eV that corresponds to real cathode plasma emission (Bugaev et al, 1975). For computational purposes, we restrict the spatial boundaries of the computational phase space to Xmin = - 2 cm, Xmax = 2 cm. It is sufficient to simulate the ionized state behavior far enough from the computational borders. The obtained results were compared with the analytical solutions of the Vlasov equations with an electric field in the "plasma approximation" (8) (Dorozhkina & Semenov, 1998).
Figure 1 shows the comparisons of the number densities distributions of electrons and ions at a time point of t = 500 ns for the cases of plasma expansion consisting of different initial number of particles. A quasi-neutral plasma distribution profile is shown with a black line in accordance with the analytical solution in the "plasma approximation" (Dorozhkina & Semenov, 1998). The results of the numerical calculations (without the "plasma approximation") - the number densities distribution profiles - are given for both electrons and ions to show the difference in their spatial distributions.
The first two plots in Figure 1 correspond to the decay of ionized states with a small number of particles - N0 = 107 and N0 = 109. Such initial distributions of charged particles cannot be called "plasma" due to the fact that the Debye length (Chen, 1984) is much larger than the characteristic scale of the bunch, i.e. Ad >> Xc. In the first case, the Debye length is much greater than xc, in the second case it has the same order of magnitude. In both cases, the numerical calculation shows a violation of the initially given quasi-neutrality conditions and a significant deviation from the "plasma approximation" profiles.
For a denser plasma (N0 = 1011), a different situation is observed. This case corresponds to true plasma decay Ad << Xc. The electron and ions number density distribution profiles in the Vlasov-Poisson model are close to the quasi-neutral profile obtained from the "plasma approximation". For the specified time point (t = 500 ns), quasi-neutrality is not significantly disturbed over the entire length of the plasma bunch. This situation is most fully described by the "plasma approximation": the plasma bunch expands
CO <1>
O
O >
CM
of
UJ
a:
ZD O O
_J
<
o
X
o
LU
I—
>-
Ql <
I—
< -J
CD >o
X LU I—
O
O >
O
with thermal velocities without significant loss of the initial quasi-neutrality. At the very beginning of the process, the more mobile and thermalized electronic component of the plasma is displaced with regard to the electrically neutral state, which leads to the appearance of a weak electric field close to the field in the "plasma approximation" in this case. Thus, plasma tends to remain quasi-neutral: if the ions move, then the electrons will follow them, and the electric field is adjusting to maintain the neutrality in accordance with the displacement of electrons and ions.
Finally, the most important case considered here represents the decay of a dense plasma bunch. The corresponding result of these calculations is given at the fourth plot in Figure 1. For a plasma bunch with the total number of particles equal to No = 1013, other regularities are observed. First of all, one can find a similarly to the previous case No = 1011. For No = 1013 the plasma bunch expands in time while preserving quasi-neutrality. But resulting from the calculations without the "plasma approximation", plasma decays faster. The profiles of a quasi-neutral plasma in the "plasma approximation" and that one obtained as a numerical solution of the Vlasov-Poisson equations (without the "plasma approximation") noticeably differ. In comparison with the decay of the ionized state (in No = 109), where the quasi-neutrality is noticeably violated, the decay of the dense plasma is more intense with regard to the reference calculation (black line in Figure 1).
Since the bunch quasi-neutrality is not violated (for N0 = 1013), it can be assumed that in this case other factors affect the plasma decay process, which are not taken into account by the "plasma approximation" (Dorozhkina & Semenov, 1998). As we have already mentioned previously, the kinetic formulation of the "plasma approximation" requires an additional "currentless" approximation (7). The "currentless" approximation is introduced independently of the quasi-neutrality approximation. It leads to a linearization of kinetic equations (3), so that the plasma number densities (black lines in Figure 1) for different number of particles have the same characteristic width and differ only by the scale factor of the curve magnitude. In reality, the system of Vlasov-Poisson equations is essentially non-linear and its solutions for various parameters do not scale. The calculations show that condition (7) is violated for a dense plasma. It leads to faster plasma expansion due to the nonlinearly increasing influence of electron and ion currents without affecting the total bunch quasi-neutrality.
Figure 1 - Comparative plasma number density distributions obtained from numerical calculations without the "plasma approximation" (red line and blue points) and from exact analytical formulas in the "plasma approximation" (black line)
Рисунок 1 - Сравнительные распределения концентрации компонентов плазмы, полученные из численных расчётов без "плазменного приближения" (цветные кривые) и по точным аналитическим формулам в "плазменном приближении" (черная линия)
Слика 1 - Упоредне расподеле концентраци^е компоненти плазме доби^ене из нумеричких прорачуна без „плазма апроксимаци]е"(обо]ене криве) и из тачних аналитичких формула у „плазма апроксимаци]и " (црна линц'а)
Conclusions
The results presented in this paper represent a counterexample convincingly showing the groundlessness of the "plasma approximation" for solving particular non-stationary plasma physics problems. The existing contradictions in the use of the "plasma approximation" can be formulated as the following theoretical statements:
• In all of the existing "plasma approximation" formulations, the electric field at t = 0 is nonzero, while the initial plasma is quasineutral, so the electric field initially is E = 0;
• Accurate numerical simulation (without the "plasma approximation") shows that for some cases of the quasi-neutral plasma bunch decay a local violation of plasma electrical neutrality appears. This leads to the electric field redistribution and affects further plasma dynamics, while in the "plasma approximation" plasma quasi-neutrality is preserved every time;
• Finally, the known "plasma approximations" are ambiguous. They require additional physical approximations depending on the theoretical approach in use. Such approximations in some cases make the "plasma approximation" less accurate.
Following the comparisons given in this paper, it can be argued that the "plasma approximation" can be used to study plasma decay without an external electric field only for a restricted range of dense plasma parameters. In other cases, the initial quasi-neutrality conditions do not guarantee the preservation of electrical neutrality during the whole ionized state decay process making the use of the "plasma approximation" unacceptable.
References
Anders, A. 1997. Ion charge state distributions of vacuum arc plasmas: The origin of species. Physical Review E, 55(1), pp.969-981. Available at: https://doi.org/10.1103/physreve.55.969.
Baitin, A.V. & Kuzanyan, K.M. 1998. A self-similar solution for expansion into a vacuum of a collisionless plasma bunch. Journal of Plasma Physics, 59(1), pp.83-90. Available at: https://doi.org/10.1017/s0022377897005916.
Bugaev, S.P., Litvinov, E.A., Mesyats, G.A. & Proskurovskii, D.I. 1975. Explosive emission of electrons. Soviet Physics Uspekhi, 18(1), pp.51-61. Available at: https://doi.org/10.1070/pu1975v018n01abeh004693.
Chen, F.F. 1984. Introduction to Plasma Physics and Controlled Fusion. New York, NY: Springer. Available at: https://doi.org/10.1007/978-1-4757-5595-4.
Dorozhkina, D.S. & Semenov, V.E. 1998. Exact Solution of Vlasov Equations for Quasineutral Expansion of Plasma Bunch into Vacuum. Physical Review Letters, 81(13), pp.2691-2694. Available at:
https://doi.org/10.1103/physrevlett.81.2691.
Gurevich, A.V., Pariiskaya, L.V. & Pitaevskii, L.P. 1966. Self-similar motion of rarefied plasma. Soviet Phisics JETP, 22(2), pp.449-454 [online]. Available at: http://jetp.ras.ru/cgi-bin/dn/e_022_02_0449.pdf [Accessed: 5 April 2022].
Kozhevnikov, V., Kozyrev, A., Kokovin, A. & Semeniuk, N. 2021. The Electrodynamic Mechanism of Collisionless Multicomponent Plasma Expansion in Vacuum Discharges: From Estimates to Kinetic Theory. Energies, 14(22), art.ID:7608. Available at: https://doi.org/10.3390/en14227608.
Knorr, G., Joyce, G. & Marcus, A. 1980. Fourth-order Poisson solver for the simulation of bounded plasmas. Journal of Computational Physics, 38(2), pp.227236. Available at: https://doi.org/10.1016/0021-9991(80)90054-6.
Nishikawa, K. & Wakatani, M. 1990. Basic Properties of Plasma. In: Plasma Physics, 8, pp.6-13. Springer, Berlin, Heidelberg: Springer Series on Atoms+Plasmas. Available at: https://doi.org/10.1007/978-3-662-02658-8_2.
Vlasov, A.A. 1968. The vibrational properties of an electron gas. Soviet Physics Uspekhi, 10(6), pp.721-733. Available at: https://doi.org/10.1070/pu1968v010n06abeh003709.
Zubarev, N.M., Kozhevnikov, V.Y., Kozyrev, A.V., Mesyats G.A., Semeniuk, N.S., Sharypov, K.A., Shunailov, S.A. & Yalandin, M.I. 2020. Mechanism and dynamics of picosecond radial breakdown of a gas-filled coaxial line. Plasma Sources Science and Technology, 29(12), art.ID:125008. Available at: https://doi.org/10.1088/1361-6595/abc414.
КИНЕТИЧЕСКОЕ МОДЕЛИРОВАНИЕ РАСШИРЕНИЯ ВАКУУМНОЙ ПЛАЗМЫ ЗА ПРЕДЕЛАМИ "ПЛАЗМЕННОГО ПРИБЛИЖЕНИЯ"
Василий Ю. Кожевников, корреспондент, Андрей В. Козырев, Александр О. Коковин, Наталия С. Семенюк
Институт сильноточной электроники, Лаборатория теоретической физики, г. Томск, Российская Федерация
РУБРИКА ГРНТИ: 29.27.03 Общие свойства плазмы
29.27.47 Численные методы в физике плазмы ВИД СТАТЬИ: оригинальная научная статья
Резюме:
Введение/цель: Одним из ключевых подходов к решению целого класса задач современной физики плазмы является так называемое "плазменное приближение". Наиболее общее определение "плазменного приближения" - это теоретический подход к вычислению электрического поля системы зарядов в
условиях их электрической квазинейтральности. Целью данной работы является сравнение результатов численного моделирования кинетических процессов распада сгустка квазинейтральной плазмы с аналитическим самоподобным решением аналогичной кинетической модели, полученным в условиях "плазменного приближения".
Методы: Приведенные результаты получения детерминистического моделирования, основанные на вычисленном обнаружении системы обнаружения Власова-Пуассона.
Результаты: Сравнение аналитических выражений решения кинетических уравнений в "плазменном приближении" и численных решений системы уравнений Власова-Пуассона убедительно показывают ограниченность использования "плазменного приближения" в ряде случаев рассматриваемой задачи о распаде плазменного образования.
Выводы: Теоретические результаты данной работы имеют большое значение для понимания недостатков "плазменного приближения", которые могут проявляться в практических приложениях вычислительной физики плазмы.
Ключевые слова: физическая кинетика, вакуумная плазма, разлет плазмы.
КИНЕТИЧКА СИМУЛАЦША ЕКСПАНЗШЕ ВАКУУМСКЕ ПЛАЗМЕ ИЗВАН ГРАНИЦА „ПЛАЗМА АПРОКСИМАЦШЕ"
Басили] ^ Кожевников, аутор за преписку, Андре] В. Козирев, Александар О. Коковин, Наталща С. Семен|ук
Институт за високостру]ну електронику, Лаборатор^а за теорийку физику, Томск, Руска Федерац^а
ОБЛАСТ: физика (општа сво]ства плазме) ВРСТА ЧЛАНКА: оригинални научни рад
Сажетак:
Увод/цил: Jедан од клучних приступа решаваъу читаве класе проблема модерне физике плазме }есте такозвана апроксимаци}а плазме. Нарпшт^е се дефинише као теорирки приступ израчунаваъу електричног пола наелектрисаних честица под условом електричне квази-неутралности. Цил> овог рада рсте да упореди резултате нумеричке симулацир кинетичких процеса ширена квазинеутралне згуснуте плазме са аналитичким реше^ем сличног кинетичког модела али у апроксимацир плазме.
Методе: Резултати су доби}ени методама детерминистичког моделована заснованим на нумеричком решету Власов-Поасоновог система ]едначина.
Резултати: Представлено поре^ене аналитичких израза за решаване кинетичких ]едначина у апроксимаци]и плазме, као и нумеричка решена Власов-Поасоновог система ]едначина, убедл>иво показу}у ограничена апроксимаци]е плазме у неким важним случа}евима разматраног проблема распадана формаци]а плазме.
Закъучци: Теорирки резултати овог рада од великог су знача}а за разумеване недостатака апроксимаци}е плазме до щих може доЬи приликом практичне примене комп}утерске физике плазме.
Къучне речи: физичка кинетика, вакуум плазма, експанзи}а плазме.
Paper received on / Дата получения работы / Датум приема чланка: 08.04.2022. Manuscript corrections submitted on / Дата получения исправленной версии работы / Датум достав^а^а исправки рукописа: 22.06.2022.
Paper accepted for publishing on / Дата окончательного согласования работы / Датум коначног прихвата^а чланка за об]ав^ива^е: 24.06.2022.
© 2022 The Authors. Published by Vojnotehnicki glasnik / Military Technical Courier (www.vtg.mod.gov.rs, втг.мо.упр.срб). This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/rs/).
© 2022 Авторы. Опубликовано в «Военно-технический вестник / Vojnotehnicki glasnik / Military Technical Courier» (www.vtg.mod.gov.rs, втг.мо.упр.срб). Данная статья в открытом доступе и распространяется в соответствии с лицензией «Creative Commons» (http://creativecommons.org/licenses/by/3.0/rs/).
© 2022 Аутори. Об]авио Во^отехнички гласник / Vojnotehnicki glasnik / Military Technical Courier (www.vtg.mod.gov.rs, втг.мо.упр.срб). Ово jе чланак отвореног приступа и дистрибуира се у складу са Creative Commons лиценцом (http://creativecommons.org/licenses/by/3.0/rs/).
