Precision of a FDTD method to simulate cold magnetized plasmas

The finite difference time domain (FDTD) method is applied to describe the propagation of the transverse electromagnetic waves through the magnetized plasmas. The numerical dispersion relation is obtained in a cold plasma approximation. The accuracy of the numerical dispersion is calculated as a fun...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:Вопросы атомной науки и техники
Datum:2014
Hauptverfasser: Pavlenko, I.V., Melnyk, D.A., Prokaieva, A.O., Girka, I.O.
Format: Artikel
Sprache:Englisch
Veröffentlicht: Національний науковий центр «Харківський фізико-технічний інститут» НАН України 2014
Schlagworte:
Online Zugang:https://nasplib.isofts.kiev.ua/handle/123456789/81191
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Zitieren:Precision of a FDTD method to simulate cold magnetized plasmas / I.V. Pavlenko, D.A. Melnyk, A.O. Prokaieva, I.O. Girka // Вопросы атомной науки и техники. — 2014. — № 6. — С. 37-40. — Бібліогр.: 7 назв. — англ.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
_version_ 1859651251455655936
author Pavlenko, I.V.
Melnyk, D.A.
Prokaieva, A.O.
Girka, I.O.
author_facet Pavlenko, I.V.
Melnyk, D.A.
Prokaieva, A.O.
Girka, I.O.
citation_txt Precision of a FDTD method to simulate cold magnetized plasmas / I.V. Pavlenko, D.A. Melnyk, A.O. Prokaieva, I.O. Girka // Вопросы атомной науки и техники. — 2014. — № 6. — С. 37-40. — Бібліогр.: 7 назв. — англ.
collection DSpace DC
container_title Вопросы атомной науки и техники
description The finite difference time domain (FDTD) method is applied to describe the propagation of the transverse electromagnetic waves through the magnetized plasmas. The numerical dispersion relation is obtained in a cold plasma approximation. The accuracy of the numerical dispersion is calculated as a function of the frequency of the launched wave and time step of the numerical grid. It is shown that the numerical method does not reproduce the analytical results near the plasma resonances for any chosen value of time step if there is not a dissipation mechanism in the system. It means that FDTD method cannot be applied straightforward to simulate the problems where the plasma resonances play a key role (for example, the mode conversion problems). But the accuracy of the numerical scheme can be improved by introducing some artificial damping of the plasma currents. Although part of the wave power is lost in the system in this case but the numerical scheme describes the wave processes in an agreement with analytical predictions. Метод конечных разностей во временной области (КРВО) применяется для описания распространения поперечных электромагнитных волн через магнитоактивную плазму. Численное дисперсионное уравнение получено в приближении холодной плазмы. Точность численной дисперсии посчитана как функция частоты распространяющейся волны и величины временного шага численной сетки. Показано, что при отсутствии механизмов диссипации в системе численный метод не воспроизводит аналитические результаты возле плазменных резонансов для любого выбранного временного шага. Это значит, что метод КРВО нельзя применять прямолинейно для моделирования задач, в которых плазменные резонансы играют ключевую роль (например, в задачах по изучению конверсии мод). Но точность численной схемы можно улучшить через включение механизма затухания для плазменных токов. Хотя в этом случае часть энергии волны теряется в системе, но численная схема описывает волновые процессы в соответствии с аналитическими предсказаниями. Метод кінцевих різниць у часовій області (КРЧО) застосовується для вивчення поширення поперечних електромагнітних хвиль через магнітоактивну плазму. Числове дисперсійне рівняння отримано в наближенні холодної плазми. Точність числової дисперсії обрахована як функція частоти електромагнітної хвилі та величини часового кроку числової сітки. Показано, що за відсутністю механізмів дисипації в системі числовий метод не відтворює аналітичні результати біля плазмових резонансів для будь- якого вибраного часового кроку. Це означає, що метод КРЧО не можна застосовувати прямолінійно для моделювання задач, у яких плазмові резонанси відіграють ключову роль (наприклад,у задачах з вивчення конверсії мод). Але точність числової схеми можна покращити через включення механізму загасання плазмових струмів. Хоча в цьому випадку частина енергії хвилі втрачається в системі, але числова схема описує хвильові процеси у відповідності до аналітичних передбачень.
first_indexed 2025-12-07T13:33:33Z
format Article
fulltext ISSN 1562-6016. ВАНТ. 2014. №6(94) PROBLEMS OF ATOMIC SCIENCE AND TECHNOLOGY. 2014, №6. Series: Plasma Physics (20), p. 37-40. 37 PRECISION OF A FDTD METHOD TO SIMULATE COLD MAGNETIZED PLASMAS I.V. Pavlenko, D.A. Melnyk, A.O. Prokaieva, I.O. Girka V.N. Karazin Kharkiv National University, Kharkiv, Ukraine The finite difference time domain (FDTD) method is applied to describe the propagation of the transverse electromagnetic waves through the magnetized plasmas. The numerical dispersion relation is obtained in a cold plasma approximation. The accuracy of the numerical dispersion is calculated as a function of the frequency of the launched wave and time step of the numerical grid. It is shown that the numerical method does not reproduce the analytical results near the plasma resonances for any chosen value of time step if there is not a dissipation mechanism in the system. It means that FDTD method cannot be applied straightforward to simulate the problems where the plasma resonances play a key role (for example, the mode conversion problems). But the accuracy of the numerical scheme can be improved by introducing some artificial damping of the plasma currents. Although part of the wave power is lost in the system in this case but the numerical scheme describes the wave processes in an agreement with analytical predictions. PACS: 52.65.-y, 52.25.Xz INTRODUCTION The FDTD method [1,2] is widely used to simulate a propagation of the electromagnetic waves through the materials with various dielectric and magnetic properties which may vary in time and space. In classical case it solves the Maxwell’s equations in partial differential form. Particularly it is frequently applied in an antenna theory to design new antennas and in a waveguide theory to optimize energy transfer. Sometime the conductive properties of the materials depend on the external electromagnetic impact. In this case the Maxwell’s equations are supplemented by the equations which describe an interaction of the material media with the electromagnetic waves. For example, if medium contains the free charges it will be the motion equations for the free charges in the electromagnetic field. At present the FDTD method is used in plasma physics mainly to solve the Maxwell’s equations together with the plasma current equations for all plasma species [3]. It does not allow to reproduce the hot plasma effects which require solving the kinetic equation but gives a space- time distribution of the wave electromagnetic field in cold plasma approach. This approach can be used as a simplified approximation to calculate wave power transmission, reflection and damping (particularly for the problems of antenna- plasma coupling). Also application of the FDTD method in plasma physics can give some advantages to study a wave-particle interaction including the nonlinear processes but it requires a development of this method for the anisotropic and nonuniform media. It was shown [4, 5] that FDTD method can be used to describe the fast wave propagation through the cold magnetized plasmas including some phenomena of the mode conversion. But description of the plasma resonances which could appear along wave propagation through a nonuniform plasma by the FDTD method is under a question. The problem is that a wavelength tends to zero near the plasma resonances. As a result a fixed space step of the numerical grid does not allow to resolve correctly the small scale wave processes. Just the mode conversion problem operates with plasma resonances and cutoffs. Therefore the additional analysis of the FDTD scheme accuracy in simulations of the magnetized plasmas is required. 1. SET OF EQUATIONS AND DISCRETIZATION SCHEME The cold magnetized plasma is described by the set of Maxwell’s equations coupled by the plasma current equations for plasma species: , , 4 41 , 1 2       JEJJ t JJ ct E c B t B c E p (1) here J  describes an external current into plasma, is the collisional frequency of the plasma species, p and  are plasma and cyclotron frequencies respectively. A standard Yee cell algorithm [1] is applied to solve the set of equations. It means that Faraday equation and Ampere equation (Maxwell-Ampere equation without a displacement current) are discretized in the way that electric E  and magnetic B  fields are calculated at different space points shifted from each other by a half of the space numerical grid step. Moreover, the different components of the fields are calculated at different space points. It allows to present the curl numerically in a convenient form (the sense of Yee cell application). For convenient treatment of the boundary conditions a discretization model has been chosen with the components of electric field at the edges of the Yee cell and the components of the magnetic field at the faces of the Yee cell. 38 ISSN 1562-6016. ВАНТ. 2014. №6(94) In such a way the components of the electromagnetic field will be calculated at following space points: ),2/1,,(),,2/1,(),,,2/1( lkiElkiElkiE zyx .), 2 1 ,, 2 1 (), 2 1 , 2 1 ,( etclkiBlkiB yx The time discretization is built also in the way that the electric E  and magnetic B  fields are calculated at different time moments shifted from each other by a half of the time numerical grid step. The current equation requires to know all components of the current and the electric field at the same space points. We have chosen the space points to calculate the currents as a center of Yee cell. In such a way all components of the current density will be calculated at the same space points: .), 2 1 , 2 1 , 2 1 (), 2 1 , 2 1 , 2 1 ( etclkiJlkiJ yx As a result the current data from the centers will be averaged to the edges to calculate the electric field from the Maxwell- Ampere equation. And the electric field data from the edges will be averaged to the centers to calculate the currents from third equation of the set (1). The FDTD scheme has been built with E  - J  collocation in time. Then the Maxwell- Faraday equation is explicit: )( )2/1,2/1(1 n nn E t B c   . (2) The mixed explicit/implicit scheme has been chosen for the Maxwell- Ampere and current equations: .][ 4 ,)( 41 ),1(),1( 2 ),1(),1( ),1()2/1( ),1( nnnnpnnnn nnn nn JEJJ t JJ c B t E c    (3) The set of the Maxwell’s equations with the current equations (1) has been discretized for a harmonic signal: all components of the vector values are proportional to ))(exp( ntmzklykixkj zyx . In our consideration the z axis is directed along the external magnetic field. Direct calculations issue the discretized dispersion relation: ,0)( ))()( )()(( )( ~ 2 21134 4 222 2 22 11 22 31 22 312 2 2 3 2 1 2 1 2 yx yxyx yx zy y zx x zy y x x c b dddd dddd c b dddk (4) where 2/ )2/sin( t t b is the discretized wave frequency, 2/ )2/sin( x xk d x x , 2/ )2/sin( y yk d y y , 2/ )2/sin( z zk d z z are the components of the discretized wave vector k ~ and 2222~ zyx dddk , c is speed of light. The components of the discretized electromagnetic tensor can be presented in convenient form: , ) ~~ ( ~~ 1 22 222 1 z pzyx bb baa (5) , ) ~~ ( ~~ 1 22 222 1 z pzxy bb baa (6) , ) ~~ ( ~~ 22 22 2 z zpzyx bb aaa (7) ,~ ~ 1 222 3 bb aa pyx (8) where the coefficients )2/cos( xka xx , )2/cos( yka yy , )2/cos( zka zz appear due to averaging the electric field values from the edges to the centers in the current equations and averaging the current density values from the centers to the edges in the Maxwell- Ampere equation, App ~ is “numerical” plasma frequency, Azz ~ is “numerical” cyclotron frequency, Ajbb ~ and )2/cos( tA is a multiplier coming from a level of the explicit/implicit presentation. Let’s consider the transverse propagation of the electromagnetic waves through the magnetized plasmas. The transverse propagation means that 0zk . Moreover, to resolve the shortest wavelengths of the problem the steps of the space grid should be enough small to provide xx kd , yy kd , zz kd , 1zyx aaa which mean 111 yx . Finally the general dispersion relation (4) is reduced to well known analytical form [6] but, please, pay an attention to the “numerical” definitions of the usual notations: ,0)())(( 2 2 2 134 4 2 2 2 131 2 2 2 4 1 c b k c b k (9) which issues the ordinary 2 2 3 2 c b k and the extraordinary 1 2 2 2 1 2 2 2 )( c b k waves. 2. ACCURACY OF THE NUMERICAL DISPERSION FOR THE EXTRAORDINARY WAVE The numerical dispersion of the studied FDTD scheme is compared with the analytical dispersion according to the motivation in [7]. The accuracy of the numerical scheme is defined as: )Re( )Re()Re( 2 22 an numan n nn , (10) where 2 ann is square of the analytical refraction index and 2 numn is square of the numerical refraction index (9) for the extraordinary wave from. The hydrogen homogeneous plasma with typical fusion parameters has been chosen for the accuracy tests. The densities of the electrons and ions are equal to 3·10 13 cm -3 , the external magnetic field value is 3.4 T. The accuracy is built in Fig. 1. Although the requested accuracy can be achieved in the wide frequency range but it is not achievable near the lower hybrid frequency. ISSN 1562-6016. ВАНТ. 2014. №6(94) 39 This region is zoomed in Fig. 2. As a sequence the FDTD scheme could issue a wrong result for the waves with the frequencies which are close to the lower hybrid one. And it is correct for any value of time step regardless of the selected (enough small to resolve the shortest wavelength) space step. Of course, there is a condition of the numerical stability for the studied scheme which provides the Fig. 1. The accuracy (10) of the numerical dispersion for the extraordinary wave in the low frequency range as a function of the normalized wave frequency and the normalized time step. Please note that for the convenience the large accuracy values are shown as 1 (absolutely wrong result) Fig. 2. The same as in Fig. 1 but with zoomed area near the lower hybrid frequency. The red area reaches the horizontal axis: numerical dispersion does not reproduce the analytical result for any chosen time step acceptable relation between the space and time steps of the numerical grid. This condition could be even stronger than the condition on the numerical accuracy but it was not strictly defined yet for the magnetized plasmas. Also, although only the lower hybrid resonance has been considered here but the similar conclusion can be made for other plasma resonances: the nondissipative FDTD scheme cannot treat correctly the wave processes near the plasma resonances. But it does not mean that the FDTD scheme can not be applied at all to simulate the processes near the plasma resonances (particularly near the lower hybrid resonance). A way to solve this problem is introducing some artificial damping of the plasma currents. In our case the collisional frequency of the plasma particles has been introduced to the current equations. The changes in the numerical accuracy can be seen from a comparison of the Fig. 2 and Fig. 3. Although the accuracy near the lower hybrid frequency depends still strongly on the time step value but there are time step values which keep the numerical accuracy in the selected limits. Fig. 3. The same as in Fig. 2 but with the artificial damping provided by the collisional frequency which is equal to 0.01 part of the wave frequency. The red area does not reach the horizontal axis: numerical dispersion reproduces the analytical result for some range of the time step values with the acceptable accuracy CONCLUSIONS The numerical dispersion of the transverse extraordinary wave for one of the widely used FDTD scheme is compared with the analytical dispersion of this wave in the low frequency range. The accuracy of the numerical dispersion is calculated as function of the wave frequency and the time step of numerical grid. It is shown that numerical dispersion issues wrong result in the frequency range near the lower hybrid frequency for any value of time step even if the space step is enough small for well resolution of the shortwavelength perturbations. Intuitively this result follows from vacuum Courant condition [2] of the numerical stability. But numerical accuracy and numerical stability are important but different characteristics of the FDTD numerical scheme. In general one of them does not follow from another. Moreover there is not a proof for Courant condition of the magnetized plasma: many authors suppose and motivate that vacuum condition is still valid for the magnetized plasmas but the conclusion is based on statistical data only. Here we have proved that the accuracy of scheme without current dissipation is not enough to simulate 40 ISSN 1562-6016. ВАНТ. 2014. №6(94) waves near the plasma resonances regardless of the numerical stability. The accuracy of the FDTD scheme near the plasma resonances can be improved by introducing some artificial dissipation (for example, the artificial damping through the collisions). In this case there is a possibility to choose the time step value to reach the selected numerical accuracy but the price of this possibility is a nonconservative system. Amount of power which is “lost” in the system due to this artificial damping is the subject of separate discussion, especially for the scenarios with the mode conversion processes. So, the artificial collisions allow to treat the plasma resonances by the FDTD method reproducing the space- time distribution of the electromagnetic field of the different plasma modes but the power conservation low for these modes should be corrected on the “lost” power. REFERENCES 1. K. Yee. A numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media // IEEE Transactions on Antennas and Propagation. 1966, v. 14, № 3, p. 02-307. 2. A. Taflove, S.C. Hagness. Computational Electromagnetics: Finite-Difference Time-Domain Method. Norwood: “Artech House”, 2005, p. 852. 3. Y. Yu, J.J. Simpson. An E-J collocated 3-D FDTD model of electromagnetic wave propagation in magnetized cold plasma // IEEE Transactions on Antennas and Propagation. 2010, v. 58, № 2, p. 469-478. 4. D.N. Smithe. Finite-difference time-domain simulation of fusion plasmas at radiofrequency time scales // Physics of Plasmas. 2007, v. 14, p. 056104. 5. W. Tierens, D.De Zutter. An unconditionally stable time- domain discretization on Cartesian meshes for the simulation of nonuniform magnetized plasma // Journal of Computational Physics. 2012, v. 231, p. 5144-5166. 6. A.I. Akhiezer, I.A. Akhiezer, R.V. Polovin, A.G. Sitenko, K.N. Stepanov. Plasma electrodynamics. v. 1:Linear theory. Oxford: “Pergamon Press”, 1975, p. 428. 7. S.A. Cummer. An analysis of new and existing FDTD methods for isotropic cold plasma and a method for improving their accuracy // IEEE Transactions on Antennas and Propagation. 1997, v. 45, № 3, p. 392-400. Article received 18.10.2014 ТОЧНОСТЬ МЕТОДА КОНЕЧНЫХ РАЗНОСТЕЙ ВО ВРЕМЕННОЙ ОБЛАСТИ ПРИ МОДЕЛИРОВАНИИ ХОЛОДНОЙ МАГНИТОАКТИВНОЙ ПЛАЗМЫ И.В. Павленко, Д.А. Мельник, А.О. Прокаева, И.А. Гирка Метод конечных разностей во временной области (КРВО) применяется для описания распространения поперечных электромагнитных волн через магнитоактивную плазму. Численное дисперсионное уравнение получено в приближении холодной плазмы. Точность численной дисперсии посчитана как функция частоты распространяющейся волны и величины временного шага численной сетки. Показано, что при отсутствии механизмов диссипации в системе численный метод не воспроизводит аналитические результаты возле плазменных резонансов для любого выбранного временного шага. Это значит, что метод КРВО нельзя применять прямолинейно для моделирования задач, в которых плазменные резонансы играют ключевую роль (например, в задачах по изучению конверсии мод). Но точность численной схемы можно улучшить через включение механизма затухания для плазменных токов. Хотя в этом случае часть энергии волны теряется в системе, но численная схема описывает волновые процессы в соответствии с аналитическими предсказаниями. ТОЧНІСТЬ МЕТОДУ КІНЦЕВИХ РІЗНИЦЬ У ЧАСОВІЙ ОБЛАСТІ ПРИ МОДЕЛЮВАННІ ХОЛОДНОЇ МАГНІТОАКТИВНОЇ ПЛАЗМИ І.В. Павленко, Д.О. Мельник, А.О. Прокаєва, І.О. Гірка Метод кінцевих різниць у часовій області (КРЧО) застосовується для вивчення поширення поперечних електромагнітних хвиль через магнітоактивну плазму. Числове дисперсійне рівняння отримано в наближенні холодної плазми. Точність числової дисперсії обрахована як функція частоти електромагнітної хвилі та величини часового кроку числової сітки. Показано, що за відсутністю механізмів дисипації в системі числовий метод не відтворює аналітичні результати біля плазмових резонансів для будь- якого вибраного часового кроку. Це означає, що метод КРЧО не можна застосовувати прямолінійно для моделювання задач, у яких плазмові резонанси відіграють ключову роль (наприклад,у задачах з вивчення конверсії мод). Але точність числової схеми можна покращити через включення механізму загасання плазмових струмів. Хоча в цьому випадку частина енергії хвилі втрачається в системі, але числова схема описує хвильові процеси у відповідності до аналітичних передбачень.
id nasplib_isofts_kiev_ua-123456789-81191
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
issn 1562-6016
language English
last_indexed 2025-12-07T13:33:33Z
publishDate 2014
publisher Національний науковий центр «Харківський фізико-технічний інститут» НАН України
record_format dspace
spelling Pavlenko, I.V.
Melnyk, D.A.
Prokaieva, A.O.
Girka, I.O.
2015-05-13T15:11:34Z
2015-05-13T15:11:34Z
2014
Precision of a FDTD method to simulate cold magnetized plasmas / I.V. Pavlenko, D.A. Melnyk, A.O. Prokaieva, I.O. Girka // Вопросы атомной науки и техники. — 2014. — № 6. — С. 37-40. — Бібліогр.: 7 назв. — англ.
1562-6016
PACS: 52.65.-y, 52.25.Xz
https://nasplib.isofts.kiev.ua/handle/123456789/81191
The finite difference time domain (FDTD) method is applied to describe the propagation of the transverse electromagnetic waves through the magnetized plasmas. The numerical dispersion relation is obtained in a cold plasma approximation. The accuracy of the numerical dispersion is calculated as a function of the frequency of the launched wave and time step of the numerical grid. It is shown that the numerical method does not reproduce the analytical results near the plasma resonances for any chosen value of time step if there is not a dissipation mechanism in the system. It means that FDTD method cannot be applied straightforward to simulate the problems where the plasma resonances play a key role (for example, the mode conversion problems). But the accuracy of the numerical scheme can be improved by introducing some artificial damping of the plasma currents. Although part of the wave power is lost in the system in this case but the numerical scheme describes the wave processes in an agreement with analytical predictions.
Метод конечных разностей во временной области (КРВО) применяется для описания распространения поперечных электромагнитных волн через магнитоактивную плазму. Численное дисперсионное уравнение получено в приближении холодной плазмы. Точность численной дисперсии посчитана как функция частоты распространяющейся волны и величины временного шага численной сетки. Показано, что при отсутствии механизмов диссипации в системе численный метод не воспроизводит аналитические результаты возле плазменных резонансов для любого выбранного временного шага. Это значит, что метод КРВО нельзя применять прямолинейно для моделирования задач, в которых плазменные резонансы играют ключевую роль (например, в задачах по изучению конверсии мод). Но точность численной схемы можно улучшить через включение механизма затухания для плазменных токов. Хотя в этом случае часть энергии волны теряется в системе, но численная схема описывает волновые процессы в соответствии с аналитическими предсказаниями.
Метод кінцевих різниць у часовій області (КРЧО) застосовується для вивчення поширення поперечних електромагнітних хвиль через магнітоактивну плазму. Числове дисперсійне рівняння отримано в наближенні холодної плазми. Точність числової дисперсії обрахована як функція частоти електромагнітної хвилі та величини часового кроку числової сітки. Показано, що за відсутністю механізмів дисипації в системі числовий метод не відтворює аналітичні результати біля плазмових резонансів для будь- якого вибраного часового кроку. Це означає, що метод КРЧО не можна застосовувати прямолінійно для моделювання задач, у яких плазмові резонанси відіграють ключову роль (наприклад,у задачах з вивчення конверсії мод). Але точність числової схеми можна покращити через включення механізму загасання плазмових струмів. Хоча в цьому випадку частина енергії хвилі втрачається в системі, але числова схема описує хвильові процеси у відповідності до аналітичних передбачень.
en
Національний науковий центр «Харківський фізико-технічний інститут» НАН України
Вопросы атомной науки и техники
Нагрев плазмы и поддержание тока
Precision of a FDTD method to simulate cold magnetized plasmas
Точность метода конечных разностей во временной области при моделировании холодной магнитоактивной плазмы
Точність методу кінцевих різниць у часовій області при моделюванні холодної магнітоактивної плазми
Article
published earlier
spellingShingle Precision of a FDTD method to simulate cold magnetized plasmas
Pavlenko, I.V.
Melnyk, D.A.
Prokaieva, A.O.
Girka, I.O.
Нагрев плазмы и поддержание тока
title Precision of a FDTD method to simulate cold magnetized plasmas
title_alt Точность метода конечных разностей во временной области при моделировании холодной магнитоактивной плазмы
Точність методу кінцевих різниць у часовій області при моделюванні холодної магнітоактивної плазми
title_full Precision of a FDTD method to simulate cold magnetized plasmas
title_fullStr Precision of a FDTD method to simulate cold magnetized plasmas
title_full_unstemmed Precision of a FDTD method to simulate cold magnetized plasmas
title_short Precision of a FDTD method to simulate cold magnetized plasmas
title_sort precision of a fdtd method to simulate cold magnetized plasmas
topic Нагрев плазмы и поддержание тока
topic_facet Нагрев плазмы и поддержание тока
url https://nasplib.isofts.kiev.ua/handle/123456789/81191
work_keys_str_mv AT pavlenkoiv precisionofafdtdmethodtosimulatecoldmagnetizedplasmas
AT melnykda precisionofafdtdmethodtosimulatecoldmagnetizedplasmas
AT prokaievaao precisionofafdtdmethodtosimulatecoldmagnetizedplasmas
AT girkaio precisionofafdtdmethodtosimulatecoldmagnetizedplasmas
AT pavlenkoiv točnostʹmetodakonečnyhraznosteivovremennoioblastiprimodelirovaniiholodnoimagnitoaktivnoiplazmy
AT melnykda točnostʹmetodakonečnyhraznosteivovremennoioblastiprimodelirovaniiholodnoimagnitoaktivnoiplazmy
AT prokaievaao točnostʹmetodakonečnyhraznosteivovremennoioblastiprimodelirovaniiholodnoimagnitoaktivnoiplazmy
AT girkaio točnostʹmetodakonečnyhraznosteivovremennoioblastiprimodelirovaniiholodnoimagnitoaktivnoiplazmy
AT pavlenkoiv točnístʹmetodukíncevihríznicʹučasovíioblastíprimodelûvanníholodnoímagnítoaktivnoíplazmi
AT melnykda točnístʹmetodukíncevihríznicʹučasovíioblastíprimodelûvanníholodnoímagnítoaktivnoíplazmi
AT prokaievaao točnístʹmetodukíncevihríznicʹučasovíioblastíprimodelûvanníholodnoímagnítoaktivnoíplazmi
AT girkaio točnístʹmetodukíncevihríznicʹučasovíioblastíprimodelûvanníholodnoímagnítoaktivnoíplazmi